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Thejtimely  design  and  development  of  the  next  generation  of  tilt- 
rotor  aircraft  4*fVXj/>depend  heavily  on  ih^in-depth  understanding  of 
existing  XV-15  dynamics  and  the  availability  of  fully  validated  simula¬ 
tion  models.  Previous  studies  have  considered  aircraft  and  simulation 
trim  characteristics,  but  analyses  of  basic  flight  vehicle  dynamics  have 
been  limited  to  qualitative  pilot  evaluations.  /'  The  present  study  has 
the  following  objectives: 

r.  Documentation  and  evaluation  of  XV-15  bare-airframe  dynamics/ 

2.  Comparison  of  aircraft  and  simulation  responses;  a-r’a 

3.  Development  of  a  validated  transferjfunction  description  of  the 
XV-15  needed  for  future  studies 


A  nonparametric  frequency-response  approach  is  used  which  does  not 
depend  on  assumed  model  order  or  structure.  Transferjfunction  represen¬ 


tations  are  subsequently  derived  which  fit  the  frequency  responses  in 
the  bandwidth  of  greatest  concern  for  piloted  handling- qualities  and 
control-system  applications. 


This  study  involved  the  planning  and  execution  of  flight  tests  on 
the  XV-15  aircraft  and  piloted-simulation  for  four  flight  conditions 
from  hover  to  cruise.  Improved  test  techniques  and  pilot-training 
procedures  were  devised.  Analytical  software  tools  were  developed  (or 
adapted)  which  allow  the  identification  of  high-resolution  spectral 
responses  and  the  derivation  and  validation  of  multi-input/multi-output 


iv 


'  transfer ^function  models.  These  technics  were  applied  in  an  extensive 
evaluation  of  the  open-loop  flight /dynamics  of  the  XV- 15  aircraft  and 
simulation  mathematical  models.  ^Deficiencies  in  the  mathematical  models 
were  exposed  and  documented.  Finally,  a  new,  fully  validated  transfer- 
function  model  was  derived  for  the  hover  and  oruise  flight  conditions. 
The  methods  developed  in  this  study  have  subsequently  been  applied  in  a 
number  of  other  flight-test  programs  and  have  been  included  in  the  U.S. 
Army's  updated  helicopter  hand ling -qualities  specification. 

Approved  for  publication: 

By 
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Dean  of  Graduate  Studies 
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Chapter  1 
INTRODUCTION 


1 . 1  Background  and  Objectives 

The  tilt-rotor  concept  combines  the  hovering  advantages  of  the 
helicopter  with  the  cruise  advantages  of  a  fixed-wing  aircraft. 
Rotor/engine  nacelles  at  the  wing  tips  are  rotated  to  the  vertical 
position  for  hbvering  flight  and  to  the  horizontal  position  for  cruising 
flight.  The  XV-15  research  aircraft  (Fig.  1.1)  was  jointly  developed  by 
the  U.S.  Army,  NASA,  and  the  Navy  to  demonstrate  tilt-rotor  technol¬ 
ogy.  A  key  objective  of  this  project  was  to  achieve  good  piloted 
handling-qualities  characteristics  in  hovering  flight  by  using  an 
advanced  stability  and  control  augmentation  system  (SCAS).  This  objec¬ 
tive  was  emphasized  because  of  serious  handling-qualities  deficiencies 
in  hovering  flight  which  were  encountered  with  the  original  tilt-rotor 
demonstrator — the  XV-3  (Ref.  1).  Two  XV-15  aircraft  were  developed 
under  contract  to  Bell  Helicopter  Company  and  delivered  to  Ames  Research 
Center  (ARC)  in  1980.  One  aircraft  (N703)  was  retained  at  Ames  for 
research  and  development  testing;  the  other  (N702)  was  leased  back  to 
the  contractor  for  operational  testing. 

Comprehensive  real-time  and  nonreal-time  simulation  codes  (Refs.  2, 
3)  were  developed  to  support  the  design  and  testing  of  the  XV-15. 
Moving-base  simulation  facilities  at  ARC  were  extensively  used  for  pilot 
training  before  the  first  flight  tests,  and  subsequently  for  advanced 
automatic  flight  control  system  (AFCS)  development.  The  XV-15  simula¬ 
tion  code  covers  the  entire  operating  envelope,  with  a  full  nonlinear 


Fig.  1.1.  The  XV-15  Tilt-Rotor  Aircraft,  (a)  Hover  Configuration; 
(b)  Cruise  Configuration. 

representation  of  the  aircraft.  Wing/body/tail  aerodynamics  are  deter¬ 
mined  from  extensive  look-up  tables  of  full-scale  wind-tunnel  data 
obtained  in  the  NASA  40-  by  80-Foot  Wind  Tunnel  facility.  Rotor  cal¬ 
culations  assume  quasi-steady  flapping  and  are  based  on  modified  Bailey 
equations  with  uniform  rotor  inflow.  The  aerodynamic  interactions 
between  the  two  rotors,  and  the  rotor  interference  with  the  other 
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aircraft  elements  are  modeled  in  detail.  Also  modeled  are  numerous  sub¬ 
system  dynamics  such  as  the  engine  drive  train  and  governor.  The  XV-15 
mathematical  model  was  the  most  complex  ever  developed  to  that  time  for 
real-time  piloted  simulation  at  ARC  (Ref.  4).  The  nonreal-time  version 
is  routinely  used  to  support  control-system  development  and  flight-test 
planning.  The  XV-15  simulation  mathematical  models  have  been  exten¬ 
sively  correlated  with  static  trim  and  performance  flight  data;  the 
comparison  is  generally  excellent  (Ref.  5).  However,  dynamic  checks 
have  been  very  limited  (Ref.  4),  with  most  of  the  validation  in  this 
area  centered  on  pilot  subjective  comparison  of  the  aircraft  and 
motion-based  simulator  response. 

The  author,  as  a  staff  member  of  an  Army  research  team  responsible 
for  simulation  technology,  started  in  1983  to  conduct  a  comprehensive 
study  to  validate  the  open-loop  dynamic  response  fidelity  of  the 
piloted-simulation  mathematical  model.  An  in-depth  understanding  of 
XV-15  dynamics  and  the  availability  of  fully  validated  simulation  models 
were  considered  important  for  the  timely  design  and  development  of  the 
Joint  Services  Operational  Tilt-Rotor  Aircraft— the  JVX,  now  designated 
the  V-22.  To  fulfill  these  needs,  the  study  was  initiated  with  the 
following  three  major  objectives: 

1.  Document  the  open-loop  dynamic  characteristics  of  the  XV-15 
aircraft  from  flight  tests  for  several  operating  conditions  including 
hover 

2.  Compare  aircraft  and  simulation  response  characteristics  to 
identify  problem  areas  in  the  mathematical  modeling 

3.  Develop  a  validated  transfer-function  model  description  of  the 
XV-15  needed  for  future  studies 


Emphasis  was  initially  placed  on  the  hover  flight  condition,  where 
unstable  open-loop  dynamics  lead  to  the  most  critical  handling-qualities 


problems . 

A  key  consideration  in  planning  this  study  was  the  selection  of  an 
appropriate  dynamics  identification  method. 


1.2  Dynamics  Identification  Methods 

Dynamics  identification  methods  generally  fall  into  two  cate- 

* 

gories:  frequency  domain  and  time  domain.  Each  approach  has  its  inher¬ 
ent  strengths  and  weaknesses  which  make  it  best  suited  for  particular 
applications. 

In  time-domain  (maximum-likelihood)  identification  (Fig.  1.2),  the 
aircraft  dynamics  are  modeled  by  a  set  of  differential  equations 


Fig.  1.2.  Time-Domain  Identification  Method. 
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describing  the  external  forces  and  moments  in  terms  of  state  and  control 
variables.  The  unknown  coefficients  in  the  equations  are  the  stability 
derivatives,  which  are  identified  by  least-squares  fitting  of  the  mea¬ 
sured  time-responses  (output-error  method).  Such  an  approach  allows  a 
direct  comparison  of  stability  derivatives  obtained  in  the  wind  tunnel 
and  those  of  the  actual  flight  vehicle.  Transfer  functions  and  fre¬ 
quency  responses  may  be  calculated  from  the  state-space  model.  A  key 
aspect  of  time-domain  identification  is  that  an  a  priori  model  formula¬ 
tion  must  be  assumed.  This  important  step  involves  consideration  of 
model  structure,  order,  and  important  nonlinearities.  Such  information 
is  generally  not  well  known  on  a  new  vehicle  such  as  the  XV-15,  and 
incorrect  model  formulation  can  bias  the  parameter  estimates  (Ref.  6). 
Also,  models  which  provide  a  good  fit  in  the  time-domain  do  not  neces¬ 
sarily  yield  accurate  transfer  functions,  since  time-domain  identifica¬ 
tion  techniques  weight  their  results  more  heavily  at  low  frequency  where 
most  of  the  data  points  are  concentrated. 

*'  The  frequency-domain  identification  approach  shown  in  Fig.  1.3  uses 
spectral  analysis  methods  to  extract  the  frequency  responses  between 
selected  input  and  output  pairs.  The  identification  results  are  usually 
presented  in  Bode-plot  format,  that  is,  log -magnitude  and  phase  of  the 
input-to-output  versus  log-frequency.  These  identification  results  are 
nonparametric  because  no  model  structure  is  assumed.  As  such,  they  are 
useful  for  flight-control  system  design  and  pilot-in-the-loop  handling- 
qualities  studies.  Frequency  responses  obtained  from  real-time  and 
nonreal-time  simulations  can  be  compared  directly  with  the  flight  data 
to  expose  limitations  and  discrepancies  in  the  simulator  models.  The 
fact  that  this  comparison  can  be  made  initally  without  an  a  priori 


Fig.  1.3.  Frequency-Response  Identification  Method. 


assumption  of  model  structure  or  order  is  especially  important  for 
verifying  mathematical  models  of  new  aircraft  configurations.  When  the 
model  structure  and  parametric  values  are  required,  they  may  be  obtained 
by  fitting  the  frequency-responses  with  transfer-function  models  to 
extract  modal  characteristics.  Examples  of  this  application  are  the 
testing  of  handling-qualities  specifications  given  in  lower-order  equiv¬ 
alent  system  terms,  and  the  examination  of  transfer  function-based 
control  system  designs.  Since  this  fitting  procedure  is  completed  after 
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the  frequency  response  is  extracted,  the  order  of  the  transfer  function 
can  be  selected  to  avoid  an  overparameterized  model.  Multi-input/multi- 
output  (MIMO)  frequency-response  methods  are  suitable  for  extracting  a 
transfer  matrix  which  includes  the  important  coupling  effects.  Finally, 
the  extracted  models  are  driven  with  the  flight  data  to  verify  the 
time-domain  characteristics.  Models  identified  by  frequency-domain 
techniques  are  often  most  accurate  at  mid-  and  high-frequency  (initial 
time-history  transients),  which  is  the  region  of  greatest  concern  to  the 
pilot.  The  low-frequency  and  steady-state  response  prediction  of  the 
extracted  models  is  generally  not  as  good  as  in  the  time-domain  identi¬ 
fication  approach. 

Since  the  completion  of  the  preceding  objectives  depends  on  obtain¬ 
ing  an  accurate  characterization  of  the  input-to-output  transient  dynam¬ 
ics  and  piloted  handling-qualities  of  a  new  aircraft  configuration, 
rather  than  on  obtaining  a  stability  derivative  model  (necessary  for 
example  to  validate  the  wind-tunnel  data  base),  the  frequency -domain 
approach  is  the  natural  choice. 

1.3  A  Historical  Summary  of  Research  in  Frequency-Response 
Identification 

The  earliest  reported  research  in  frequency-response  identification 
of  aircraft  dynamics  from  flight-test  data  was  conducted  at  the  Cornell 
Aeronautical  Laboratory  beginning  in  1945  (summarized  in  Ref.  7). 
Steady-state  sine-wave  inputs  were  used  to  (laboriously)  extract  the 
frequency  responses  of  the  North  American  B-25J  (fixed-wing)  aircraft. 
Then,  lower-order  transfer-function  models  were  derived  from  a  least- 
squares  fit  of  the  frequency -responses  (displayed  on  a  polar  plot). 


Fourier  transform  methods  were  subsequently  developed  (Refs.  7,  8)  to 
allow  frequency-response  identification  from  (shorter-duration) 
discrete-maneuver  data,  such  as  that  obtained  from  step  and  pulse 
inputs.  These  techniques  were  applied  in  flight  research  activities  at 
the  Air  Force  Flight  Test  Center  (Edwards  Air  Force  Base)  during  the 
1950s  (see  Ref.  9  for  a  list  of  references).  As  pointed  out  in  Ref.  9, 
all  of  these  early  efforts  in  frequency-response  identification  suffered 
from  the  lack  of  large-scale  computing  power.  The  development  of  the 
fast  Fourier  transform  (FFT)  algorithms  in  the  1960s,  and  the  signifi¬ 
cantly  improved  computing  capabilities  of  this  period  led  to  much 
greater  interest  and  success  in  frequency-response  identification. 

A  comprehensive  facility  for  multivariable  frequency-response 
(matrix)  identification  and  analysis  (FRA)  based  on  the  FFT  was  devel¬ 
oped  by  Twisdale  (1975)  at  the  Edwards  Test  Center  (complete  documenta¬ 
tion  in  Ref.  10).  One  key  feature  in  this  identification  approach  was 
incorporation  of  the  ordinary,  partial,  and  multiple  coherence  function 
calculations  which  provide  an  important  measure  of  spectral  estimation 
accuracy  (Chap.  3).  Identification  was  achieved  from  flight  data  of 
tracking  and  refueling  handling-qualities  tests  instead  of  from  data 
obtained  with  prescribed  inputs. 

Marchand  and  Koehler  (Ref.  11),  of  the  Institute  for  Flight  Mechan¬ 
ics  (DFVLR)  in  the  Federal  Republic  of  Germany,  developed  a  method  for 
extracting  the  stability -derivative  matrix  from  identified  frequency- 
responses;  flight  data  were  obtained  with  prescribed  "multi-step" 
inputs— an  outgrowth  of  research  in  optimal  input  design  (see  Ref.  6). 
Frequency-response  identification  from  flight  data  obtained  with  the 
prescribed  "frequency-sweep"  input  was  pioneered  by  Systems  Technology, 


Inc.,  (Refs.  12,  13).  Advancements  in  lower-order  transfer-function 


modeling  were  made  by  Hodgkinson  et  al.  (Ref.  14),  Bischoff  and  Palmer 
(Ref.  15),  and  Mitchell  and  Hoh  (Ref.  16)  in  support  of  the  development 
of  an  updated  handling-qualities  specification  for  military  fixed-wing 
aircraft  (Ref.  17). 

The  author  was  the  first  to  extensively  identify  frequency- 
responses  and  transfer-functions  of  rotorcraft  from  flight  tests  using 
the  frequency-sweep  input  (Ref.  18).  The  frequency-domain  approach  used 
in  this  study  and  depicted  in  Fig.  1.3  draws  heavily  on  the  preceding 
researchers'  efforts.  The  frequency-response  identification  (FRESPID) 
program  developed  for  this  effort  (Chap.  3)  was  patterned  after  the  FRA 
facility  (Ref.  10).  A  major  advancement  in  the  FRESPID  program  is  the 
incorporation  of  the  chirp  z-transform,  an  algorithm  for  obtaining 
frequency-responses  which  are  of  much  higher  quality  than  those  obtain¬ 
able  from  standard  FFT  procedures.  Also,  the  author  has  stressed  the 
importance  of  deriving  physically  consistent  transfer-function  models 
and  verifying  those  models  with  time-domain  data  not  used  in  the  identi¬ 
fication  procedure. 

1.4  Scope  of  Research 

This  study  involved  the  planning  and  execution  of  flight  tests  on 
the  XV- 15  aircraft  and  piloted-simulation  for  four  flight  conditions 
from  hover  to  cruise.  Improved  test  techniques  and  associated  pilot 
training  procedures  were  devised.  Analytical  software  tools  were  devel¬ 
oped  (or  adapted)  which  allow  the  identification  of  high-resolution 
spectral  responses  and  the  derivation  and  validation  of  MIMO  transfer- 
function  models.  These  techniques  were  applied  in  an  extensive 


evaluation  of  the  open-loop  flight  dynamics  of  the  XV- 15  and  simulation 
mathematical  models.  Mathematical  model  deficiencies  were  clearly 
exposed  and  documented.  Finally,  a  new,  fully  validated  transfer- 
function  model  was  derived  for  the  hover  and  cruise  flight  conditions. 

1.5  Extensions  and  Applications 

As  in  most  research  efforts,  the  scope  of  this  project  grew  well 
beyond  the  original  objectives  outlined  in  Sec.  1.1.  The  advantages  of 
the  frequency-domain  approach  for  documenting  the  response  characteris¬ 
tics  of  new  configurations  became  readily  apparent  in  the  XV-15  study 
and  led  to  the  use  of  this  technique  in  a  variety  of  related  helicopter 
flight  test  projects. 

Frequency-domain  testing  of  the  Bell  21^-ST  single-rotor  helicopter 
was  completed  by  the  author  in  October  1985  to  support  the  Army's  devel¬ 
opment  of  an  updated  handling-qualities  specification  (Ref.  19).  A 
discussion  of  the  testing  procedures  and  identification  results  for  the 
Bell  21H-ST  is  presented  in  Ref.  20.  This  project  demonstrated  the 
feasibility  of  frequency-sweep  testing  and  frequency-domain  identifica¬ 
tion  as  a  specification-compliance  documentation  tool.  Further,  the 
results  showed  that  very  low  order  transfer-function  models  can  accu¬ 
rately  predict  the  large  motion  time-domain  behavior  of  single-rotor 
helicopters— even  in  hover. 

The  frequency -domain  identification  method  was  also  used  to  docu¬ 
ment  the  dynamic  characteristics  of  the  CH-97B  research  aircraft 
(Refs.  6,  21).  Current  research  by  the  author  and  Kaletka  (Ref.  22) 
under  a  memorandum-of-understanding  between  the  United  States  and  the 
Federal  Republic  of  Germany  is  concerned  with  comparing  frequency  and 


time-domain  identification  results  for  the  open-loop  XV- 15  dynamics  in 
hover.  Research  by  the  author  and  Acree  (Ref.  23)  concerns  identifica¬ 
tion  of  XV- 15  structural  dynamics  using  the  frequency -domain  approach. 
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1.6  Organization  of  Dissertation 


The  central  focus  of  this  dissertation  is  the  documentation  and 
analysis  of  the  dynamics  of  a  unique  vehicle — the  XV- 15.  Chapter  2 
describes  the  general  vehicle  configuration  and  important  subsystems, 
and  reviews  the  flight  test  conditions.  The  frequency -domain  method, 
briefly  overviewed  in  the  present  chapter,  is  discussed  in  detail  in 
Chapter  3.  Chapters  4  and  5  present  the  dynamics  identification  results 
for  the  hover  and  cruise  flight  conditions.  Frequency  responses  <•  :d 
transfer-function  models  extracted  from  flight  tests  are  compared  with 
real-time  and  nonreal-time  simulation  results.  Time-history  verifica¬ 
tion  responses  show  the  ability  of  the  extracted  lower-order  transfer- 
function  models  to  accurately  predict  large-amplitude  response  charac¬ 
teristics  of  the  XV-15  in  hover  and  cruise.  Chapter  6  reviews  the  key 
conclusions  and  contributions  of  this  work.  Recommendations  for  future 
research  in  frequency-response  identification  are  discussed  in  Chap¬ 
ter  7.  Appendix  A  reviews  the  important  advantages  of  the  chirp 
z- transform,  as  compared  to  the  standard  fast  Fourier  transform,  for 
frequency-response  identification  from  flight  data.  Appendix  B  dis¬ 
cusses  the  limitations  in  identifying  open-loop  vehicle  dynamics  from 
closed-loop  flight  tests.  An  appreciation  for  these  limitations  is 
important  in  applying  the  present  identification  method  to  future  flight 
tests.  The  text  of  this  dissertation  draws  heavily  on  the  author's 
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publications  concerning  the  tilt-rotor  study  (Refs.  18,  22,  24)  and 
related  research  efforts  (Refs.  6,  20). 


Chapter  2 

XV- 15  RESEARCH  AIRCRAFT  AND  FLIGHT  TEST  CONDITIONS 


An  understanding  of  XV- 15  aircraft  dynamics  requires  an  apprecia¬ 
tion  of  this  unique  vehicle's  configuration  and  key  subsystems.  This 
Chapter  first  reviews  the  vehicle  layout,  operational  flight  envelope, 
cockpit  controls,  automatic  flight-control  system,  and  research  instru¬ 
mentation,  and  then  outlines  the  flight-test  conditions  and  test  inputs 
for  this  study. 

2.1  Aircraft  Description 

2.1.1  General  Layout 

The  XV- 15  aircraft  is  a  lateral-tandem  rotor  vehicle.  The  rotors 
are  25  ft  in  diameter  and  the  spinner-to-spinner  span  is  32  ft 
(Fig.  2.1).  The  aircraft  is  powered  by  two  Lycoming  T-53  turboshaft 
engines  (1250  SHP  each),  one  mounted  in  each  wing-tip  nacelle.  These 
nacelles  rotate  from  i^  =  0  deg  (cruise  configuration)  to  iN  =  90  deg 
(hover  configuration).  Rearward  flight  acceleration  is  enhanced  by 
rotating  the  nacelles  to  ijj  =  95  deg. 

The  rotor  system  is  a  three-bladed  prop-rotor  with  a  stiff  in-plane 
gimbal  mounting  to  the  hub.  Rotor  tip-path-plane  orientation  is  con¬ 
trolled  through  standard  cyclic  and  collective  feathering  of  the  indi¬ 
vidual  blades.  The  resultant  hub  moments  cause  the  entire  rotor  system 
to  rotate  (flap)  as  a  unit,  rather  than  each  blade  independently  as  in 
articulated  helicopters.  Cross-shafting  between  the  nacelles  synchro¬ 
nizes  the  rotors  and  provides  a  one-engine-out  capability.  An  engine 


governor  system  adjusts  the  pilot's  power  lever  commands  to  the  collec¬ 
tive  pitch  of  the  blades  to  maintain  constant  rotor  rpm. 


2.1.2  Operating  Flight  Envelope 

The  operational  envelope  shown  in  Fig.  2.2  is  from  slow  rearward 
flight  to  a  maximum  true  airspeed  of  Vfc  =  300  knots  (with  the  nacelles 


NORMAL  RATED 
POWER  LIMIT 


Fig.  2.2.  Airspeed-Altitude  Envelope. 

locked  at  i^  =  0  deg).  When  the  nacelles  are  not  in  their  locked 
position,  the  maximum  airspeed  is  restricted  to  Vfc  =  170  knots  because 
of  aeroelastic  stability  limitations.  The  "conversion  corridor"  of 
Fig.  2.3  defines  the  allowable  range  of  nacelle-angle/airspeed 
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O  FLIGHT  TEST  CONDITIONS 
FOR  PRESENT  STUDY 


Fig.  2.3.  Conversion  Corridor. 

combinations  in  the  transition  flight  regime  between  hover  and  cruise. 
The  lower  boundary  on  nacelle  angle  is  determined  by  wing-loading 
limits,  and  the  upper  boundary  is  determined  by  blade  loading  limits. 

The  nominal  weight  of  the  XV- 15  is  13,000  lb  and  the  flight  endurance  is 
approximately  1  hr. 

2.1.3  Vehicle  Control 

The  pilot's  cockpit  controls  in  the  XV-15  are  the  same  as  those 
found  in  conventional  helicopters.  A  center  stick  is  used  to  control 
the  pitch  and  roll  motions.  Pedals  are  used  to  control  yaw  motions,  and 
a  power  lever  i3  used  to  control  vertical  motion  in  hover  and  airspeed 
in  forward  flight.  The  means  with  which  these  cockpit  controls  generate 
the  required  forces  and  moments  depends  on  the  aircraft  configuration. 
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The  control  method  for  helicopter-mode  flight  is  shown  in 
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Fig.  2.4.  Rolling  moments  are  generated  with  differential  collective, 
pitching  moments  with  uniform  longitudinal  cyclic  (fore-aft  tip-path- 
plane  rotation),  yawing  moments  with  differential  longitudinal  cyclic, 
and  heave  forces  with  uniform  collective.  In  cruising  flight,  the  rotor 
controls  are  phased  out  and  control  moments  are  generated  with  standard 
aircraft  aerodynamic  surfaces;  the  ailerons  produce  rolling  moments,  the 
elevator  produces  pitching  moments,  and  the  rudder  produces  yawing 
moments . 
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Fig.  2.4.  Control  Method  for  Helicopter-Mode  Flight. 


Speed  control  is  achieved  through  the  throttle/governor  system. 

The  rotor  control  (swashplate)  and  aerodynamic  surfaces  receive  commands 
from  the  cockpit  via  mechanical  linkages  and  hydraulic  actuators.  These 


commands  are  augmented  by  signals  from  the  automatic  flight-control 
system  through  series  actuators. 

The  phase-out  of  the  rotor  controls  from  hover  to  cruise  is  auto¬ 
matically  scheduled  as  a  function  of  nacelle  angle  i^.  The  aerodynamic 
surfaces  are  actuated  throughout  the  entire  flight  envelope,  although 
they  have  no  effect  at  very  low  airspeeds.  It  is  most  expedient  to 
refer  all  of  the  open-loop  vehicle  response  characteristics  (e.g., 
frequency  responses,  transfer  functions,  etc.)  to  these  surface  deflec¬ 
tions  since,  neglecting  the  small  servo  lags,  these  are  related  to  the 
sum  of  the  pilot  and  SCAS  inputs  (total  input  to  the  aircraft)  through  a 
mixing  ratio  which  is  constant  across  the  entire  flight  envelope. 

In  the  helicopter  configuration,  the  XV-15  open-loop  dynamics  are 
typical  of  hovering  vehicles,  although  the  number  of  cross-coupling 
paths  is  much  lower  than  in  single-rotor  helicopters.  The  planar  sym¬ 
metric,  lateral-tandem  configuration  yields  vehicle  dynamics  which  are 
essentially  decoupled  between  the  longitudinal  and  lateral/directional 
degrees-of- freedom.  One-way  coupling  from  roll  input  to  yaw  response  is 
significant.  This  is  due  to  the  differential  rotor  torque  which  accom¬ 
panies  the  differential  collective  inputs  for  roll  control.  The  open- 
loop  pitch  and  roll  dynamics  in  hover  are  characterized  by  highly 
unstable  low-frequency  attitude-speed  divergences.  The  heave  and  yaw 
dynamics  are  decoupled  from  the  attitude  motions  and  are  essentially 
first-order.  The  time-constants  of  these  dynamics  are  very  long 
(  =  10  sec)  because  the  tilt-rotor  configuration  has  a  relatively  high 
disk-loading  (low  heave  damping)  and  no  tail-rotor  (low  yaw  damping). 

Limited  evaluations  of  the  SCAS-off  handling-qualities  in  hover 
have  been  conducted  (Ref.  25).  These  evaluations  indicate  that  the 


severe  attitude  instability  and  small  translational  damping  result  in 
very  poor  handling-qualities.  Precision  hover  tasks  (out  of  ground- 
effect)  were  rated  with  Level  II  handling-qualities  ("adequate  perfor¬ 
mance  achieved  with  high  level  of  pilot  workload").  Precision  transla¬ 
tion  tasks  were  rated  with  Level  III  handling-qualities  ("unsatisfactory 
performance").  Quantitative  documentation  of  the  open-loop  XV-15  dynam¬ 
ics  in  hover  is  critical  to  "convert"  these  pilot-opinion  ratings  into 
engineering  requirements  for  future  tilt-rotor  configurations. 

2.1.4  Automatic  Flight  Control  System 

A  nominal  stability  and  control  augmentation  system  (SCAS)  was 
incorporated  in  the  XV-15  to  enhance  the  poor  bare-airframe  dynamic 
characteristics.  Improved  closed-loop  handling  qualities  were  achieved 
with  an  advanced  SCAS  developed  by  Churchill  and  Gerdes  (Ref.  25) 
between  1982  and  1984.  This  advanced  system  combines  feed-forward 
shaping,  model-following,  and  disturbance  rejection  to  achieve  crisp 
first-order  responses  in  the  attitude  rates  and  improved  inherent  sta¬ 
bility.  The  closed-loop  response  time  constants  are  about  0.4  sec  in 
pitch  and  roll  and  about  0.6  sec  in  yaw.  An  attitude-retention  system 
based  on  attitude-angle  feedback  provides  hands-off  stability  when  the 
controls  are  in  their  neutral  position.  In  forward  flight,  the  control 
system  provides  constant  load  factor  responses  to  longitudinal  stick 
inputs  and  coordinated  turn  responses  to  lateral  stick  inputs.  These 
improved  characteristics  give  the  aircraft  Level  I  handling  qualities 
("performance  satisfactory  without  improvement")  for  the  precision 
piloting  tasks  (Ref.  25). 

The  SCAS  is  a  two-channel  (summed)  fail/operate  system.  Thus, 
failure  of  either  channel  results  in  reversion  to  a  single  (lower-gain) 


channel.  Failure  of  the  second  channel  reverts  the  aircraft  back  to  the 


SCAS-off  configuration.  The  SCAS  authority  is  limited  to  20%  of  maximum 
stick  deflection  in  the  pitch  axis,  and  30%  of  the  maximum  control  in 
roll  and  yaw. 

2.1.5  Research  Instrumentation 

The  XV- 15  is  heavily  instrumented  to  provide  real-time  and  post¬ 
flight  engineering  data.  Measurements  of  over  150  variables  including 
angular  rates,  attitudes,  accelerations,  surface  deflections,  cockpit 
controls,  and  structural  loads  are  recorded  on  board  at  250  Hz  on 
pulse-code-modulation  (PCM)  tapes.  Besides  this  on-board  recorder,  a 
telemetry  (TM)  link  provides  safety-related  data  which  are  monitored  by 
a  six-member  flight-test  ground  crew.  This  extensive  instrumentation 
system  is  carefully  maintained  and  regularly  calibrated.  After  the 
flight  tests,  the  PCM  tapes  are  digitized,  and  the  data  are  converted  to 
engineering  units. 

The  three-axis  gyroscope  package  has  a  bandwidth  of  about  20  Hz, 
which  is  well  beyond  the  frequency  range  of  interest  for  rigid-body 
dynamics  identification.  Analog  measurements  of  controls  and  responses 
are  conditioned  with  50  Hz  analog  pre-filters  to  reduce  digital  alias¬ 
ing,  while  maintaining  a  broad  band  of  accurate  dynamic  measurements. 
Matching  the  filters  on  the  input  and  output  signals  minimizes  the  phase 
distortions  in  the  identified  frequency  responses. 

2.2  Flight  Test  Conditions 

The  following  four  flight  test  conditions  were  selected  to  span  the 
full  operating  range  of  the  XV-15  aircraft: 


1.  Hover:  Ambient  wind  less  than  5  knots,  iN  =  90  deg,  gear  down, 
out-of-ground  effect  (altitude  =  100  ft) 

2.  Low-speed  transition:  Indicated  airspeed  V^  =  100  knots 
(calibrated  airspeed,  VQ  =  105  knots),  iN  =  70  deg 

3.  High-speed  transition:  Indicated  airspeed  V.^  =  120  knots 
(calibrated  airspeed,  VQ  =  128  knots),  i^  =  30  deg 

4.  Cruise:  Indicated  airspeed  V^  =  170  knots  (calibrated  air¬ 
speed,  Vc  =  180  knots),  ijj  =  0  deg 

These  flight  conditions  are  noted  on  the  conversion  corridor  of 
Fig.  2.3. 

Dynamics  identification  tests  were  conducted  from  February  through 
December  1983,  in  about  12  flight  hours.  During  the  same  period, 
frequency-sweep  tests  were  conducted  in  the  Vertical  Motion  Simulator 
(VMS)  to  document  the  real-time  XV- 15  mathematical  model.  Subsequent 
analyses  of  aircraft  and  simulation  flight  dynamics  concentrated  on  the 
hover  and  VA  =  170  knots  cruise  conditions. 

Frequency-sweep  (Sec.  3.1.2)  and  step  inputs  were  executed  in  each 
axis  for  all  four  flight  conditions.  In  hover,  three  frequency-sweeps 
were  conducted  to  ensure  that  sufficient  dynamic  data  were  obtained  for 
good  identification  of  this  most  important  and  (difficult  to  analyze) 
condition.  In  the  remaining  flight  conditions,  where  the  vehicle  is 
much  more  stable,  only  two  frequency-sweeps  were  required. 

The  high  degree  of  open-loop  pitch  and  roll  instability  in  the 
hover  flight  condition  dictated  that  longitudinal  and  lateral  stick 
frequency-sweeps  in  this  flight  condition  be  conducted  with  all  SCAS 
channels  ENGAGED.  Pedal-sweeps  were  conducted  with  the  yaw-SCAS  disen¬ 
gaged,  because  yaw-SCAS  failures  occurred  in  mid-run  resulting  from  the 
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relatively  large  angular  rates  which  were  generated.  The  decoupled  and 
stable  nature  of  the  yaw  (and  heave)  dynamics  allowed  SCAS-off  sweeps  to 
be  conducted  in  these  axes  without  difficulty.  Step  inputs  were  exe¬ 
cuted  in  the  SCAS-on  and  SCAS-off  configurations.  SCAS-off  step 
responses  were  completed  to  verify  that  the  extracted  transfer-functions 
reflect  the  open-loop  response  characteristics  and  not  those  of  the 
inverse  feedback  dynamics  (Chap.  3).  These  inputs  were  completed  with 
the  off-axis  channels  engaged  (e.g.,  open-loop  pitch-axis  inputs  were 
conducted  with  roll  and  yaw  channels  engaged),  to  reduce  the  level  of 
coupled,  unstable  dynamic  response.  SCAS-on  step  responses  were  also 
completed  because  the  SCAS  helps  to  steady  the  initial  conditions,  which 
is  helpful  in  exposing  small  differences  between  the  model  and  aircraft 
dynamics. 

In  the  transition  and  cruise  flight  conditions,  the  longitudinal 
frequency-sweeps  were  conducted  with  pitch  SCAS  engaged .  Frequency- 
sweeps  in  the  lateral  and  directional  axes  were  completed  with  the  roll 
and  yaw  SCAS  channels  disengaged  because  initial  SCAS-on  results  exhib¬ 
ited  unacceptable  levels  of  aileron/rudder  correlation  (Chap.  3).  Step 
inputs  in  all  axes  were  completed  for  both  SCAS-on  and  SCAS-off  config¬ 
urations  as  in  hover. 

The  next  chapter  discusses  the  dynamics-identification  method  in 


detail . 


Chapter  3 

FREQUENCY-DOMAIN  IDENTIFICATION  METHOD 

The  frequency -domain  method  of  Fig.  1.3  was  briefly  outlined  in 
Chapter  1.  In  this  Chapter,  the  individual  analysis  steps  of  the 
approach  are  discussed  in  detail. 

3.1  Flight  Test  Technique 

3.1.1  General  Requirements 

The  principal  ingredient  for  success  in  any  identification  scheme 
is  the  selection  of  an  input  signal  that  excites  the  vehicle  in  all  of 
its  dominant  modes  of  motion  in  the  frequency  range  of  interest.  In 
piloted  handling-qualities  studies,  the  important  frequency  range  is 
approximately  0. 2-6.0  rad /sec  (Refs.  17,  19).  To  achieve  good 
frequency-domain  identification  in  this  range,  the  magnitude  of  the 
input  autospectrum  should  be  nearly  constant.  From  a  piloting  stand¬ 
point,  the  selected  inputs  should  be  easily  repeatable  and  not  involve 
drastic  maneuvers  or  significant  changes  in  the  equilibrium  flight 
condition.  One  input  which  fulfills  the  preceding  requirements  and 
which  has  been  used  successfully  in  a  number  of  rotorcraft  flight  tests 
(Refs.  6,  18,  20,  21,  26)  and  non-rotorcraft  flight  tests  (Refs.  12,  13, 
27)  is  the  ''frequency-sweep." 

3.1.2  Frequency  Sweep 

A  typical  lateral-stick  ( 5lat)  frecluency"sweeP  completed  during  the 
XV-15  hover  tests  and  developed  by  the  author  is  shown  in  Fig.  3.1. 
Pilot-generated  inputs  were  used  instead  of  computer-generated  inputs. 
The  sweep  is  initiated  with  two  low-frequency  sinusoid-shaped  cycles, 
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Fig.  3.1.  Typical  Lateral-Stick  Frequency  Sweep 

with  periods  corresponding  to  the  lower  bound  of  the  desired  identifica¬ 
tion  frequency  range.  These  cycles  ensure  good  excitation  of  the  low- 
frequency  dynamics.  The  desired  lower  bound  of  “mjn  =  0.2  rad/sec 
gives  a  low-frequency  period  of  about  T__v  =  30  sec.  However,  attempts 
to  execute  this  low-frequency  input  in  flight  resulted  in  undesirably 
large  vehicle  motions.  Therefore,  it  was  necessary  to  reduce  the  period 
of  the  low-frequency  cycle  to  Tjj^  =  20  sec  (Fig.  3.1).  The  magnitude 
of  the  control  input  is  adjusted  to  keep  roll  and  pitch  attitudes,  as 
well  as  resulting  translations,  comfortable  and  reasonably  close  to 
their  trim  values.  The  lower  pitch  and  roll  damping  in  hovering  flight 
demands  smaller  pilot  inputs  than  those  which  are  used  in  transition  and 
forward  flight. 

After  the  initial  two  low-frequency  cycles,  the  control  is  moved  at 
gradually  increasing  frequency  for  another  50  sec.  This  increasing- 
frequency  (as  opposed  to  a  decreasing-frequency)  input  allows  the  tran¬ 
sients  of  the  low-frequency  modes  to  persist  for  a  few  more  cycles 
during  the  remainder  of  the  run,  thus  improving  the  low-frequency 


identification;  the  high-frequency  modes  need  less  excitation  time  to 
persist  the  same  number  of  cycles,  so  they  are  excited  toward  the  end  of 
the  input.  By  the  end  of  the  run,  the  control  is  being  driven  quite 
rapidly  (about  4  Hz  in  Fig.  3.1)  with  generally  smaller  displacements. 
The  control  is  then  returned  to  trim,  ending  the  approximately  90-sec 
test  period.  The  selection  of  the  (rather  long)  overall  90-sec  test 
period  was  found  (empirically)  to 

1 .  Allow  a  good  identification  of  the  low-frequency  modes 

2.  Result  in  an  even  pilot-excitation  of  the  vehicle  dynamics  over 
the  identification  frequency  range 

3.  Yield  a  large  quantity  of  data,  thereby  reducing  the  variance 
of  the  spectral  estimate  (see  Sec.  3.2.3) 

At  least  two  sweeps  are  performed  consecutively  for  each  control  to 
ensure  a  sufficient  amount  of  good  quality  dynamic  data  for  the  spectral 
analysis.  Frequency-sweeps  are  conducted  individually  on  each  of  the 
four  primary  aircraft  controls  (lateral  and  longitudinal  stick,  pedals, 
and  power  lever)  for  one  flight  condition,  before  other  conditions  are 
investigated. 

3.2  Frequency-Response  Identification 

The  most  important  step  in  the  frequency-domain  identification 
method  is  the  accurate  estimation  of  input-to-output  spectral  responses, 
since  all  succeeding  calculations  are  based  on  these  results.  Consider¬ 
able  effort  was  invested  at  the  beginning  of  this  study  in  evaluating 
and  implementing  various  spectral-analysis  methods.  Numerical  studies 
based  on  reconstructing  frequency  responses  from  simulated  time- 
histories  of  known  transfer  functions  were  useful  in  evaluating  the 


various  computational  methods.  However,  methods  which  proved  successful 
with  simulated  data  did  not  always  give  satisfactory  results  when 
applied  to  flight  data.  Good  quality  results  at  low  frequency  initially 
proved  very  difficult  to  obtain.  This  was  due  to  the  poor  spectral 
content  which  generally  occurs  at  low  frequency  because  of  the  lower 
pilot  input  amplitudes  (see  Fig.  3.1).  Considerable  improvements  were 
achieved  with  more  advanced  analysis  methods  developed  as  a  result  of 
identification  efforts  on  the  XV-15,  CH-47,  and  Bell  214-ST  aircraft 
(Refs.  6,  18,  20-24).  These  methods  draw  heavily  on  the  signal  process¬ 
ing  literature  (Refs.  28-32). 

Figure  3.2  shows  the  computational  procedures  which  are  performed 
by  the  FRESPID  program  to  generate  the  spectral  responses  of  selected 
input  and  output  time-histories.  This  program  was  jointly  developed  by 
the  author  and  J.  G.  M.  Leung. t  The  FRESPID  program  operates  in  a  batch 
processing  mode,  under  the  control  of  a  command  file.  As  illustrated  in 
Fig.  3-3,  this  file  contains  command  directives  and  flight-data  record 
identifiers  which  specify  the  desired  input/output  records  and  data 
processing  options.  The  next  two  sections  describe  in  more  detail  the 
computations  which  are  performed  by  the  FRESPID  program. 

3.2.1  Data  Preparation 

As  indicated  in  the  flowchart  of  Fig.  3.2,  a  large  part  of  the 
analysis  involves  data  preparation  before  applying  the  Fourier  transform 
algorithm.  These  preparation  steps  are  crucial  to  obtaining  an  accurate 
frequency-response  estimate. 

t  Electronics  Engineer;  Flight  Experiments  Support  Branch,  NASA  Ames 
Research  Center,  Moffett  Field,  Calif. 


Command  file  directives 


Remarks 


V-O.  P/DAIL  883,884  60S  WINDOW  [OC128] 
WINDOW 

MYDISK: [TISCHLER . TILTBFILE1XV15W60S 
LINK 

D645R883 . 150  V012R883.150 
D645R884 . 150  V012R884 . 150 
FILTER 

MYDISK: [TISCHLER. BFILE] AALIAS . 62H 

0.1,  10.0 

0C128 

STOP 


command  file  identifier 
key  word:  window  specification 
window  specification  file 
key  word:  file  concatenation 
input  file-1  output  file-1 
input  file-2  output  file-2 
key  word:  digital  filtering 
impulse  response  file 
Chirp  z-transform  freq.  range 
output  file  name 
program  terminator 


Fig.  3-3-  Sample  Command  File  for  FRESPID. 


The  dc  (offset)  components  and  linear  drifts  (trends)  of  the 
selected  input  and  output  time-histories  are  first  removed  to  prevent 
oscillation  in  the  spectral  calculations.  Since  the  frequency-sweep 
begins  and  ends  in  trim  and  is  a  roughly  symmetric  excitation 
(Fig.  3.1),  the  dc  and  drift  are  simply  determined  from  the  mean  offset 
and  slope  of  the  time-history  trace  (Ref.  28).  When  (nr)  multiple  runs 
of  good  spectral  quality  are  available,  they  are  concatenated  to  form 
composite  input  and  output  time-histories.  The  concatenated  time- 
histories  are  then  digitally  filtered  to  eliminate  the  jitter  (digital 
noise)  which  arises  from  asynchronous  sampling  of  the  instruments 
("sample-skewing"),  and  can  cause  distortion  in  the  spectral  estimate. 
The  filter  cutoff  frequency  is  62.5  Hz,  which  is  50%  of  the  Nyquist 
frequency.  Finite- impulse-response  (FIR)  filters  are  used  [instead  of 
recursive  infinite-impulse-response  (HR)  filters]  because  they  have 
magnitude  characteristics  that  can  be  easily  tailored,  and  they  exhibit 
a  linear  phase  shift  with  frequency  (independent  of  the  magnitude 
response)  (Ref.  28).  The  additional  computational  time  for  the  FIR 


filters  compared  to  that  of  IIR  filters  is  not  a  concern  for  off-line 
analyses. 

The  filtered  concatenated  records  are  partitioned  into  (K)  tapered 
and  overlapping  subrecords  of  (L)  discrete  data  points.  These  sub¬ 
records,  or  time-history  sections,  are  often  referred  to  as  "windows," 
because  the  data  contained  therein  are  the  only  part  of  the  time-history 
that  the  fast  Fourier  transform  (FFT,  Sec.  3.2.2)  can  "see"  at  one 
time.  Applying  a  cosine-squared  tapering  (weighting)  function  (also 
referred  to  as  "Hanning  window")  prevents  side  lobes  and  leakage  which 
would  otherwise  occur  with  a  rectangular  (nontapered)  window.  Over¬ 
lapping  the  windows  makes  more  efficient  use  of  the  data  at  the  edges  of 
the  tapered  window,  thereby  reducing  the  spectral  bias  and  variance 
compared  with  the  non-overlapped  tapered  estimate  (Ref.  28). 

The  width  of  the  window  (LAt,  sec;  At  is  the  time  increment)  is  an 
important  parameter  in  determining  the  quality  of  the  spectral  identifi¬ 
cation.  The  windows  are  optimized  for  each  run  (and  flight  condition) 
to  achieve  maximum  spectral  accuracy  over  the  desired  frequency  range. 

As  discussed  in  Sec.  3.2.3,  the  random  error  is  inversely  proportional 
to  the  square  root  of  the  number  of  (independent)  time-history  sec¬ 
tions — which  is  the  effect  of  averaging.  So  for  a  given  total  run  time, 
smaller  windows  yield  more  averaging  and  lower  spectral  variance.  How¬ 
ever,  this  reduction  in  window  length  also  reduces  the  low-frequency 
content  in  each  window  and,  accordingly,  in  the  overall  identifica¬ 
tion.  Thus,  the  selection  of  window  length  involves  a  compromise 
between  a  high  number  of  averages  and  adequate  low-frequency  signal 
content.  By  concatenating  good  repeat  runs  into  a  single  time-history, 
longer  windows  can  be  selected  without  compromising  the  number  of 


averages.  This  yields  wide-band  identification  and  low  spectral 
variance.  When  only  one  suitable  run  is  available,  shorter  windows  are 
necessary  to  maintain  a  sufficient  number  of  averages,  and  the  low- 
frequency  identification  is  compromised.  A  window  length  which  contains 
roughly  two  periods  of  the  lowest-frequency  wave  of  interest 
(LAt  =  2Tmax  =  40  sec)  is  generally  a  good  compromise. 

3.2.2  Spectral  Functions^- 

The  finite  Fourier  transform  X(f,T)  of  a  (finite  length,  T)  con¬ 
tinuous  time-history  subrecord  x(t)  (e.g.,  one  window  of  time-history 
data)  is: 


X(f,T)  =  f  x(t)e" 

no 


j2nft 


dt 


(3.1) 


When  the  time-history  data  are  from  a  discrete  sequence  xn,  the 
Fourier  transform  integral  relation  is  approximated  by  the  discrete 
Fourier  transform  (DFT): 


N-1 

X(fk)  =  X(kAf)  =  At  £  xn  exp[-J2n(kn)/N]  ;  k  =  0,  1 ,  2, . . . ,  N-1 

n=0  n 

(3.2) 


where 

X(f^)  =  Fourier  coefficients 
xn  =  x(nAt) 

At  =  time  increment 

N  =  number  of  discrete  frequency  points 

t  Background  information  on  spectral  analyses  and  the  chirp  z-transform 
is  given  in  Appendix  A. 


N  =  L  =  number  of  time-history  points  in  the  sequence,  xn  (3.3) 

so, 

T  =  LAt  =  time  duration  of  sequence  (i.e.,  "width  of  window") 

(3.4) 

Af  =  ^  =  frequency  resolution,  Hz 

The  evaluation  of  Eq.  (3.2)  requires  approximately  N  multiply-add 
operations  for  each  window  of  data. 

The  fast  Fourier  transform  (FFT)  refers  to  efficient  algorithms 
(there  are  many,  Ref.  28)  for  determining  the  DFT  of  Eq.  (3-2).  The  FFT 
procedures  require  approximately  4N  log2  N  computations,  which  gives  a 
speed  ratio  of  (Ref.  28): 

N2  • 

speed  ratio  =  4^  iQg  ^  =  160  for  the  present  analysis 

indicating  a  considerable  computational  savings.  One  particularly 
flexible  FFT  algorithm  is  the  chirp  z-transform  (also  referred  to  as 
the  "zoom  transform"  or  CZT).  This  algorithm  (Refs.  30,  32)  allows  the 
rather  severe  restriction  of  Eqs.  (3.3)  and  (3.4)  to  be  relaxed,  which 
yields  much  finer  frequency  resolution  (Af)  for  a  given  window  size 
(LAt).  Another  key  advantage  of  the  chirp  z-transform  is  that  the 
number  of  time-history  points  t  can  be  arbitrary,  whereas  L  must  be 
a  composite  integer  (an  integer  power  of  2)  in  standard  FFT  proce¬ 
dures.  These  capabilities  of  the  chirp  z-transform  allow  the  extraction 


of  high-resolution  spectra  in  a  narrow  frequency  band,  and  an  increase 
in  the  identified  dynamic  range,  especially  at  the  low-frequency  end 
(Ref.  32).  The  chirp  z-transform  has  been  found  in  this  and  in  a  number 
of  other  studies  by  the  author  to  be  well  suited  for  frequency-response 
identification  from  flight  data  (Refs.  6,  20,  21,  23). 

The  one-sided  Input  autospectral-density  estimate  for  the  sub¬ 
record  xn  is  determined  from  the  chirp  z-transform  (Fourier) 
coefficients : 


5xx(V  =fiJ  lX(fk>|2  •  k  =  0,  1,...,  N/2  (3.5) 

where  U  is  the  scale  factor  for  window  tapering  [U  =  (8/3)  for 
Hanning  window,  Ref.  28], 

The  one-sided  output  autospectral-density  estimate  is  similarly 
obtained  from  the  output  subrecord  by  replacing  y  (output  record)  for 
x  in  Eqs .  (3.1M3.5): 

Gyy(fk)  =  fjj  |Y(fk)|2  ,  k  =  0,  1,...,  N/2  (3.6) 

The  one-sided  cross-spectral  density  estimate  is  determined  byf : 

GXy(fk)  =  |j  [X*(fk)Y(fk)J  ,  k  =  0,  1 . N/2  (3.7) 


t  Asterisk  (*)  denotes  complex-conjugate. 


Finally,  the  total  spectral-function  estimates  for  the  entire 
(concatenated)  time-history  is  obtained  from  a  linear  average  of  the 
spectra  for  the  (K)  overlapped  subrecords: 

Gxx(fk)  =  averaged  input  autospectral-density  estimate 

G.,„(f|.)  =  averaged  output  autospectral-density  estimate 

yy  K 

GXy(fk)  =  averaged  cross-spectral  density  estimate 

These  quantities  are  physically  interpreted  as  giving  the 
distribution  of  the  mean-squared  response  of  the  respective  signals 
(i.e.,  xx,  yy,  xy)  as  a  function  of  frequency.  For  illustrative 
purposes,  the  spectral-density  magnitudes  are  presented  in  power- 
dB  (=  10  log10|Gxx(fj{)  | ,  etc.)  which  gives  the  distribution  of  the  root 
mean-squared  (rms)  response.  These  results  are  the  bases  for  the 
frequency-response  calculations  in  the  next  section. 

3.2.3  Single-Input/Single-Output  Frequency-Response  Calculations 
in  the  Open-Loop 

Once  the  input,  output,  and  cross-spectral  density  estimates 
[(Gxx(fk)*  Gyy(fk),  Gxy^*V»  respectively]  have  been  determined  for  a 
selected  time-history  pair  (x,  y  in  Fig.  3.4)  the  single- input/single¬ 
output  frequency  response  is  estimated  by 


Fig.  3.4.  Single-Input/Single-Output  Open-Loop  Response 
Identification . 


For  the  remainder  of  this  dissertation,  the  estimate  symbol  (*)  and 
the  discrete  frequency  dependency  (f^)  are  omitted,  but  they  are  implied 
throughout.  The  transfer-function  results  are  presented  in  standard 
Bode  form;  that  is,  a  plot  of  log-magnitude  (dB  =  20  log10|H|)  and  phase 
(deg)  of  H  versus  log-frequency  (ui  =  2n f^,  rad/sec). 

In  the  context  of  nonlinear  systems  analysis,  the  result  presented 
in  Eq.  (3.8)  is  a  describing  function  since  it  relates  that  part  of  the 
output  which  can  be  linearly  related  to  the  input  (Ref.  33).  In  other 
words,  this  calculation  results  in  an  equivalent  linear  model,  which 
minimizes  the  mean  square  difference  between  the  actual  output  signal 
and  its  approximation  by  the  fundamental  harmonic  (first  sinusoidal 
component  of  the  Fourier  series).  The  remaining  output  harmonics  are 
called  remnants.  Most  of  these  are  largely  filtered  out  by  the  low-pass 
nature  of  the  aircraft  dynamics. 

The  piloting  technique  and  time-history  concatenation  methods 
described  above  result  in  input  amplitudes  that  are  not  constant  over 
the  frequency  range.  Therefore,  the  resulting  describing  function  is 
representative  of  the  average  input  amplitude  over  the  ensemble  time- 
history.  Since  the  pilot  inputs  are  typical  of  those  experienced  during 
normal  flight  operations,  this  approximation  produces  a  satisfactory 
description  of  the  vehicle  dynamics  in  its  normal  operating  environment. 

A  good  indication  of  the  quality  of  the  first  harmonic  as  a  model 

of  the  input-to-output  dynamics  is  obtainable  from  the  coherence  func- 
2 

tion  (y  )  defined  by 
xy 


which  may  be  interpret'd  as  that  fraction  of  the  output  spectrum  that 
can  be  accounted  for  by  linear  relation  with  the  input  spectrum 
(Ref.  29).  When  the  process  under  investigation  is  perfectly  linear  and 
the  spectral  estimates  are  noise  free,  the  coherence  function  will  be 
unity  within  the  excitation  frequency-range.  (The  meaningless  result 

p 

of  y  =  1  also  occurs  when  only  one  window  section  is  used — that  is, 
no  averaging). 

A  coherence  function  of  less  than  unity  results  from  three  basic 
causes  (Refs.  28,  34).  First,  nonlinearities  in  the  system  may  produce 
remnants  which  are  not  accounted  for  by  the  first  harmonic  approxima¬ 
tion.  When  the  input  and  output  excursions  are  small,  the  system  non- 
linearities  are  less  important.  This  is  generally  true  for  the  cruise 
flight  conditions.  However,  in  hover,  where  low-frequency  inputs  induce 
large  motions  relative  to  the  trim  condition,  nonlinear  effects  become 
more  important;  a  significant  reduction  in  the  coherence  function  will 
result  when  first  harmonic  approximations  are  inadequate  for  represent¬ 
ing  these  nonlinearities.  In  the  higher-frequency  range,  where  the 
excursions  are  smaller,  the  coherence  function  again  approaches  unity. 

A  second  and  very  common  cause  of  reduced  coherence  is  input  and 
output  noise.  The  effect  of  noise  sources  differs  chiefly  depending  on 
whether  they  are  correlated  with  the  input  and  output  signals,  and 
whether  they  act  on  the  "process"  or  on  the  "measurement."  Process 
noise  (e.g.,  turbulence)  contributes  to  the  aircraft  response,  and  so 
acts  the  same  as  a  secondary  (but  unmeasurable!)  input;  the  problem  of 


I 

secondary  inputs  is  discussed  below.  Measurement  noise  does  not  con¬ 
tribute  to  the  aircraft  response,  but  it  does  contaminate  the  recorded 
|  input  and  output  signals.  Process  and  measurement  noise  are  important 

•  because  they  cause  errors  in  the  estimate  of  the  frequency  response. 

These  errors  are  categorized  as  "bias"  (systematic)  or  "random"  (nonsys- 
tematic).  Bias  error  exists  when  the  expected  value  of  the  ensemble- 
I  averaged  frequency-response  estimate  E[H(f)]  is  different  from  the 

actual  frequency  response  H(f)  or 

j  b[H(f ) ]  =  E[H(f ) ]  -  H(f) 

» 

i 

;  Random  error  exists  when  there  is  a  dispersion  of  the  frequency-response 

5  estimates  about  the  expected  value  for  a  collection  of  repeated  runs: 

) 


o [ H ( f ) )  =  ♦  ^E[(H(f)  -  E[H])2] 

Measurement  noise  is  generally  considered  to  be  uncorrelated  with 
the  measured  signal.  Uncorrelated  measurement  noise  at  the  output  [ v( t ) 
in  Fig.  3.5]  does  not  bias  the  spectral  estimate  of  Eq.  (3.8).  However, 
uncorrelated  measurement  noise  at  the  input  [u(t)  in  Fig.  3.5]  does  bias 
the  estimate  as  a  function  of  the  noise-to-signal  ratio  (Ref.  28);  so. 


Fig.  3.5.  Uncorrelated  Measurement  Noise  at  the  Input  and  Output. 


input  sensors  must  be  of  the  highest  possible  quality  (ideally  noise- 
free).  Correlated  measurement  noise  at  the  input  or  output  or  both  will 
also  bias  the  estimate.  Uncorrelated  measurement  noise  also  causes 
random  error,  which  is  reflected  by  the  drop  in  the  coherence  func¬ 
tion.  This  coherence  drop  is  often  noticeable  in  the  hover  flight 
condition,  where  the  pilot's  low-frequency  inputs  are  generally  small. 
This  also  often  occurs  at  higher  frequency  (above  about  1  Hz),  where  the 
aircraft  response  is  generally  small.  For  a  given  level  of  coherence,  a 
significant  improvement  in  the  spectral  accuracy  is  obtained  when  multi¬ 
ple  runs  are  concatenated,  since  the  random  (normalized)  error 
(ep  =  o[H]/H)  is  inversely  proportional  to  the  number  of  independent 
time-history  averages  (Ref.  34): 


B 


c 


r 


(t) 


(3.10) 


and 


nd 


LAt 


=  number  of  independent  time-history  averages 


where 

Tp  =  Total  concatenated  run  length,  sec 
L  =  Number  of  time-history  points  in  the  window 
At  =  Time  increment 


t  The  factor  of  . 55  accounts  for  the  45 %  reduction  in  spectral 
variance  due  to  the  50 %  window  overlapping  used  in  this  study  (Ref.  31 


). 
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Based  on  the  average  window  length  2Tmax  =  40  sec,  and  two  concatenated 

90-sec  runs,  a  random  error  of  er  <  0.1  (less  than  1 0% )  is  achieved  for 

2 

a  coherence  function  of  y  >0.85. 

xy 

The  third  source  of  reduced  coherence  results  from  secondary 
inputs;  these  are  defined  as  inputs  (other  than  the  control  which  is 
undergoing  a  frequency-sweep,  x  in  Fig.  3.6)  which  excite  the  system. 


SECONDARY 

INPUTS 

Fig.  3.6  Effect  of  Secondary  Inputs  on  Single-Input/Single-Output 
Identification. 


Secondary  inputs  occur  from  gusts  and  turbulence  (process  noise),  and 
from  off-axis  control  activity  of  the  pilot  or  SCAS.  The  gust  and 
turbulence  inputs  cannot  be  measured,  they  contribute  to  the  measured 
response  of  the  system,  and  they  generally  evoke  correlated  regulatory 
inputs  from  the  pilot  and  the  SCAS;  therefore,  the  extracted  frequency- 
responses  are  biased  (see  Sec.  3.2.4).  Furthermore,  the  random  error 
increases.  Therefore,  frequency-sweeps  must  be  conducted  in  conditions 
of  minimum  wind  and  turbulence.  The  author's  experience  indicates  that 
hover  tests  must  be  conducted  with  steady  winds  that  are  no  greater  than 
5  knots  and  that  cruise  tests  must  be  conducted  with  turbulence  that  is 
no  greater  than  ±1  knot.  The  residual  effects  of  small  levels  of  gusts 
and  turbulence  are  minimized  by  the  concatenation  of  multiple  frequency- 


sweep  records.  Pilot  or  SCAS  inputs  or  both  for  regulation  against 


excursions  which  occur  in  the  coupled  degrees-of-freedom  also  contribute 
to  the  system  response;  however,  these  inputs  are  measurable.  If  the 
secondary  and  primary  inputs  are  not  fully  correlated,  their  effects  can 
be  separated  by  using  the  multi-input/multi-output  identification  method 
described  in  Sec.  3.2.5. 

3. 2. 4  Identification  of  the  Open-Loop  Single-Input/Single-Output 
Frequency  Response  From  Closed-Loop  Tests 
When  the  bare-airframe  characteristics  of  the  vehicle  under  study 
are  highly  unstable,  as  are  those  of  the  XV- 15  in  hover,  execution  of  a 
90-sec  frequency-sweep  on  the  open-loop  configuration  is  not  practi¬ 
cal.  The  single- input/single-output  bare-airframe  frequency-response 
can  be  extracted  from  the  closed-loop  flight  data  using  Eq.  (3.8)  sub¬ 
ject  to  an  important  qualification.  Referring  to  Fig.  3.7,  the  closed- 
loop  system  responds  to  commanded  inputs  ( AT )  and  Process  noise  (n). 
Identification  of  the  open-loop  frequency  response  (p/«a)  requires  the 
measurement  of  the  output  (p)  and  the  surface  deflection  (6_) 

<4 


(SCAS) 


Fig.  3.7.  Single  Degree-of-Freedom  Closed-Loop  Roll-Response  Model. 


signals.  This  total  surface  deflection  is  made  up  of  components  from 
the  pilot  ( 5LAT ^  and  the  ^CAS  feedback  (6^  ).  Thus  the  identified 
frequency  response  of  p/6_  is  also  expressible  as 

cl 


When, 


pGF  >>  6LAT> 
P 


(3.11) 


(3.12) 


which  is  the  troublesome  result  of  identifying  the  inverse  feedback 
frequency-response  in  closed-loop  pilot  tracking  tasks  without 

a  significant  known  external  input  (Ref.  35).  The  limiting  condition  of 
Eq.  (3.12)  is  reachable  only  in  the  case  of  system  excitation  by  process 
noise  (n)  (e.g.,  turbulence)  with  no  pilot  inputs  ($LAT  =  0).  However, 
biases  in  the  open-loop  frequency-response  estimate  obtained  from 
Eq.  (3-8)  will  occur  even  with  non-zero  pilot  input  («Lat  *  whenever 
process  noise  is  present;  this  occurs  since  the  feedback  loop  causes  the 
surface  deflection  to  be  partially  correlated  with  the  process  noise 
(Appendix  B).  The  amount  of  bias  error  depends  on  the  noise-to-signal 
rms  ratio  o(uj)  at  each  frequency;  the  numerical  study  presented  in 
Appendix  B  shows  that  the  (normalized)  bias  errors  become  significant 

A 

(eb  =  b[H]/H  >  0.10;  greater  than  10% )  when  the  system  rms  excitation 
that  is  due  to  process  noise  is  greater  than  33/S  of  that  that  is  due  to 
pilot's  inputs  [o(u)  >  0.33].  Total  error  (composed  of  normalized 
random  and  bias  error)  becomes  significant  (greater  than  10J)  at  lower 

40 


noise-to-signal  rms  ratios  [0(01)  >  0.12)  depending  on  the  level  of 
feedback  gain. 

Therefore,  the  key  requirement  for  identifying  the  open-loop  fre¬ 
quency  response  from  closed-loop  flight  data  is  that  the  total  surface 
deflection  (6a)  must  contain  a  significant  component  from  the  pilot 
(6laj)  which  is  uncorrelated  with  the  output  in  the  frequency  range  of 
interest.  The  "programmed”  frequency-sweep  input  has  good  input  power 
across  the  entire  frequency  range  of  concern,  and  persists  throughout 
the  duration  of  the  time-history;  as  a  result,  it  is  especially  good  for 
closed-loop  testing.  Multi-step  inputs  are  not  particularly  well  suited 
for  open-loop  identification  in  closed-loop  tests  because  the  excitation 
occurs  principally  at  the  beginning  of  the  run.  Most  of  the  long-term 
response  is  excited  by  fully-correlated  inputs  from  the  SCAS.  The  use 
of  the  frequency-sweep  input  and  conducting  the  flight  tests  in  low- 
turbulence  conditions  ensures  that  an  accurate  open-loop  identification 
is  achieved.  SCAS-off  step  responses  of  the  vehicle  are  compared  with 
the  extracted  transfer-function  model  responses  to  verify  that  the  true 
(unbiased)  bare  airframe  frequency  response  is  identified. 

Identification  of  the  unstable  bare-airframe  dynamics  is  thus 
greatly  simplified  since  the  pilot  is  flying  the  closed-loop  (stable) 
aircraft.  When  turbulence  is  present,  the  remaining  (secondary) 
degrees-of-freedom  will  be  excited.  For  example,  pitch  and  yaw  motions 
will  generally  be  present  during  lateral  stick  frequency-sweeps.  How¬ 
ever,  if  the  roll  response  is  completely  decoupled  from  the  pitch  and 
yaw  responses,  the  identified  roll  frequency-response  and  coherence 
functions  will  not  be  affected.  For  such  aircraft,  the  pilot  or 
stability-augmentation  systems  should  be  used  to  regulate  against 


introducing  disturbances  into  the  secondary  degrees-of- freedom.  A 
problem  arises,  though,  when  coupled  aircraft  dynamics  are  considered. 
3.2.5  Multi-Input/Multi-Output  Frequency  Response  Identification 
in  the  Closed-Loop 

Consider  the  case  of  coupled  roll-yaw  dynamics  with  SCAS  feedbacks 

in  both  channels,  as  illustrated  in  Fig.  3.8.  Following  the  preceding 

discussion,  we  assume  that  the  total  surface  deflections  (6_  and  6  ) 

d  r 

contain  substantial  components  from  the  pilot  inputs  (6^^  and  fipgp) 

which  are  uncorrelated  with  the  output  responses. 

Now  suppose  that  a  lateral-stick  sweep  is  being  executed  to  obtain 

the  open-loop  transfer  function  for  yaw  response  to  aileron  (r/6a). 

Unlike  the  previous  decoupled  case,  aileron  inputs  induce  coupled 

responses  in  the  yaw  degree-of-freedom.  If  SCAS  and  pilot  yaw  feedbacks 

(Gf  )  act  to  restrain  these  coupled  yaw  excursions,  the  identified 
rr 

transfer  function  r/$a  will  be  a  greatly  distorted  indication  of  the 
actual  bare-airframe  transfer  function. 

In  a  general  sense,  yaw-rate  feedback  decouples  the  yaw  and  roll 
dynamics^;  this  may  be  desirable  for  assessing  the  open-loop  roll 
response  (p/60)  using  the  simplified  one  degree-of-freedom  model 

a 

(Fig.  3.4).  However,  when  the  interaxis-coupling  characteristics  are 
sought  (p/6„,  r/6_),  artificial  decoupling  is  not  satisfactory.  In 
hovering  rotorcraft,  rotor  flapping  dynamics  induce  interaxis  coupling 
effects  which  may  degrade  piloted  handling-qualities  (Ref.  36);  there¬ 
fore,  the  identification  of  interaxis  coupling  responses  is  a  specific 
concern . 

t  This  assumes  that  L,  is  small,  as  it  is  in  the  present  XV-15  case. 
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Fig.  3.8.  Coupled  Roll/Yaw  Response  Model. 


The  two-input/single-output  identification  problem  illustrated  in 
Fig.  3.8  is  a  special  case  of  the  more  general  multi- input/single-output 
analyses  discussed  in  Refs.  10,  28,  and  29.  The  present  analysis  is 
limited  to  the  two-input  problem  since  the  side-by-side  rotor  configura¬ 
tion  of  the  XV- 15  aircraft  induces  coupling  between  the  roll  and  yaw 
degrees-of-freedom  only  (Chap.  2).  (Identification  of  open-loop 
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dynamics  for  single-rotor  helicopters  may  require  consideration  of  more 
than  two-control  inputs . ) 


The  components  of  the  yaw-rate  spectrum  resulting  from  the  aileron 
and  rudder  inputs  are  determined  from 
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(3.13) 


(3.14) 


The  desired  frequency  responses  are  yaw-rate  response  to  aileron  (r/6a) 
and  yaw-rate  response  to  rudder  (r/6r).  These  are  given  by  the  simulta¬ 
neous  solution  of  Eqs.(3.13)  and  (3.14): 


r_ 

6 


G6  r[1  ■  (GS  6  G6  r/G6  6  G6  r)] 
a _ a  r  °r  r  r  a 

,  <’  -  4  4  > 
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(3.15) 


r"  -  <GS  «  G5  r/G«  5  GS  r>> 
r  r  raa  aar 
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(3.16) 


C6  6  (1  -  Y6  6  > 
r  r  a  r 


It  is  clear  from  the  above  equations  that  the  aileron  and  rudder 

inputs  must  not  be  fully  correlated  since  this  would  result  in  an  input 
2 

cross-coherence  (y,  .  )  of  unity,  thereby  rendering  the  desired  fre- 

a  r 

quency  responses  undefined  (the  "multipath  propagation  problem", 

Ref.  34).  Then,  only  the  (assumed)  decoupled  roll-rate  response  to 
aileron  can  be  identified.  Note  also  that  if  the  inputs  are  fully 
uncorrelated.  Eqs.  (3.15)  and  (3.16)  reduce  to  the  single-input/single¬ 
output  relations  from  Eq.  (3.8). 


Equations  (3.15)  and  (3  16)  can  be  written  more  compactly  by 


defining 


Gf  ,  =  G  fl  -  (C.  .  G,  /G .  .  G„  )] 

6  r-6  6  r  6  6  6r  66  6  r 
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Each  quantity  defined  in  Eqs.  (3- 17)~(3. 20)  has  important  physical 

significance.  For  example,  G  is  the  cross-spectrum  from  aileron 

ar  r 

input  to  yaw-rate  response  (6a  to  r)  with  the  linear  effect  of  yaw-rate 
due  to  rudder  excitation  removed.  This  quantity  is  termed  the  "condi¬ 
tioned  (or  residual)  cross-spectrum"  because  the  "ordinary"  cross¬ 
spectrum  (Gg  )  has  been  conditioned  to  remove  the  biasing  effects  of 
°ar 

partially-correlated  rudder  excitation.  Similarly,  the  quantity 

G  is  termed  "conditioned  input  autospectrum,"  because  the  ordi- 

a  a  r 

nary  input  autospectrum  (G.  ,  )  has  been  conditioned  to  remove  the 

6  6 
a  a 

partially-correlated  contribution  of  6p  inputs.  The  remaining  quanti¬ 
ties  of  Eqs.  (3.19)  and  (3.20)  have  analogous  interpretations. 

The  desired  frequency  responses  of  Eqs.  (3.15)  and  (3.16)  are  then 
expressed  as: 


6  r  •  6 

- —  s  r/6  -6 
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(3.21) 
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which  are  in  the  form  of  the  single-input/single-output  relation  of 
Eq.  (3.8)  with  cross-spectrum  replaced  by  conditioned  cross-spectrum, 
input-autospectrum  replaced  by  conditioned  input-autospectrum,  and 
frequency  response  replaced  by  conditioned  frequency-response. 

A  generalization  of  the  ordinary  coherence  function  concept  pre¬ 
sented  above  allows  the  quality  of  the  conditioned  frequency-response 
estimates  to  be  assessed.  The  partial  coherence  function  is  defined  as 
that  fraction  of  the  output  that  can  be  linearly  related  to  a  single 
input  with  the  linear  effects  of  all  remaining  inputs  removed.  It  is 
interpreted  in  exactly  the  same  way  as  the  ordinary  coherence  functions 
for  single-input/single-output  systems  described  earlier  Eq.  (3.9),  and 
has  the  same  form.  For  the  coupled  system  of  Fig.  3.6,  the  partial 

coherence  of  yaw  response  to  aileron  with  the  linear  effects  of  the 

2 

rudder  inputs  removed  (y^  ^  )  is  calculated  by 

a  r 
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where  G.  .  and  G.  _  ,  are  defined  in  Eqs.  (3.17)  and  (3.18), 
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respectively;  and  where 
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is  termed  the  "conditioned  output  autospectrum"  because  it  is  the  ordi¬ 
nary  output  autospectrum  (Grr)  with  the  linear  effect  of  partially 
correlated  rudder  excitation  removed.  The  corresponding  partial  coher¬ 
ence  of  yaw-rate  response  to  rudder  with  the  linear  effects  of  aileron 

inputs  removed  (y.  .  )  can  be  obtained  by  interchanging  and  fi,  in 

or • o_  r  a 

Eq.  (3.23). 

Another  useful  metric  is  the  multiple  coherence  function.  This  is 

the  fraction  of  the  output  that  may  be  linearly  correlated  to  all  of  the 

inputs  being  considered. 

( Yg  6  )  is  given  by* 

a  r  * 

2 

Vr 

The  multiple  coherence  will  be  less  than  unity  when  input  and  output 
noise  or  nonlinearities  are  present,  or  when  additional  secondary  con¬ 
trol  inputs  contribute  to  the  aircraft  response  but  are  not  included  in 
the  multi-input  analysis.  As  before,  the  identified  transfer  functions 
will  incorrectly  appear  decoupled  when  these  additional  secondary  inputs 
are  not  taken  into  account. 

The  two-input/single-output  spectral  analysis  calculations  of 
EqS.  (3- 15)— (3-25)  may  be  performed  using  a  computer  program  (TISOSA) 
developed  by  the  author.  The  inputs  to  the  program  are  the  single¬ 
input/single-output  spectral-response  results  of  FRESPID,  and  a  command 
file  which  contains  the  execution  directives.  The  program  outputs  are 
the  conditioned  spectral-response  plots  and  a  data  file. 


For  the  two-input  case,  the  multiple  coherence 


(r^X  r  .  (r/s  )G*  p 
a  r 


(3.25) 


t  Asterisk  (*)  denotes  complex-conjugate. 


3.3  Transfer-Function  Model in 


With  the  flight-test  frequency  responses  in  hand,  the  next  step  in 
the  identification  method  of  Fig.  1.3  is  the  derivation  of  transfer- 
function  models.  This  step  involves  three  important  considerations 
which  are  treated  in  the  next  sections: 

1.  Order  and  form  of  transfer-function  model 

2.  Fitting  range  and  relative  weighting  of  the  magnitude  and  phase 
curves 

3.  Multi-degree-of-freedom  models 

3.3.1  Selection  of  Model  Order  and  Form 

Analytical  transfer-function  forms  are  selected  for  each  input- 
output  pair.  These  are  based  on  configuration  and  flight  condition 
considerations  and  must  have  a  sufficient  number  of  free  parameters  to 
adequately  fit  the  identified  frequency  response  over  the  frequency 
range  of  concern.  In  the  hover  flight  condition,  vehicle  (rigid  body) 
dynamics  are  dominated  by  the  hovering  cubic  and  decoupled  heave  and  yaw 
modes.  In  wing-borne  flight,  the  conventional  longitudinal  and  lateral 
quartic  equations  dominate.  Thus,  transfer-function  models  which  are 
appropriate  for  the  hover  flight  condition  are  not  necessarily  appli¬ 
cable  to  forward-flight  conditions.  Obviously,  if  a  model  of  high 
enough  order  is  selected,  the  parameters  can  be  adjusted  to  accommodate 
each  flight  condition.  However,  such  transfer-function  models  no  longer 
retain  the  physical  significance  of  the  classical  lower-order  parameters 
and  are  often  not  unique  (Ref.  6).  Also,  higher-order  transfer- 
functions  models  tend  to  be  strongly  tuned  to  the  specific  inputs  which 
are  used  in  the  identification  procedures  (e.g.,  a  frequency-sweep)  and 
can  be  very  poor  predictors  of  responses  to  other  test  inputs  (e.g., 


step  and  pulse  inputs)  and  responses  at  nearby  flight  conditions. 
Therefore,  higher-order  models  are  not  desirable.  The  appropriate 
transfer-function  order  is  the  minimum  order  which  can  satisfactorily 
fit  the  frequency  responses,  with  the  upper  limit  taken  as  the  physical 
order  of  the  system. 

3.3.2  Transfer-Function  Fitting  Range  and  Relative  Weighting 

The  maximum  frequency  range  usable  for  transfer-function  fitting 
corresponds  to  the  region  of  acceptable  coherence.  Often,  however,  the 
frequency  range  of  suitable  coherence  extends  beyond  the  range  of  appli¬ 
cability  of  the  selected  transfer-function  models.  For  example,  in  the 
XV-15  cruise  identification  (Chap.  5),  a  simple  short-period  model  is 
assumed,  since  it  adequately  characterizes  the  rigid-body  pitch  dynamics 
in  the  broad  frequency  range  of  0.3-10  rad/sec.  At  the  lowest  input 
frequencies  (0. 1-0.3  rad/sec),  the  vehicle  response  is  dominated  by  the 
lightly  damped  phugoid  dynamics;  these  dynamics  are  not  very  important 
for  piloted  handling-qualities  related  to  attitude  control  tasks  and  so 
can  be  ignored  for  the  present  study.  By  restricting  the  fitting  range 
to  0.3-10  rad/sec,  where  the  short-period  model  is  appropriate,  the 
physical  significance  of  the  parameters  is  maintained.  In  other  cases, 
structural  resonances  may  be  present  at  the  high-frequency  end  of  the 
range  of  good  coherence.  If  a  rigid  body  transfer-function  model  is 
assumed,  the  fitting  range  must  be  restricted  to  exclude  the  effects  of 
such  resonances  (Ref.  6). 

The  interactive  computer  program  NAVFIT,  originally  developed  at 
McDonnell  Aircraft  Company  by  Hodgkinson  and  Buckley  (Ref.  37)  and 
enhanced  by  the  author  to  access  the  output  file  from  the  FRESPID  pro¬ 
gram,  fits  the  frequency-response  data  with  a  selected  transfer-function 


model.  The  numerator  and  denominator  orders  are  arbitrary,  and  a  time 
delay  can  be  included  to  approximate  the  effects  of  unmodeled  high- 
frequency  dynamics.  Least-squares  fitting  is  completed  using  up  to 
50  points  spaced  linearly  across  the  logarithmic  frequency  scale  in  the 
desired  frequency  range.  The  relative  weighting  of  the  magnitude  and 
phase  errors  is  arbitrary,  but  a  commonly  selected  value  for  handling- 
qualities  applications  is  1  dB  magnitude  error  to  7  deg  phase  error^ 
(Ref.  17).  The  outputs  from  the  NAVFIT  program  are  closed-form 
transfer-function  models  and  plotted  frequency-response  comparisons  of 
the  flight  data  and  model. 

3.3.3  Fitting  Multi-Degree  of  Freedom  Frequency  Responses 
A  good  review  of  the  practical  considerations  involved  in  transfer- 
function  fitting  is  given  by  Bischoff  and  Palmer  (Ref.  15).  They  found 
that  satisfactory  identification  of  the  transfer-function  parameters 
generally  requires  the  simultaneous  consideration  of  at  least  two 
degrees-of- freedom.  By  constraining  the  transfer-function  models  to 
contain  mutual  denominator  factors  (as  they  do  physically),  the  location 
of  numerator  factors  becomes  much  better  defined.  In  the  cruise  condi¬ 
tion,  for  example,  the  denominators  of  the  longitudinal  transfer  func¬ 
tions  relating  elevator  inputs  to  normal-acceleration  and  pitch-rate 
outputs  (az/6g  and  q/6e)  are  the  same.  Simultaneously  fitting  these 
degrees-of-freedom  maintains  the  commonality  of  the  denominator  factors 
(the  short-period  mode),  and  avoids  the  "galloping  La"  problem  (this 
refers  to  the  problem  of  a  poorly  defined  pitch-rate  numerator  location 

^Professor  Stephen  P.  Boyd  notes  that  this  is  the  natural  weighting  that 
is  obtained  when  the  fit  is  performed  on  a  real- imaginary  plot  of  the  20  log1Q 
values  of  the  complex  transfer-function. 


in  single-input/single-output  q/«e  identification,  Ref.  17).  Two  com¬ 
puter  programs  (LONFIT  and  LATFIT)  were  developed  at  McDonnell  Douglas 
Aircraft  Company  by  Givan,  LaManna,  and  Hodgkinson  (Refs.  38,  39)  and 
modified  by  the  author  to  allow  multiple-degree-of-freedom  fitting  of 
the  frequency-response  output  files  from  FRESPID.  A  new  option  incorpo¬ 
rated  by  the  author  into  LONFIT  allows  the  identification  of  the  instan¬ 
taneous  center-of-pitch  rotation,  xntT  (Sec.  5.2.2). 

eg 

3.4  Time-Domain  Verification  of  Transfer-Function  Models 

Time-domain  verification  involves  driving  the  extracted  transfer- 
function  models  with  the  flight  data  for  inputs  not  used  in  the  identi¬ 
fication  procedure  and  comparing  predicted  and  actual  time  histories. 
This  analysis  is  necessary  to  check  a  number  of  important  assumptions 
and  steps  in  the  frequency -domain  method: 

Linearity  assumption:  The  adequacy  of  linearized  lower-order 
models  in  describing  the  time-domain  dynamics  of  nonlinear  systems 
is  always  a  key  question--especially  for  helicopters  in  hovering 
flight.  Since  the  reference  condition  is  zero  velocity,  even  small 
translational-velocity  perturbations  can  produce  significant  changes  in 
the  dynamic  response  characteristics.  Lateral  and  longitudinal  fre¬ 
quency  sweeps  can  involve  translational  velocities  of  up  to  40  ft/sec  at 
low  frequencies,  so  a  unique  lower-order  linear  model  may  not  be  achiev¬ 
able  for  the  entire  frequency  range  (Ref.  22).  As  discussed  earlier, 
the  extracted  transfer-function  model  is  actually  a  describing  function 
that  is  strictly  valid  only  for  the  input  amplitudes  which  are  used  in 
the  flight  test  experiment.  However,  since  the  motions  encountered 
during  frequency-sweep  testing  are  typical  of  those  in  near-hover 


maneuvering  flight,  the  extracted  describing  functions  are  generally 
quite  suitable.  Time-domain  verification  checks  this  assumption. 

Transfer-function  fitting  mismatch:  Sometimes,  apparently  small 
mismatches  in  the  frequency-domain  fit  can  produce  surprisingly  large 
discrepancies  in  time-domain  correlations,  especially  when  these  mis¬ 
matches  occur  in  the  frequency  range  near  the  dominant  vehicle  modes. 

If  the  identified  model  poorly  predicts  the  time-domain  verification 
response,  the  transfer-function  model  order  and  structure  is  adjusted  to 
improve  the  correlation. 

Overtuned  models:  Identification  procedures  tend  to  be  tuned  to 
the  type  of  test  inputs  used  (e.g.,  a  frequency-sweep),  and  their  accu¬ 
racy  in  predicting  responses  to  other  classes  of  inputs  (e.g.,  steps 
and  doublets)  is  often  uncertain.  The  possibility  of  overly-tuned 
models  increases  as  the  order  of  the  transfer  function  is  increased 
(Ref.  40).  By  testing  the  transfer-function  model  with  other  input 
forms,  over tuning  is  avoided. 

Identification  of  the  open-loop  frequency  response  from  closed-loop 
flight  data:  Although  a  complete  90-sec  frequency-sweep  cannot  be 
generated  on  the  open-loop  unstable  XV- 15  aircraft  in  hover,  shorter 
duration  open-loop  step  and  pulse  inputs  are  possible.  Comparisons  of 
these  SCAS-off  responses  with  the  extracted  transfer-function  models 
verifies  that  the  true  (unbiased)  open-loop  vehicle  response,  instead  of 
the  inverse  feedback  transfer-function  (at  worst),  has  been  identified. 

3.4.1  Time-Domain  Verification  Method 

Time-domain  verification  may  be  done  by  converting  the  transfer 
functions  to  canonical  form: 


x(t)  =  Fx( t)  ♦  GUp(t) 


(3.26) 


y(t)  =  Cx(t) 


where 


and  aQ,  a.,  .  .  an-1  are  the  denominator  (ascending  coefficients) 
and  bQ,  .  .  . ,  bm  are  the  numerator  (ascending)  coefficients  of 
the  transfer-function  model  being  evaluated. 

The  control  input  to  the  system  Up  is  the  low-pass  filtered  time- 
history  of  the  surface  deflection  signal  (e.g.,  «a  in  Fig.  3.7).  The 
filter  cutoff  frequency  is  selected  to  eliminate  spectral  components  of 
the  test  inputs  (steps  and  pulses)  which  are  above  the  frequency  range 
of  applicability  of  the  derived  transfer-function  models  (e.g.,  the 
upper  bound  of  the  least-squares  fitting  range,  Sec.  3.3.2)  (Ref.  6). 

State  derivative  biases  and  output  reference  values  are  included  as 
unknown  constants  in  the  perturbation  model  of  Eq.  (3.19),  representing 
drifts  and  zero  shifts,  to  allow  a  comparison  with  the  flight  responses: 


These  bias  and  reference  values  are  determined  from  a  least-squares  fit 
of  the  model  response  y(t)  and  the  low-pass  filtered  flight  data  yp(t) 
[with  x(0)  =0  to  avoid  overparameterizing  the  model].  The  more 
common  practice  of  simply  subtracting  the  initial  values  of  control  and 
output  from  the  time-histories  can  be  used  when  the  vehicle  dynamics  are 
stable,  since  drifts  and  reference  errors  decay  with  the  time-constant 
of  the  transfer-function  poles.  However,  when  the  identified  transfer 
function  is  highly  unstable,  as  it  is  in  the  pitch  and  roll  dynamics  of 
the  XV-15  in  hover,  small  errors  in  the  biases  and  reference  values 
quickly  propagate  and  the  responses  diverge  even  for  a  highly  accurate 
model,  ft  computer  program  (TFTHISTORY)  was  developed  by  the  author  to 
estimate  state-derivative  biases  and  output  reference  values  using  a 
least-squares  procedure  (a  simpler  but  still  satisfactory  alternative  to 
a  Kalman  filter  or  smoother  estimate),  generate  the  model  time-history, 
and  compare  the  result  with  the  filtered  flight  time-history  data. 

This  completes  the  discussion  of  the  frequency-domain  identifica¬ 


tion  method.  Identification  results  for  the  hover  and  cruise  flight 
conditions  are  presented  in  Chapters  4  and  5. 


Chapter  4 
HOVER  DYNAMICS 


This  Chapter  discusses  the  identification  and  analysis  of  the  XV-15 
open-loop  dynamics  in  hovering  flight.  The  transfer  functions  of  pri¬ 
mary  interest  for  handling-qualities  are  associated  with  "inner-loop" 
regulation  of  the  on-axis  responses:  pitch-rate  response  to  elevator 
(q/6e),  vertical-acceleration  response  to  power  lever  (az/6c),  yaw-rate 
response  to  rudder  (r/6_),  and  roll-rate  response  to  aileron  (p/6_). 

The  important  coupling  transfer  function  is  yaw-rate  response  to  aileron 
(r/6a).  To  illustrate  the  identification  method  described  in  Chapter  3, 

1  the  lateral/directional  on-axis  and  coupled  dynamics  are  discussed  in 

detail.  Flight-extracted  frequency-responses  are  compared  with  piloted- 
|  simulation  responses.  The  simulator  responses  are  from  the  model 

f 

states,  so  the  drive  dynamics  of  the  vertical  motion  simulator  (VMS) 
motion-base  system  are  not  included.  Comparisons  of  the  flight 
responses  and  the  current  nonreal-time  simulation  model  (Ref.  3)  are 
also  included  to  show  the  substantial  modeling  improvements  that  have 
been  made  in  the  3  years  since  the  VMS  simulation  tests  were  conducted. 
Closed-form  transfer-function  models  of  the  flight  responses  are 
obtained  using  least-squares  fits.  Finally  step-response  time-history 
matching  is  used  to  verify  the  predictive  capability  of  the  identified 
transfer  functions. 

The  results  presented  in  this  Chapter  are  considerably  improved 
over  the  original  hover  flight  condition  results  presented  in  Ref.  18. 
Subsequent  experience  on  the  XV-15  and  related  identification  efforts 


(Refs.  6,  21-24)  led  to  refinements  in  the  analysis  methods  and  insight 
into  techniques  for  extracting  higher-quality  results. 

4.1  Lateral/Directional  Dynamics 

The  frequency-sweeps  in  roll  were  the  most  difficult  to  execute 
because  of  significant  interaxis  coupling  between  lateral  stick  inputs 
and  yaw  response.  Since  the  yaw  SCAS  employs  a  washed-out  yaw-rate 
feedback,  low-frequency  yaw-attitude  excursions  required  continuous 
pilot  regulation  during  the  lateral  stick  frequency-sweeps.  However,  a 
two-input  analysis  of  the  roll-axis  response  showed  very  little  effect 
of  the  pedal  inputs  on  the  single-input  aileron  transfer  function.  This 
suggests  that  rudder  inputs  produce  largely  decoupled  yawing  motions; 
this  is  a  typical  assumption  for  hovering  aircraft  (Ref.  41).  There¬ 
fore,  the  following  analysis  of  the  roll  dynamics  considers  only  the 
single- input/single-output  transfer  function  p/6a . 

4.1.1  On-Axis  Responses 

The  lateral-stick  inputs  for  two  good  frequency-sweeps  were  concat¬ 
enated  to  form  the  composite  input  time-history  of  Fig.  4.1.  The  20-sec 

FIRST  SECOND 


Fig.  4.1.  Two  Lateral-Stick  Frequency-Sweeps  (Slat)  in  Hover> 


low-frequency  cycles  result  in  large  lateral-velocity  variations 
(10-20  knots).  At  the  higher  frequencies,  stick  deflections  are  larger, 
roll  rates  are  smaller,  and  velocity  perturbations  are  less  than 
5  knots.  Both  runs  are  roughly  90  sec  in  duration,  but  the  distribution 
of  frequency  content  in  each  is  seen  to  be  quite  different.  The  first 
frequency-sweep  contains  more  mid-frequency  power,  whereas  the  second 
contains  more  low-frequency  power.  The  input  autospectrum  for  the 
concatenated  time-history  is  shown  in  Fig.  4.2.  Maximum  input  power  is 


Fig.  4.2.  Lateral-Stick  Input  Autospectrum  {G*  A  ). 
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generated  at  a  frequency  of  u>  =  0.4  rad/sec,  which  is  slightly  higher 
than  that  corresponding  to  the  maximum  period  of  Tmax  =  20  sec.  The 
input  power  remains  strong  out  to  w  =  6  rad/sec,  and  then  falls  off 
rapidly  at  the  higher  frequencies.  The  frequency  content  which  exists 
below  the  minimum  even-input  frequency  (1/20  Hz)  and  above  the  maximum 
even-input  frequency  (4  Hz)  (not  shown  in  Fig.  4.2)  is  due  to  the 
various  nonsinusoidal  components  which  are  apparent  in  the  pilot  input 
(Fig.  4.1). 


The  corresponding  aileron  surface  deflection  time-history  is  shown 
in  Fig.  4.3.  This  signal  reflects  inputs  from  the  pilot  (and  shaped  by 
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Fig.  4.3.  Aileron  Surface  Deflections  (6„)  during  the  Lateral- 

a 

Stick  Frequency  Sweeps,  Reflecting  Contributions  From  Pilot  and  SCAS 
Inputs. 

SCAS  feed-forward  elements)  and  from  the  SCAS  feedbacks.  Once  again, 
the  first  frequency-sweep  has  more  mid-  and  high-frequency  power,  and 
the  second  frequency-sweep  has  more  low-frequency  power.  The  input 
autospectrum  for  the  concatenated  surface  deflection  signal  is  shown  in 
Fig.  4.4.  The  autospectrum  peaks  at  u>  =  0.3  rad/sec  (=  low-frequency 


Fig.  4.4.  Aileron  Surface  Input  Autospectrum  (Gg  g  ). 

ct  cl 

input)  and  is  roughly  constant  up  to  6  rad/sec.  The  fact  that  both  the 
lateral-stick  input  autospectrum  and  surface-deflection  input  auto¬ 
spectrum  are  strong  indicates  that  the  surface  deflection  contains  a 


significant  component  from  the  pilot's  stick  inputs  (thus  avoiding  the 
-1/Gf  problem  discussed  in  Chap.  3). 


The  roll-rate  responses  for  the  two  lateral  frequency-sweeps  are 
shown  in  Fig.  4.5.  The  roll-rate  amplitudes  vary  from  5-10  deg/sec, 
depending  on  input  frequency.  The  output  autospectrum  (Fig.  4.6)  peaks 


Fig.  4.5.  Roll-Rate  Response  (p)  during  Lateral-Stick 
Frequency-Sweeps . 


Fig.  4.6.  Roll-Rate  Output  Autospectrum  (Gpp). 


between  0.3  and  0.5  rad/sec,  because  of  the  low-frequency  input  cycles 
and  the  presence  of  the  dominant  roll  modes  in  this  frequency  range. 
Although  the  input  surface  autospectrum  remains  roughly  constant  for 
w  >  1  rad/sec,  the  output  spectrum  falls  off  rapidly  because  of  the 


rigid-body  (inertial)  response  of  the  aircraft.  The  cross-spectrum 
between  aileron  and  roll  rate  (Fig.  4.7)  is  quite  similar  in  form  to  the 
roll-rate  output  autospectrum  (Fig.  4.6)  because  the  aileron  input 
autospectrum  (Fig.  4.4)  is  roughly  constant  with  frequency. 


FREQUENCY,  rad/sec 

Fig.  4.7.  Magnitude  of  Cross-Spectrum  between  Aileron  Surface 


The  frequency  response  of  p/6a  is  shown  in  Fig.  4.8.  Once  again, 
the  magnitude  response  peaks  between  0.5  and  1  rad/sec,  a  result  of  the 
dominant  roll  modes  which  are  in  this  frequency  range.  For  frequencies 
greater  than  1  rad/sec,  the  Bode  plot  shows  the  classic  -20  dB/decade 
roll-off  caused  by  the  K/s  (inertial)  roll-rate  response  of  the  air¬ 
craft.  The  phase  increase  in  the  vicinity  of  the  magnitude  peak 
(w  =  0.3-0. 6  rad/sec)  indicates  that  the  dominant  roll  dynamics  are 
unstable.  At  higher  frequencies,  the  phase  is  roughly  -90  deg  because 
of  the  inertial  response  and  the  small  high-frequency  lags. 

The  coherence  function  for  the  roll-rate  response  identification  is 
3hown  in  Fig.  4.9.  Excellent  coherence  is  achieved  in  the  frequency 
ranges  of  0.2-0. 3  rad/sec  and  0.6-9  rad/sec.  The  coherence  function 
drops  sharply  in  the  frequency  range  in  which  the  input  and  output 
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Fig.  4.8.  Roll-Rate  Response  to  Aileron  (p/-«a).  (a)  Magnitude; 

(b)  Phase. 
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Fig.  4.9.  Coherence  Function  for  Roll-Rate  Response  Identification 


autospectrum  are  both  very  strong  and  the  transfer  function  peaks 
(0.3-0. 6  rad/sec).  This  indicates  that  nonlinear  behavior  is  present  in 
the  low-frequency  unstable  dynamics.  Current  research  by  the  author  and 
Kaletka  (Ref.  22)  shows  that  a  major  source  of  this  nonlinearity  is 
rolling-moment  dependency  on  higher  powers  of  lateral  velocity  [L(v^), 
L(v^)].  For  frequencies  below  0.2  rad/sec  and  above  9.0  rad/sec,  the 
coherence  function  drops  rapidly  because  of  the  lack  of  excitation 
outside  of  the  frequency  range  of  dominant  pilot  inputs. 

The  spectral  results  presented  thus  far  are  impressive  in  terms  of 
their  smoothness,  resolution,  and  coherence.  The  chirp  z-transform, 
described  earlier,  is  a  significant  improvement  over  the  standard 
(Cooley-Tukey )  FFT  methods,  which  were  found  to  produce  much  poorer 
frequency-response  resolution  and  associated  erratic  coherence-function 
behavior  at  low  frequency.  This  is  because  resolution  of  the  standard 
FFT  is  limited  by  the  length  of  the  window,  whereas  the  chirp 
z-transform  has  arbitrary  resolution  (Sec.  3.2.2).  In  addition,  the 
concatenation  of  good  repeat  runs  is  seen  to  produce  a  reasonably  con¬ 
stant  input  (surface)  autospectrum  over  the  frequency  range  of  interest, 
and  provides  a  large  number  of  sections  for  the  averaging  process.  This 
ensures  minimal  random  error  and  satisfactory  accuracy  of  the  transfer- 
function  estimates,  even  when  the  magnitude  of  the  coherence  function  is 
low  (for  example  near  w  =  0.4  rad/sec  in  Fig.  4.9). 

The  hover  frequency-sweep  tests  and  analysis  described  above  for 
the  flight  vehicle  were  repeated  in  the  VMS.  The  extracted  aileron-to- 
roll-rate  open-loop  transfer  function  is  shown  in  Fig.  4.8  for  compari¬ 
son  with  the  flight  results.  The  associated  coherence  function  shown  in 


Fig.  4.9  reflects  excellent  simulator  identification  for  the  entire 
range  of  0.1-10  rad/sec. 

The  correlation  of  the  VMS  simulation  and  the  flight  data  is  very 
good  for  the  frequency  range  1.0-10.0  rad/sec,  indicating  that  the  roll- 
control-response  sensitivity  is  accurately  modeled.  The  low-frequency 
phase  comparison  suggests  that  the  damping  ratio  of  the  unstable  roll 
mode  is  slightly  overestimated  by  the  simulation  model  (more  unstable), 
and  the  natural  frequency  is  accurate. 

An  interesting  demonstration  of  the  effects  of  nonlinearities  in 
the  low-frequency  roll  dynamics  (p/Sa)  is  presented  in  Fig.  4.10.  This 
figure  compares  the  simulator  frequency-sweep  result  with  the  small- 
perturbation  transfer  function  (obtained  by  numerically  linearizing  the 
VMS  simulation  model).  The  perturbation  transfer-function  matches  the 
sweep  results  well  in  the  high-frequency  range  ( oj  >  3.0  rad/sec)  where 
motions  are  small;  however,  considerable  deviation  occurs  in  the  low- 
frequency  end,  clearly  a  result  of  the  large-perturbation  nonlinear 
effects.  (The  large  phase  shift  at  low  frequency  is  due  to  low- 
frequency,  right-half  plane  poles  and  zeroes  in  the  simulator  perturba¬ 
tion  transfer-function.)  For  forward-flight  conditions,  the  differences 
were  found  to  be  much  less  significant  since  rotor  nonlinear i ties  are 
overshadowed  by  the  linear  aerodynamics  of  the  wing  and  tail,  and 
flight-path  excursions  about  the  steady-state  trim  condition  are 
smaller. 

The  fact  that  the  simulation  (and  flight  data)  coherence  function 
remains  strong  at  low  frequencies  (Fig.  4.9)  leads  to  the  important 
conclusion  that  the  dominant  nonlinear  effects  can  be  well  approximated 
by  linear  (but  not  small-perturbation)  describing  functions.  By  fitting 
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Fig.  4.10.  Comparison  of  Small-Perturbation  Transfer-Function  and 
Frequency-Sweep  Responses  Obtained  from  the  Simulation;  Data  Presented 
are  for  the  Roll-Rate  Response  to  Ailerons,  p/-6_.  (a)  Magnitude; 

d 

(b)  Phase. 

these  describing  functions  with  standard  transfer-function  forms,  the 
physical  significance  of  the  small-perturbation  model  parameters  is 
retained,  and  the  numerical  values  of  these  parameters  reflect  the  gross 
nonlinear  effects.  This  "equivalent  low-order  modeling"  approach  has 
been  useful  in  analyzing  and  predicting  the  handling  qualities  of  many 
aircraft  (Refs.  15,  17). 

The  low-frequency  errors  in  the  simulation,  which  are  shown 
in  Fig.  4.8,  were  also  reported  during  the  pilot's  qualitative 


H 
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evaluation.  Subsequent  analysis  of  the  real-time  mathematical  model 
indicated  problems  in  the  representation  of  rotor  flapping  for  large 


lateral-velocity  changes  (such  as  occurs  in  the  low-frequency  range  of 
the  frequency-sweep).  There  is  strong  sensitivity  of  the  numerically 
linearized  perturbation  transfer-functions  to  the  size  of  the  perturba¬ 
tion  increments.  Corrections  to  the  mathematical  modeling  in  these  two 
areas  were  included  later  in  the  nonreal-time  XV-15  simulation.  The 
comparison  between  the  current  nonreal-time  mathematical  model,  the  VMS 
original  model,  and  the  flight  data  is  shown  in  Fig.  4.8;  it  indicates 
the  significant  improvements  that  have  been  achieved  during  the  last  3 
years  of  development. 

The  identified  yaw-rate  response  to  rudder  inputs  r/6p  is  pre¬ 
sented  in  Fig.  4.11.  The  associated  coherence  function  (Fig.  4.12) 
shows  that  the  identification  is  accurate  from  0. 1-3.0  rad/sec;  above 


this  range,  the  erratic  coherence  behavior  suggests  a  less  accurate 
estimate. 

The  yaw-rate  response  follows  the  K/s  characteristic  in  magnitude 
and  phase  for  frequencies  greater  than  0.5  rad/sec,  which  means  that 
rudder  inputs  generate  constant  yaw-rate  accelerations.  The  fact  that 
the  phase  does  not  rise  above  the  -90  deg  value  except  at  very  low 
frequencies  suggests  an  associated  low  value  of  yaw  damping  (Np)  as 
compared  to  conventional  helicopters.  This  is  due  to  the  lack  of  a  tail 
rotor  for  the  tilt-rotor  configuration. 

A  comparison  of  the  simulator  yaw-rate  response  with  the  flight- 
test  data  is  presented  in  Figs.  4.11  and  4.12.  As  in  the  roll  axis,  the 
magnitude  response  generally  matches  very  well  for  higher  frequencies 
a)  >  1.0  rad/sec.  The  parallel  (vertical)  shift  in  the  magnitude  curve 
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Fig.  4.11.  Yaw-Rate  Response  to  Rudder  (r/6r).  (a)  Magnitude; 

(b)  Phase. 
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Fig.  4.12.  Coherence  Function  for  Yaw-Rate  Response  Identification 


indicates  that  the  control  sensitivity  is  overestimated  by  about  30J. 
Subsequent  analyses  as  a  result  of  this  finding  indicate  that  ballast, 
which  was  added  in  the  tail  (for  c.g.  control),  was  not  accounted  for  in 
the  simulation  model.  This  added  inertia  reduces  the  control  sensitiv¬ 
ity  and,  therefore,  could  contribute  to  the  observed  high-frequency 
discrepancy . 

The  yaw-mode  damping  is  significantly  underestimated  in  the  simula¬ 
tor,  as  can  be  seen  from  the  low-frequency  phase  drop  (which  actually 
indicates  an  unstable  response!)  compared  with  the  (stable)  flight 
response  (Fig.  4.11b).  This  discrepancy  was  verified  in  step-response 
comparisons  between  the  aircraft  and  the  simulator.  A  likely  cause  of 
error  in  the  yaw  damping  is  the  assumption  of  uniform  lateral  distribu¬ 
tion  of  the  rotor  inflow.  Wing-on-rotor  interference  on  the  inboard 
half  of  the  rotor  disks  probably  causes  a  significant  distortion  in  the 
shed  wake  and  consequently  in  the  lateral  inflow  distribution  in  hover¬ 
ing  flight.  The  assumption  of  uniform  inflow  causes  an  underestimation 
in  the  longitudinal  flapping  response  to  longitudinal  velocity  changes 
(Ref.  42),  which  is  the  primary  source  of  yaw  damping  for  the  hover 
flight  condition.  Even  minor  coupling  to  the  (unstable)  roll  dynamics 
could  then  result  in  the  observed  unstable  yaw  simulation  response. 

With  the  flight  and  simulation  frequency-responses  now  identified 
and  compared,  the  next  step  in  the  method  of  Chapter  3  is  to  derive 
lower-order  equivalent  transfer  functions.  Examination  of  the  roll  and 
yaw  frequency-responses  of  Figs.  4.8  and  4.11  shows  that  they  exhibit 
classic  characteristics  of  hovering  vehicles.  The  yaw  response  is 
decoupled  and  first-order  in  nature;  the  roll  response  is  dominated  at 
low-frequency  by  the  unstable  complex  mode  and  has  one  excess  pole  at 


67 


high  frequency.  Thus,  the  yaw-rate  response  to  rudder  inputs  is  modeled 
by 


r(s) 

Vs> 


N5  e 
r 


-t  ,s 
i|) 


(s  +  1/T  ) 


(4.1) 


where  r(s)/6p(s)  is  the  Laplace-transformed  yaw-rate  response  to 

rudder,  in  degrees  per  second  per  degree,  and  s  is  the  Laplace- 

transform  variable.  The  transfer-function  high-frequency  gain  (N.  )  is 

°r 

the  yaw-control  response  sensitivity,  that  is,  the  initial  yaw- 
acceleration  response  (deg/sec  )  to  a  1-deg  step  input  in  rudder.  The 
inverse  yaw-mode  time-constant  (1/Ty)  is  the  negative  value  of  damping 
stability  derivative  (Np).  The  effective  time-delay  (t^)  accounts  for 
phase  distortion  caused  by  unmodeled  high-frequency  dynamics.  The 
numerical  values  for  these  parameters  are  obtained  from  the  NAVFIT 
program  (Sec.  3-3.2)  using  a  fitting  range  of  0. 1-3.0  rad/sec: 


N,  =  0.619  deg/sec^/deg-rudder 
r 

=  5.45  deg/sec^/in.-ped 

(4.2) 

1/Ty  =  0.102  rad/sec 

=  0.0210  sec 

As  shown  in  Fig.  4.13,  this  model  fits  the  flight  data  well  at  mid- 


and  high-frequencies,  but  exhibits  noticeable  discrepancies  at  low 
frequencies,  where  the  simple  first-order  time-constant  approximation  is 


FREQUENCY,  rad/sec 

Fig.  4.13.  Transfer-Function  Model  Identification  for  Yaw-Rate 
Response  to  Rudder  (r/6r).  (a)  Magnitude;  (b)  Phase. 

inadequate.  The  low-frequency  mismatch  is  not  very  important,  as  is 
demonstrated  in  the  model  verification  study  (Sec.  4.3.1).  The  small 
effective  time-delay  (t^  =  21  msec)  indicates  that  the  unmodeled  high- 
frequency  dynamics  (e.g.,  rotor  and  actuators)  are  not  very  important  in 
the  identification  range  (0. 1-3.0  rad/sec). 


The  standard  roll-response  transfer-function  model  for  hovering 


vehicles  is  given  in  Ref.  41  as^ 


L.  s(  1/T  ) ( 1/T  )e 

p(s)  _  A  ^2 


'Tos 


6a(s) 


(1/Ty)(1/Tr)[cpf«r: 


(4.3) 


where  p(s)/5  (s)  is  the  Laplace-transformed  roll-rate  response  to 

ci 

aileron  surface  deflection,  deg/sec/deg-aileron ;  L.  is  the  aileron 

°a 

roll  sensitivity:  1/T,  ,  1/T.  are  first-order,  numerator  inverse  time- 
J  4>  i  4>2 

constants;  1/T  ,  is  the  first-order  roll  mode  inverse  time-constant;  cp, 
u>r  are  the  second-order  roll-mode  damping  and  natural  frequency, 
respectively;  and  is  the  effective  time-delay.  The  decoupled  yaw 

mode  1/T  is  fixed  at  its  previously  identified  [Eq.  (4.2)]  value 
(1/T  =  0.102  rad/sec).  The  remaining  numerator  and  denominator  parame¬ 

ters  are  varied  to  obtain  a  good  fit  in  the  range  of  0. 2-9.0  rad/sec. 

The  resulting  roll  transfer-function  parameters  are 


t  Shorthand  notation:  [c,  u>]  implies  s2  ♦  2cws  +  u>2,  5  =  damping  ratio; 
w  =  undamped  natural  frequency  (rad/sec);  and  (1/T)  implies  s  +  (1/T), 
rad/sec. 


=  -3.71  deg/secc/deg-aileron 

=  15.1  deg/sec2/in.-6Lafc 

1/T .  =  -0.107  rad/sec 

♦l 

1/T.  =  0.412  rad/sec 

$2 

1/Ty  =  0.102  rad/sec 

(4.4) 

1/Tr  =  1.23  rad/sec 
Cr  =  -0.418 
wr  =  0.447  rad/sec 

t.  =  0.0313  sec 

As  shown  in  Fig.  4.14,  the  transfer-function  model  is  generally  a  very 
good  representation  of  the  aircraft  roll  response.  However,  a  close 
examination  of  the  model-fit  and  flight  data  in  the  frequency-range  near 
the  dominant  mode  (w  =  0.5  rad/sec)  shows  the  inability  of  the  linear 
model  to  simultaneously  match  both  the  magnitude  and  phase  responses. 
Based  on  the  magnitude  response  alone,  a  somewhat  more  unstable  (higher 
negative  damping  ratio)  mode  is  indicated;  but,  based  on  the  phase 
response  alone,  a  slightly  less  unstable  (smaller  negative  damping 
ratio)  model  is  indicated.  This  inconsistency  reflects  the  significant 
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Fig.  4.14.  Transfer-Function  Model  Identification  for  Roll-Rate 
Response  to  Ailerons  (p/-6a).  (a)  Magnitude;  (b)  Phase. 

nonlinearities  which  are  in  this  frequency  range.  The  linear  model  of 
Eq.  (4.4)  is  thus  a  compromise,  and  reflects  the  current  relative 
weighting:  1  dB  magnitude  error:  7  deg  phase  error.  Increasing  the 

weighting  of  the  phase  errors  resulted  in  a  less  unstable  damping  ratio 
(cr),  better  phase  matching,  and  poorer  magnitude  matching  as 
expected.  The  current  weighting  was  retained  since  it  is  consistent 
with  practice  in  the  handling-qualities  community  (Ref.  37),  and  as 
shown  in  Sec.  4.3  it  leads  to  transfer-function  models  with  good  predic¬ 
tive  capability.  The  capability  to  specify  relative  magnitude  and  phase 
error  weighting  is  an  advantage  of  the  frequency-domain  approach. 
Time-domain  identification  methods  inherently  weight  the  phase  errors 
much  greater  than  the  magnitude  errors,  since  small  shifts  in  phase 


72 


between  the  flight-data  and  the  model-response  generate  large  changes  in 
the  time-domain  cost  function  (Ref.  22). 


Some  important  conclusions  can  be  drawn  from  the  transfer-function 
model  parameters  of  Eq.  (4.4),  and  the  location  of  the  poles  and  zeros 
in  the  complex  plane  (Fig.  4.15).  First,  the  open-loop  roll  response 


Fig.  4.15.  Locations  of  Roll-Response  Transfer-Function  Parameters 
in  the  s-Plane. 

dynamics  are  dominated  by  a  very  unstable  oscillation  in  hover,  with  an 
associated  time-to-double  amplitude  of  3.7  sec.  Reference  43  specifies 
an  unstable  damping  ratio  limit  of  z,  =  -0.22,  for  Level  I  handling 
qualities,  which  is  shown  on  Fig.  4.15.  The  fact  that  the  XV-15  roll 
mode  [-0.418,  0.447]  lies  beyond  this  boundary  suggests  that  the  SCAS- 
off  vehicle  has  Level  II  handling  qualities  (at  best).  This  is  cor¬ 
roborated  by  the  SCAS-off  handling-qualities  evaluation  discussed  in 
Sec.  2.1.3.  The  maximum  roll  acceleration  capability  of  1.27  rad/sec^ 
(based  on  the  lateral-stick  sensitivity  and  maximum  deflection)  provides 
a  sufficient  margin  for  SCAS-off  roll-attitude  control  (AGARD  specifica¬ 
tion,  Ref.  43). 
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The  second  important  conclusion  is  that  the  pole-zero  combination 


(1/T.  )/(1/Ty)  is  at  very  low  frequency  and  nearly  cancels  out 
(Fig.  4.15).  This  reveals  that  yaw  coupling  does  not  noticeably  affect 
the  roll-response  characteristics,  again  a  common  assumption  for  hover¬ 
ing  aircraft  (Ref.  41).  Third,  the  low-frequency  numerator  factor 

(1/T.  )  associated  with  coupled  lateral  translation/yaw  damping  is  also 

♦l 

the  limiting  location  of  the  dominant  lateral  translation  mode  when  the 
inner  attitude  is  tightly  pilot-regulated  (Ref.  44).  The  fact  that  this 
parameter  is  marginally  unstable  ( time-to-double  amplitude  =  6.8  sec) 
suggests  that  lateral  station-keeping  will  be  difficult.  This  is  typi¬ 
cal  of  unaugmented  helicopters  (Ref.  45)  and  is  corroborated  by  the 
handling-qualities  evaluation  (Sec.  2.1.3).  Fourth,  the  effective 
time-delay  (t^)  is  small,  again  suggesting  that  unmodeled  high-frequency 
dynamics  are  not  significant  in  the  frequency  range  of  identification 
(0. 2-9.0  rad/sec). 

4.1.2  Lateral/Directional  Coupling 

As  discussed  previously,  the  dominant  source  of  interaxis  coupling 
for  the  hovering  flight  condition  is  yaw-rate  response  to  aileron 
inputs .  Identification  of  this  coupling  characteristic  is  important  for 
assessing  the  pilot's  work  load  in  open-loop  lateral  station-keeping 
tasks . 

Clearly  the  pilot  and  SCAS  feedbacks  to  the  pedals  will  have  a 
first-order  effect  on  the  identification  of  the  open-loop  yaw-rate 
response  to  aileron.  Therefore,  the  application  of  the  two- input/ 
single-output  analysis  method  is  essential.  This  case  will  be  presented 
in  detail,  since  it  is  an  excellent  example  of  how  coupled  bar e-air frame 
characteristics  can  be  extracted  from  a  multi-loop  system. 


The  yaw-rate  response  to  the  lateral-stick  sweeps  of  Fig.  4.1  is 
shown  in  Fig.  4.16.  Excursions  generally  do  not  exceed  ±5  deg/sec. 
This  is  due  to  the  combined  effect  of  the  pilot  and  stability- 
augmentation  stabilization  of  the  yaw  degree-of-freedom  during  the 
lateral-stick  sweep.  In  essence,  the  pilot  is  attempting  to  eliminate 
the  yaw-rate  response  to  aileron— that  which  is  now  being  identified. 


Fig.  4.16.  Yaw-Rate  Response  (r)  during  Lateral-Stick  Frequency  Sweeps. 


The  autospectrum  for  aileron  surface  inputs  (Fig.  4.17)  is  very 
flat  across  the  entire  frequency  spectrum,  as  desired.  The  aileron 
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Fig.  4.17.  Aileron  Surface  Input  Autospectrum. 
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autospectrum  G.  (Fig.  4.17)  is  not  the  same  as  that  in  Fig.  4.4 

Va 

(although  the  time-history  data  are  identical)  because  the  optimum 
window  length  (Sec.  3.2.1)  was  found  to  be  longer  for  the  single-axis 
identification  of  p/«a  than  for  the  coupled-system  identification.  In 
this  case,  the  output  is  the  yaw-rate  signal,  which  has  an  autospectrum 
that  is  not  very  flat  (Fig.  4.18).  The  cross-spectrum  between  aileron 


FREQUENCY,  rad/sec 


Fig.  4.18.  Yaw-Rate  Output  Autospectrum. 


input  and  the  yaw-rate  response  (Fig.  4.19)  is  relatively  constant  in 
the  high-frequency  range,  and  peaks  at  low  frequency  (u»  =  0.3  rad/sec), 
as  does  the  output  autospectrum  (Fig.  4.18). 

The  frequency  response  of  aileron  to  yaw-rate,  with  rudder  inputs 
present  (r/6a:«r)  is  shown  in  Fig.  4.20.  The  frequency  response  is 
relatively  flat  in  the  mid-frequency  range,  indicating  a  fixed  relation 
between  yaw-rate  and  aileron  when  the  (pilot  and  SCAS)  yaw-axis  stabili 
zation  is  active.  The  ordinary  coherence  between  yaw-rate  and  aileron, 
shown  in  Fig.  4.21,  has  a  very  low  value  over  most  of  the  frequency 
range.  This  indicates  the  presence  of  other  inputs  which  significantly 
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Fig.  4.19.  Magnitude  of  Cross-Spectrum  between  Aileron  Surface 
Inputs  and  Yaw-Rate  Response. 
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Fig.  4.20.  Yaw-Rate  Response  to  Ailerons,  (a)  Magnitude; 


(b)  Phase. 
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Fig.  4.21.  Coherence  Function  for  Yaw-Rate  Response  to  Aileron 
2 

Inputs  ( Yg  p). 

°ar 

contribute  to  the  observed  yaw-rate  output  signal;  clearly,  these  are 
the  rudder  inputs,  which  are  being  supplied  by  the  pilot  and  the  yaw- 
stabilization  system. 

Rudder  surface  inputs  during  the  lateral  frequency-sweep  are  shown 
in  Fig.  4.22.  The  deflections  are  generally  within  ±5  deg,  with  most  of 
the  low-frequency  activity  coming  from  pilot  inputs.  The  coherence 


Fig.  4.22.  Rudder  Surface  Deflections  (6p)  during  the  Lateral- 


function  for  the  cross-cor relation  between  rudder  and  aileron  surface 


inputs  presented  in  Fig.  4.23  is  very  low  over  most  of  the  frequency 


Fig.  4.23.  Coherence  Function  for  Cross-Correlation  between 

2 

Aileron  and  Rudder  Surface  Inputs  (y,  ,  ). 

oo 

a  r 

range.  This  shows  that  the  inputs  are  not  fully  correlated,  the  key 
requirement  for  two-input/single-output  identification  (Sec.  3-2.5). 

For  other  flight  conditions  in  which  similar  analyses  were  conducted, 
the  cross-coherence  between  rudder  and  aileron  inputs  was  generally  much 
greater  because  pilot-stabilization  requirements  in  the  secondary  axis 
were  less  severe.  This  increased  cross-coherence  caused  an  overall 
degradation  in  the  cross-coupling  identification  for  the  forward-flight 
cases.  Therefore,  when  possible,  the  frequency-sweeps  should  be  con¬ 
ducted  without  engaging  the  stability-augmentation  systems  in  the  sec¬ 
ondary  degrees-of- freedom. 

The  conditioned  aileron  input  autospectrum  (G.  ),  that  is, 

6aV6r 

with  the  effect  of  rudder  inputs  eliminated,  is  presented  in 
Fig.  4.17.  The  comparison  of  results  in  this  figure  shows  that  the 
input  autospectrum  is  not  substantially  changed.  This  is  because  the 
cross-coherence  between  inputs  is  small.  A  similar  characteristic  is 
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apparent  for  the  conditioned  output  autospectrum  (G  )  shown  in 

rr*6r 

Fig.  4.18.  The  conditioned  cross-spectrum  between  aileron  and  yaw-rate 

(Gfi  )  presented  in  Fig.  4.19  is  generally  more  flat  than  before,  a 
a  r 

result  of  the  elimination  of  the  yaw-rate  feedback  effect. 

The  conditioned  frequency-response  of  aileron  to  yaw-rate  (r/6  *6  ) 

3  r 

is  shown  in  Fig.  4.20.  Notice  that  the  shape  of  the  magnitude  and  phase 
responses  is  significantly  different  than  it  was  before  the  rudder 
inputs  were  taken  into  account.  The  response  from  1.0-10.0  rad/see  now 
exhibits  the  familiar  -20  dB/decade  roll-off.  This  suggests  that  the 
interaxis  coupling  is  dominated  by  the  aileron-to-yawing  moment 
control-coupling  derivative  (Nfi  ),  as  expected.  The  low-frequency 

H 

dynamics  are  unstable,  as  seen  earlier  in  the  roll  response. 

2 

The  partial  coherence  function  y.  „  of  Fig.  4.24  shows  that 

6  r  •  6 
a  r 

satisfactory  identification  of  the  coupling  dynamics  has  been  achieved 


Fig.  4.24.  Partial  Coherence  Function  for  Yaw-Rate  Response  to 
2 

Aileron  Inputs  (y{  p  &  ). 

a  r 

in  the  frequency  ranges  of  0. 1-0.4  rad/sec  and  0.75-7.0  rad/sec.  The 

p 

multiple  coherence  plot  y  of  Fig.  4.25  indicates  that  the  com- 

6_ • r 
a  r 

bined  linear  effects  of  aileron  and  rudder  account  for  most  of  the 
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Fig.  4.25.  Multiple  Coherence  Function  for  Yaw-Rate  Response  to 

2 

Combined  Aileron  and  Rudder  Surface  Inputs  (y5  >r). 

a  r* 

yaw-rate  response  at  high  frequencies  with  a  drop-off  at  lower  frequen¬ 
cies  because  of  nonlinearities  and  turbulence,  as  was  seen  previously  in 
the  single-input,  roll-response  analysis. 

The  following  transfer-function  model  adequately  characterizes  the 
aileron-to-yaw  rate  coupling: 


jin 

6a(s)  -  (1/Ty)(1/Tr)tcr,u»r] 


(4.5) 


where  r(s)/6a(s)  is  the  Laplace- transformed  yaw-rate  response  to 
aileron  surface  deflection,  deg/sec/deg-aileron;  Nfi  is  the  aileron  yaw 
sensitivity;  1/T^  is  the  first-order  numerator  inverse  time-constant; 

and  are  the  second -order  numerator  damping  and  natural  frequency, 
respectively;  and  is  the  effective  time-delay.  The  denominator 

factors  are  the  previously  identified  lateral  modes  [Eq.  (4.4)].  These 


denominator  factors  are  held  constant,  and  the  remaining  parameters  are 


iteratively  varied  using  NAVFIT  to  fit  the  data  in  the  frequency  range 
0. 1-7.0  rad/sec.  The  resulting  transfer-function  parameters  are 

Ng  =  0.344  deg/sec2/deg-aileron 

=  -1.40  deg/s2/in.-6Lat 

1/T^  =  -0.345  rad/sec 

=  0.868 

=  0.487  rad/sec  (4.6) 

1/T  =  0.102  rad/sec 

1 /Tr  =  1.23  rad/sec 

Cr  =  -0.418 

wr  =  0.447  rad/sec 

t.  =  0.00900  sec 
*  2 

The  yaw-rate  responses  to  both  aileron  and  rudder  are  dominated  by 
the  "K/s  characteristic"  in  the  pilot's  operating  (crossover)  frequency 
region  for  attitude  control  [generally  taken  as  1. 0-3.0  rad/sec 
(Ref.  41)].  Therefore,  the  pilot  cross-feed  (G^p)  that  is  needed  to 


maintain  a  fixed  heading  during  lateral  manuevering  is  approximately  a 
constant : 


CCF  =  ^  =  0.556  deg-rud/deg-ail  =  0.489  *-ped/%-lat  (4.7) 
6 

r 

This  shows  that  pedal  inputs  must  be  about  50%  of  the  lateral  stick 
inputs  to  maintain  a  constant  heading  in  the  piloted  crossover  frequency 
region.  This  level  of  coupling  is  quite  severe,  and  probably  contrib¬ 
utes  to  the  poor  SCAS-off  handling-qualities  in  lateral  translation 
tasks  (Sec.  2.1.3). 

The  comparison  between  the  transfer-function  model  fit  and  the 
flight- test  results  for  this  case  is  shown  in  Fig.  4.26.  The  magnitude 
response  correlates  well  with  somewhat  poorer  agreement  for  the  phase 
response.  However,  the  nature  of  the  coupling  characteristics  is  well 
represented. 

In  the  simulated  lateral-stick  sweeps,  pilot  regulation  of  yaw 

excursions  was  not  very  aggressive.  Therefore,  the  rudder  inputs  were 

dominated  by  the  SCAS,  causing  a  strong  correlation  between  the  rudder 

and  aileron  responses.  This  precluded  the  application  of  the  two-input 

method.  Instead,  the  flight-test  results  are  compared  in  Fig.  4.26  with 

the  real-time  mathematical  model  (numerical)  perturbation  transfer- 

function.  As  expected,  the  high-frequency  correlation  is  excellent, 

since  it  depends  only  on  the  cross-coupling  moment  derivative  N.  .  The 

°a 

agreement  at  low  frequency  is  much  less  satisfactory  because  nonlinear 
yaw  and  roll  dynamics  are  not  well  represented  by  the  small  perturbation 
transfer-function  (Sec.  4.1.1). 
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Fig.  4.26.  Comparison  of  Flight-Test  and  Simulator  Results  for 
Yaw-Rate  Response  to  Aileron  (r/6  ).  (a)  Magnitude;  (b)  Phase. 


4.2  Longitudinal  Dynamics 

Power  lever  inputs  produce  essentially  decoupled  first-order 
vertical-acceleration  (heave)  responses.  So,  these  can  be  considered 
independently  from  the  pitch  responses.  The  vertical  acceleration 
response  to  power-lever  inputs  (az/6c)  is  shown  in  Fig.  4.27.  For 
low-frequency  inputs,  the  first-order  heave  mode  dominates  the  response, 


-  FLIGHT  DATA 

. SIMULATOR  DATA 

- MODEL 


Fig.  4.27.  Vertical-Acceleration  Response  to  Power  Lever 
(-az/6c).  (a)  Magnitude;  (b)  Phase. 

causing  the  phase  rise  seen  in  Fig.  4.27b.  The  relatively  slight  change 
in  the  magnitude  response  (Fig.  4.27a)  suggests  that  this  mode  is  at  a 
very  low  frequency,  as  is  verified  later.  In  the  frequency  range  of 
1.0-10.0  rad/sec,  the  phase  curve  is  nearly  constant,  so  the  aircraft 
vertical-acceleration  response  to  power-lever  inputs  is  instantaneous, 
and  the  high-frequency  engine-governor  dynamics  are  not  significant. 

The  coherence  function  is  shown  in  Fig.  4.28,  and  indicates  an  excellent 
identification  of  the  vertical  response  characteristics  in  the  frequency 
range  of  0. 1-3.0  rad/sec. 

The  simulation-response  characteristics  are  also  shown  in 
Figs.  4.27  and  4.28.  In  general,  the  agreement  is  good  in  the  form 
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Fig.  4.28.  Coherence  Function  for  Vertical-Acceleration  Response 
2 

Identification  (y  ). 

o  a„ 
c  z 

of  the  amplitude  and  phase  characteristics  for  frequencies 
a)  >  0.3  rad/sec,  suggesting  that  the  modeling  of  initial  thrust  response 
to  collective  inputs,  and  governor  dynamics,  is  satisfactory .  The 
constant  amplitude  error  (vertical  shift)  may  be  due  to  a  mismatch  in 
the  assumed  vehicle  mass  or  rotor  lift  curve  slope.  The  low-frequency 
agreement  is  not  as  good,  which  appears  to  be  due  to  an  underestimate  of 
the  heave-mode  time  constant  by  the  simulation. 

The  vertical-velocity  response  (w  =  a2/s)  to  collective  is  identi¬ 
cal  in  form  to  the  yaw-rate  response  to  rudder  presented  earlier.  In 
fact,  these  two  decoupled  degrees-of-freedom  are  directly  analogous. 
Typically,  they  both  have  very  long  time-constants  for  tilt-rotor-class 
VTOL  vehicles,  in  contrast  to  single-rotor  helicopters,  so  that  their 
response  is  basically  an  acceleration;  the  control  sensitivity  is  a 
function  of  the  applied  control  force  (or  moment)  normalized  by  the 
respective  mass  characteristic.  A  first-order  model  fit  to  the 
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vertical-acceleration  response  in  the  frequency  range  0. 1-3.0  rad/sec 
yields 


a  (s) 

= 


-0.00980s  e 


-0.00740s 


(s  +  0.105) 


g/% 


(4.8) 


where  the  heave-mode  time-constant  has  about  the  same  value  as  the 
yaw-mode  time-constant  presented  earlier,  and  the  time  delay  is  again 
small.  The  control  sensitivity  (-0.0098  g/%)  is  equivalent  to 
-0.1  g/in.  of  power  lever;  this  is  adequate,  based  on  the  data  of 
Ref.  43.  The  very  low  value  of  the  flight-extracted  inverse  time- 
constant  (1/Th  =  0.105  rad/sec)  relative  to  typical  hovering  helicopters 
is  due  to  the  high  disc-loading  of  the  tilt-rotor  configuration. 

Although  the  frequency-sweeps  in  roll  were  more  difficult  to  per¬ 
form  than  those  in  pitch,  the  identification  of  the  pitch  results  proved 
to  be  the  most  difficult.  Stick  inputs  tended  to  be  very  small  at  low 
frequency,  which  resulted  in  poor  uncorrelated  input  power.  Of  the 
three  longitudinal  frequency-sweeps  which  were  conducted  (and  used  in 
the  original  analysis  of  Ref.  18),  only  one  run  subsequently  proved  to 
be  of  satisfactory  quality  (based  on  input  autospectrum  and  coherence) 
for  the  pitch-response  identification  in  the  mid-  and  low-frequency 
range.  In  the  high  frequency  range  (w  >  2  rad /sec),  all  sweeps  are  of 
good  spectral  quality,  so  concatenation  of  the  time-histories  is  useful 
to  obtain  minimum  variance  results  there.  However,  to  achieve  good 
identification  over  the  entire  frequency  range  only  the  single 
frequency-sweep  is  used.  The  possibility  of  combining  spectral  results 
which  are  optimiaed  for  particular  frequency-ranges  has  not  been  fully 


investigated,  but  may  enhance  the  overall  quality  of  the  identification 
results. 

The  elevator  to  pitch-rate  frequency  response  is  shown  in 
Fig.  4.29.  As  before  in  the  roll  axis,  the  response  peak  is  due  to  the 
dominant  pitch  modes  which  are  at  low  frequency.  For  frequencies 
ai  >  2  rad/sec,  the  magnitude  response  rolls  off  with  the  K/s  rigid- 
body  characteristic,  as  expected.  The  phase  rise  near  the  low-frequency 
response  peak  reflects  an  unstable  pitch-mode  associated  with  the 
speed/attitude  divergence,  as  in  the  roll  case.  At  high-frequency,  the 
pitch-rate  and  elevator  traces  are  just  over  90  deg  out  of  phase  because 
of  the  rigid-body  response  and  the  small  high-frequency  lags. 
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Fig.  4.29.  Pitch-Rate  Response  to  Elevator  (q/-6g). 


(a)  Magnitude;  (b)  Phase. 


The  coherence  function  shown  in  Fig.  4.30  indicates  a  strong  iden¬ 
tification  in  the  frequency  range  0. 2-6.0  rad/sec.  The  drop  in  coher¬ 
ence  between  0. 5-1.0  rad/sec  suggests  that  nonlinearities  are  important 
in  low-frequency  dynamics,  as  would  be  expected  from  the  similar  roll- 
response  behavior. 
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Fig.  4.30.  Coherence  Function  for  Pitch-Rate  Response 
2 

Identification  (yr  ). 
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The  simulation  pitch-response  characteristics  are  also  shown  in 

Figs.  4.29  and  4.30.  The  parallel  (vertical)  shift  in  the  magnitude 

plots  at  high  frequency  suggests  differences  in  the  high-frequency  gain 

parameter,  .  As  in  the  yaw-case,  the  discrepancy  (about  25%)  is 
°e 

probably  a  result  of  ballast  in  the  tail  of  the  aircraft,  which  was  not 
included  in  the  simulation  model.  Considerable  differences  between  the 
simulation  and  flight-extracted  responses  are  apparent  in  the  frequency 
range  of  0.3-0. 8  rad/sec.  The  simulation  model  appears  to  underestimate 
both  the  frequency  and  the  unstable  damping  ratio  of  this  mode.  Com¬ 
parisons  of  time-domain  identification  results  (Ref.  22)  with  the 


nonreal-time  simulation  indicates  that  "speed-stability  derivative" 


^  is  underestimated  in  the  simulator  by  about  70% 

^^u^FLIGHT  =  °'00^22  versus  (Mu)siMUl  ATION  =  °*00200  rad/sec2/ft/sec] 
which  would  contribute  to  the  observed  frequency  and  damping 
discrepancies  (Ref.  42).  A  probable  source  of  this  error  in  the  Mu 
derivative  is  the  rotor  longitudinal  flapping  angle  predictions  based  on 
the  uniform  lateral  inflow  approximation,  as  was  also  indicated  in  the 
yaw-rate  response  (Sec.  4.1.1). 

The  transfer-function  model  for  pitch-rate  response  to  elevator  has 
the  same  form  as  the  roll-response  model  of  Eq.  (4.3): 


Vs(1/T.  )( 1/T_  )S 
q(s)  .  6  91  62 

«e(s)  -  (1/Th)(1/Tp)[tp,u,p] 


(4.9) 


where  q(s)/6e(s)  is  the  Laplace- transformed  pitch-rate  response  to 

elevator  surface  deflection,  deg/sec/deg-elevator ;  M.  is  the  elevator 

°e 

pitch  sensitivity;  1/T0^,  1/T0^  are  first-order  numerator  inverse 
time-constants;  1/Tp  is  first-order  pitch-mode  inverse  time-constant; 
Cp  and  iup  are  the  second-order  pitch-mode  damping  and  natural  fre¬ 
quency,  respectively;  and  t0  is  the  effective  time-delay.  As  before, 
the  decoupled  mode,  (in  this  case  1/Th),  is  fixed  at  its  previously 
identified  value  (1/Th  =  0.105  rad/sec).  The  remaining  numerator  and 
denominator  parameters  are  varied  to  obtain  a  good  fit  in  the  frequency 
range  0. 2-5.0  rad/sec: 


f  * 

V 


a 


M.  =  -2.66  deg/sec^/deg-elevator 
°e 

=  -13. 2  deg/s€c2/in.-6Lon 

1/T„  =  -0.271  rad/sec 
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1/T0  =  0.508  rad/sec 

e2 

1/Th  =  0.105  rad/sec 

(4.10) 

1/Tp  =  1.32  rad/sec 
tp  =  -0.463 
u>p  =  0.579  rad/sec 
t q  =  0.0656  sec 

A  comparison  of  the  lower-order  model  with  the  flight-test  results  is 
shown  in  Fig.  4.29.  Overall,  the  model  is  a  good  approximation  to  the 
flight  data  in  the  fitting  range.  However,  as  in  the  roll  case,  the 
presence  of  nonlinearities  in  the  low-frequency  unstable  dynamics  pre¬ 
cludes  an  entirely  representative  fit  in  both  magnitude  and  phase  near 
the  response  peak. 

Once  again,  some  important  conclusions  can  be  drawn  from  the 
lower-order  transfer-function  model.  First,  the  open-loop  (bare- 
airframe)  configuration  is  very  unstable  in  the  pitch  axis,  with  nearly 
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the  sane  damping  ratio  and  natural  frequency  as  in  the  roll  axis.  As 
before,  poor  SCAS-off  handling-qualities  are  indicated  and  have  been 
reported  in  piloted  evaluations  (Sec.  2.1.3).  Second,  the  pole-zero 


pair  associated  with  heave  motions  ( 1/T0^)/( 1/Th)  is  widely  separated, 
implying  that  significant  heave-pitch  coupling  results  from  pitch- 
control  inputs.  This  separated  pole-zero  pair  precludes  the  use  of  the 


standard  third-order,  pitch  transfer-function  model  for  hover.  Third, 


the  numerator  factor  (1/T„  )  associated  with  the  translational  mode 

e1 

dynamics  for  constrained  attitude  control  (Ref.  44),  is  very  unstable 
(time-to-double  amplitude  =  2.6  sec),  and  probably  contributes  to  the 
poor  translational-flight  handling-qualities  (Sec.  2.1.3).  Fourth,  the 
effective  time-delay  (t0)  is  very  small,  as  it  is  in  the  roll  axis, 
suggesting  that  unmodeled  high-frequency  dynamics  are  not  significant  in 
the  frequency  range  of  identification  (0. 2-5.0  rad/sec). 


4.3  Verification  of  Open-Loop  Transfer-Function  Models  for  Hover 

Verification  of  all  identified  open-loop  transfer  functions  should 
be  completed  with  open-loop  flight- test  data.  However,  the  duration  and 
size  of  the  verification  test-inputs  is  severely  restricted  by  the 
inherent  stability  and  response  characteristics  for  each  axis  and  flight 
condition.  The  decoupled  stable  first-order  nature  of  the  open-loop  yaw 
and  heave  responses  allows  large  and  very  long  duration  step-inputs 
(about  10  sec)  to  be  executed.  Then,  the  modes  with  frequencies  of 
oi  >  0.1  rad/sec  contribute  at  least  one  time-constant  of  their  inherent 
response  to  the  overall  behavior,  which  checks  the  yaw  and  heave  trans¬ 
fer-functions  [Eq.  (4.2)  and  (4.8),  respectively]  over  the  entire  iden¬ 
tified  frequency  range.  However,  as  previously  discussed,  the  open-loop 


pitch  and  roll  dynamics  are  dominated  by  strong  instabilities,  with 
time-to-double  amplitudes  of  about  3  sec.  This  severely  restricts  the 
allowable  size  and  duration  of  the  verification  test  inputs  in  these 
axes. 

To  reduce  the  level  of  coupled  divergence,  verification  of  the 
open-loop  pitch  response  is  completed  with  the  roll  and  yaw  SCAS  axes 
engaged.  Similarly,  verification  of  the  open- loop  roll  response  is 
completed  with  the  pitch  and  yaw  SCAS  axes  engaged.  Even  so,  maximum 
achieveable  step-input  durations  are  1-2  sec,  which  allows  the  dominant 
modes  with  frequencies  of  about  u>  =  0.5  rad/sec  to  experience  only 
one  time  constant.  Therefore,  the  verification  step-inputs  do  not  give 
a  complete  test  of  model  fidelity  at  very  low  frequencies.  However, 
ensuring  model  accuracy  at  mid-  and  high-frequency  is  of  principal 
importance  for  handling-qualities  (and  control)  applications  since 
the  pilot  is  mostly  concerned  with  the  response  in  the  first  few 
seconds.  SCAS-on  step  responses  are  often  useful  for  verification  when 
quiescent  initial  test  conditions  cannot  be  achieved  in  the  open-loop 
configuration. 

4.3.1  Roll  Response  to  Aileron  Inputs 

The  (filtered)  roll-rate  response  of  the  aircraft  to  a  step  aileron 
input  with  the  roll  SCAS  disengaged  is  shown  in  Fig.  4.31.  The  input  is 
held  for  2  sec  and  the  peak  response  is  10  deg/sec,  which  is  roughly  the 
same  as  was  achieved  during  the  lateral  frequency-sweeps  (Fig.  4.5). 

When  the  transfer-function  model  of  Eq.  (4.4)  is  driven  with  the  same 
(filtered)  input  flight  data,  the  response  co-plotted  in  the  dashed 
curve  of  Fig.  4.31  is  obtained.  The  model  and  flight  data  responses 
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Fig.  4.31.  Roll-Response  Model  Verification  Using  SCAS-off  Flight 
Data,  (a)  Aileron  Input;  (b)  Roll-Rate. 


generally  compare  very  well,  although  there  is  a  slight  underestimation 
of  the  peak  roll  rate,  probably  a  result  of  nonlinear  effects. 

4.3.2  Yaw  Response  to  Rudder  Inputs 

The  yaw-rate  responses  of  the  aircraft  and  model  Eq.  (4.2)  to  an 
open-loop  step-rudder  input  in  hover  are  shown  in  Fig.  4.3 2.  The  con¬ 
trol  deflection  represents  about  40>  of  the  maximum  control  power,  a 
fairly  large  amplitude  input.  The  maximum  yaw-rate  response  is  about 
35  deg/sec,  which  is  40£  greater  than  the  peak  values  obtained  in  the 
frequency  sweeps.  The  model  and  flight  responses  generally  compare  very 
well  over  the  entire  run — even  for  this  fairly  large  input  amplitude 
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Fig.  4.32.  Yaw-Response  Model  Verification  Using  SCAS-off  Flight 
Data,  (a)  Rudder  Input;  (b)  Yaw-Rate. 


case.  The  small  discrepancies  are  probably  associated  with  the  mid¬ 
frequency  mismatch  between  the  flight-extracted  frequency  response  and 
transfer-function  fit  (Fig.  4.13),  or  large  amplitude  nonlinear  effects. 

4.3.3  Yaw  Response  to  Aileron  Inputs 

The  verification  of  the  identified  model  for  yaw-rate  response  to 
aileron  (Eq.  (4.6)]  is  shown  in  Fig.  4.32.  The  SCAS  configuration  is 
roll,  SCAS-off;  and  yaw,  SCAS-on.  Therefore,  some  rudder  activity  is 
also  present  (Fig.  4.33b).  The  roll-rate  match  is  very  good 
(Fig.  4.33c)  as  was  seen  in  the  previous  comparison  (Fig.  4.31).  The 
observed  yaw-rate  (Fig.  4.33d)  is  composed  of  responses  to  rudder  and 
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aileron  inputs  in  about  equal  proportions.  The  total  yaw-rate  response 
has  the  proper  form  although  the  peak  is  underestimated;  this  is  prob¬ 
ably  a  result  of  nonlinearities  in  the  coupled  response  characteristics 
(or  a  gust),  as  well  as  the  slight  underestimation  in  the  yaw-rate 
response  to  rudder  (Fig.  4.32). 

4.3.4  Vertical  Acceleration  Response  to  Collective 
The  verification  data  for  the  step  inputs  in  the  collective  axes 
are  very  noisy.  Even  so,  the  correlation  of  the  model  and  flight  data 
(Fig.  4.34)  indicates  that  a  proper  identification  has  been  achieved. 
The  comparison  is  better  at  the  beginning  of  the  run  where  the  inputs 
are  smooth  as  compared  to  the  end  of  the  run  where  the  inputs  are  more 

- FLIGHT  DATA 

- MODEL 


Fig.  4.34.  Vertical-Acceleration  Response  Model  Verification, 
(a)  Power  Lever  Input;  (b)  Vertical  Acceleration. 
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sharp.  This  is  because  the  sharp  change  in  collective  (rotor  thrust) 
excites  the  rotor  inflow  dynamics  which  are  not  included  in  the  simple 
first-order  model  (Ref.  6). 

**•3.5  Pitch-Rate  Response  to  Elevator 

The  aircraft  and  model  responses  to  a  2-sec  step-input  in  elevator 
with  the  pitch  SCAS  disengaged  are  shown  in  Figs.  4.35.  The  relatively 


FLIGHT  DATA 


Fig.  4.35.  Pitch-Response  Model  Verification  Using  SCAS-off  Flight 
Data,  (a)  Elevator  Input;  (b)  Pitch-Rate. 

unsteady  initial  conditions  cause  the  shift  between  the  flight  data  and 
transfer-function  models  at  the  beginning  of  the  run  (t  =  6-8  sec),  but 
the  comparison  improves  considerably  during  the  large  amplitude  recovery 
phase  (t  =  8-15  sec).  Somewhat  better  initial  conditions  are  achieved 
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with  the  pitch-SCAS  engaged.  The  elevator  surface-deflection  response 
to  a  step  input  of  longitudinal  stick  is  a  doublet,  as  seen  in 
Fig.  4.36a.  As  expected  from  the  control  system  configuration 
(Chap.  2),  a  steady  pitch-rate  is  achieved  within  0.5  sec  of  the  longi¬ 
tudinal  stick  input.  The  model  prediction  is  seen  to  be  very  good  for 
the  input  and  recovery  (t  =  9-13  sec). 
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Fig.  4.36.  Pitch-Response  Model  Verification  Using  SCAS-on  Flight 
Data,  (a)  Elevator  Input;  (b)  Pitch-Rate. 


The  verification  results  confirm  the  following: 

1.  The  hover  dynamics  are  essentially  linear. 

2.  Mismatches  between  the  extracted  frequency-responses  and 
transfer-function  model  fits  are  not  significant,  except  perhaps  in  the 
yaw-rate  response  to  rudder. 

3.  The  derived  transfer-function  models  are  not  overtuned  or 
over-parameterized . 

4.  The  unbiased  open-loop  vehicle  response  has  been  identified 
from  the  closed-loop  flight  data. 

The  derived  transfer-function  models  for  hover  are  summarized  in 
Table  4.1.  The  verification  results  show  the  utility  of  these  fairly 
simple  models  in  representing  the  transient  responses  to  relatively 
large  and  varied  inputs.  Therefore,  the  models  presented  in  Table  4.1 
should  be  used  for  future  XV-15  control  system  and  handling-qualities 
analyses  for  the  hover  flight  condition. 

The  next  chapter  discusses  the  dynamics  for  the  cruise  flight 
condition. 


Chapter  5 


CRUISE  DYNAMICS 

This  chapter  discusses  the  open-loop  dynamics  of  the  XV- 15  for  the 
cruise  flight  condition  (flight  condition  4,  Sec.  2.2). 

As  in  hovering  flight,  the  longitudinal  and  lateral  flight  dynamics 
are  essentially  fully  decoupled.  The  primary  longitudinal  bare-airframe 
transfer  functions  of  interest  for  "inner-loop"  regulation  are  pitch- 
rate  and  normal-acceleration  (at  the  center-of-gravity )  responses  to 
elevator  deflection  ( q/6e )  and  (az/6e)  respectively.  The  important 
lateral-directional  transfer  functions  are  roll-rate  response  to  aileron 
( p/S  )  and  sideslip  (at  the  center-of-gravity)  response  to  rudder 

cl 
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Frequency-sweep  tests  were  executed  much  more  easily  in  the  cruise 
flight  condition  than  in  the  hover  condition.  Cruise  dynamics  are 
inherently  very  stable,  and  transient  motion  is  small  relative  to  the 
large  steady  airspeed  (Vi  =  170  knots).  This  allowed  good  flight  data 
to  be  rapidly  obtained  without  requiring  a  large  number  of  repeti¬ 
tions.  Also,  high-quality  spectral  identification  was  achieved  without 
considerable  adjusting  of  the  various  processing  options — as  was  neces¬ 
sary  in  hover.  Only  two  frequency-sweeps  were  needed  in  each  axis 
because  the  spectral  quality  of  all  of  the  individual  runs  was  satisfac¬ 
tory,  and  none  needed  to  be  discarded  (as  was  required  in  the  hover 
case).  The  chirp  z-transform  algorithm  yielded  excellent  frequency- 
response  identification  in  both  axes  over  a  wide  frequency  range. 

Frequency-sweep  data  from  the  (VMS)  piloted  simulation  were  not 
available  for  the  (V^  =  170  knots)  cruise  flight  condition,  so  a  direct 


comparison  was  not  possible.  However,  the  nonreal-time  mathematical 
model  (Ref.  3)  for  this  condition  is  essentially  identical  to  the  real¬ 
time  piloted-simulation  (Ref.  2).  Since  the  relative  amplitude  of 
motion  about  the  high  equilibrium  flight  speed  is  generally  small,  and 
the  dynamics  in  the  fixed-wing  configuration  are  basically  linear,  the 
perturbation  transfer-functions  obtained  from  the  nonreal-time  version 
are  satisfactory  for  mathematical  model  validation  in  this  flight 
condition. 


5.1  Lateral/Directional  Dynamics 

5.1.1  Frequency-Response  Identification 

The  lateral  stick  inputs  for  two  concatenated  frequency-sweeps  in 
the  hover  flight  condition  are  shown  in  Fig.  5.1a.  The  input  form  is 
seen  to  be  much  more  regular  than  in  hover  (compare  with  Fig.  4.1).  The 
two  runs  are  quite  similar,  with  maximum  amplitudes  of  about  20%  of  the 
full  control  authority.  Since  these  tests  were  conducted  with  the 
lateral/directional  SCAS  axes  disengaged,  the  surface  deflection  (exci¬ 
tation  input)  is  related  to  the  cockpit  signal  by  a  constant  (mechanical 
compliance  and  hydraulic  actuators  contribute  less  than  10  msec  of 
effective  time  delay).  The  roll-rate  response  for  these  two  frequency- 
sweeps  is  shown  in  Fig.  5.1b.  The  maximum  response  amplitude  of 
10  deg/sec  is  similar  to  that  experienced  in  hover  (Fig.  4.5). 
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Fig.  5.1.  Two  Lateral  Stick  Frequency-Sweeps  (6LftT)  in  Cruise, 
(a)  Lateral  Stick  Inputs;  (b)  Roll-Rate. 


The  (open-loop)  roll-rate  response  to  ailerons,  p/6_,  is  shown  in 

a 

Fig.  5.2.  This  is  a  classical  first-order  response  in  which  aileron 
inputs  produce  a  constant  roll  acceleration  at  high  frequency  and  a 
constant  roll  rate  at  low  frequency.  The  corner  frequency  is  at  about 
id  =  1  rad/sec,  with  an  associated  phase  lag  of  roughly  *  =  -^5  deg,  as 
expected.  The  dominant  time-constant  is  thus  about  1  sec,  implying 
about  2-3  sec  to  reach  a  steady-state  roll  rate.  Higher-order  dynamics 
are  not  significant  in  this  response  since  the  phase  lag  is  nearly 
constant  for  frequencies  above  3  rad/sec.  The  coherence  function  shown 
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Fig.  5.2.  Roll-Rate  Response  to  Aileron  (p/-«a).  (a)  Magnitude; 

(b)  Phase. 

in  Fig.  5.3  is  very  strong  over  the  frequency  range  0. 1-8.5  rad/sec, 
indicating  excellent  identification  of  the  roll-rate  dynamics. 

The  perturbation  transfer-function  response  from  the  nonreal-time 
simulation  is  shown  with  the  flight  data  in  Fig.  5.2.  The  magnitude 
response  comparison  is  excellent  at  high  and  low  frequencies,  and  shows 
only  a  slight  error  at  mid-frequencies.  The  phase  comparison 
(Fig.  5.2b)  indicates  that  the  simulation  model  overpredicts  the  damping 
ratio  of  the  dominant  roll-response  mode  (u>  =  1  rad/sec).  The  phase 
error  at  high  frequency  is  due  to  the  omission  of  the  high-frequency 
actuator  dynamics  in  the  simulation  model.  These  various  small  discrep¬ 
ancies  in  Fig.  5.2  are  insignificant  because  it  is  not  possible  to 
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Fig.  5.3.  Coherence  Function  (y,  )  for  Roll-Rate  Response 
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perfectly  match  the  test  conditions  of  the  aircraft  with  the  simula¬ 
tion.  Clearly,  though,  the  simulation  response  characteristics  are  a 
very  good  representation  of  the  actual  vehicle  response. 

Aerodynamic  sideslip  measurements  are  obtained  from  a  sideslip 
indicator  located  about  18  ft  ahead  of  the  aircraft  center-of-gravity 
(c.g.).  The  sideslip  at  the  c.g.  ( 6-,, )  is  calculated  by  correcting  the 

eg 

measured  signal  for  position  error  based  on  yaw-rate  and  airspeed  (no 
flow-distortion  corrections  were  made).  Also,  no  corrections  are  made 
for  sensor  dynamics,  because  these  are  not  felt  to  be  significant  within 
the  identification  bandwidth. 

The  sideslip  response  (at  the  c.g.)  to  rudder  inputs,  8Cg/6p,  is 
shown  in  Fig.  5.4.  The  response  is  characterized  by  a  lightly  damped 
second-order  mode  with  a  frequency  of  about  1.6  rad/sec.  High-frequency 
rudder  inputs  yield  a  constant  sideslip  acceleration  with  no  appreciable 
time-delay,  and  low-frequency  rudder  inputs  yield  a  constant  sideslip 
angle.  The  coherence  function  (Fig.  5.5)  is  strong  over  the  frequency 
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Fig.  5.4.  Sideslip  Response  to  Rudder  (Bcg/-6r).  (a)  Magnitude 

(b)  Phase. 


Identification. 


range  0. 1-3.5  rad/sec,  falling  off  sharply  at  higher  frequencies.  This 
rapid  and  relatively  low-frequency  decline  of  the  coherence  function 
results  because  the  sideslip  response  falls  off  at  a  rate  of 
-40  dB/decade,  in  contrast  to  the  -20  dB/decade  roll-off  for  the  other 
response  variables.  However,  since  the  yaw  mode  of  interest  has  a 
natural  frequency  of  about  1.6  rad/sec,  the  falling  coherence  for  fre¬ 
quencies  greater  than  3.5  rad/sec  is  not  a  severe  limitation. 

The  simulation  frequency  response  of  8„_/6r  is  compared  with  the 

eg  r 

flight  data  in  Fig.  5.4.  As  in  the  roll  case,  the  comparison  is  gener¬ 
ally  excellent.  The  parallel  (vertical)  shift  in  the  magnitude  curves, 
suggests  that  the  rudder  yaw-response  sensitivity  is  overestimated  in 
the  mathematical  model  by  about  30%.  This  matches  the  rudder-response 
sensitivity  discrepancy  in  hover  and  so  further  indicates  a  yaw  inertia 
error.  The  slightly  more  shallow  phase  roll-off  in  the  simulation 
response  at  high  frequency  is  again  due  to  the  omission  of  actuator 
characteristics  and  perhaps  a  small  overestimation  in  the  damping  ratio 
of  the  Dutch  roll  mode.  Once  again,  however,  these  differences  are  not 
considered  to  be  significant  and  could  be  attributed  to  differences  in 
flight  and  loading  conditions  between  the  aircraft  and  the  simulation 
model . 

5.1.2  Transfer-Function  Models 

Equivalent  system  fitting  for  fixed-wing  aircraft  using  decoupled 
lateral/directional  models  has  been  considered  in  detail  by  Bischoff  and 
Palmer  (Ref.  15).  With  the  decoupled  model,  the  responses  are  fit 
independently,  using  a  first-order  roll  rate  and  a  second-order  sideslip 
transfer  function.  The  coupled  model  approach  is  based  on  simultaneous 
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fitting  of  the  roll  and  sideslip  responses  to  aileron  and  rudder, 
respectively,  i.e.,  we  obtain  the  following  fourth-order  transfer 
functions: 


6a(s)  -  (t/Ts)(1/Tr)[cd,  «d] 


(5.1) 


where  p(s)/6_(s)  is  the  Laplace-transformed  roll-rate  response  to 
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aileron  surface  deflection  (deg/sec/deg-aileron) ,  Lg  is  the  aileron 
roll  sensitivity,  and  are  the  equivalent  Dutch  roll-mode  damping 
and  natural  frequency,  respectively,  and  is  the  effective  time 

delay,  and 
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where  8cg(s)/6r(s)  is  the  Laplace-transformed  sideslip  response  to 

rudder  surface  deflection  (deg/deg-rudder) ;  is  the  rudder  sideslip 

r 

sensitivity;  1/Tfl  ,  1/T*  ,  and  1/T0  are  the  first-order  numerator 

D 1  02 

inverse  time  constants;  and  the  denominator  parameters  are  identical  to 
those  of  Eq.  (5.1).  The  effective  time  delay  for  the  sideslip  response 
is  Tg.  The  simultaneous  fitting  approach  is  consistent  with  that  used 
in  the  longitudinal  axis  and  allows  the  identification  of  the 
"coupling  effect"  in  the  roll  response  [Eq.  (5.1)],  which  is  important 
in  handling-qualities  assessments. 

The  coherence  function  results  of  Figs.  5.3  and  5.5  show  a  satis¬ 
factory  identification  of  both  roll  and  sideslip  responses  in  the 
frequency  range  0. 1-3.5  rad/sec.  Based  on  this  fitting  range,  the 
computer  program  LATFIT  is  used  to  obtain  the  parameters  of  Eqs.  (5.1) 


satisfactory  over  the  wider  range  of  0. 1-8.5  rad/sec,  this  degree  of 
freedom  is  refitted  alone,  holding  the  denominator  factors  constant  at 
the  values  obtained  from  the  simultaneous  solution.  This  procedure 
optimizes  the  values  of  the  high-frequency  gain  and  numerator  parameters 
in  the  roll  transfer  function.  The  final  results  for  the  parameters  of 
Eqs.  (5.1)  and  (5.2)  are 
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The  lateral/directional  transfer  function  fits  are  plotted  as 
dashed  lines  for  comparison  with  the  flight  data  of  Figs.  5.6  and  5.7. 
The  matching  of  the  roll-response  magnitude  (Fig.  5.6a)  is  excellent 
over  the  entire  frequency  range.  Only  a  slight  anomaly  in  the  phase 
matching  (Fig.  5.6b)  is  apparent  for  low-frequency  inputs.  The  agree¬ 
ment  between  the  sideslip  fit  and  the  flight  data  is  also  excellent 
(Fig.  5.7). 

As  a  result  of  the  simultaneous  fitting  of  the  roll  and  yaw 
responses,  the  transfer-function  parameters  of  Eqs.  (5.3)  and  (5.4)  have 
clearly  retained  their  physical  significance.  The  inverse  roll  mode 
time-constant  1/Tp  =  1.09  rad/sec  is  roughly  equal  to  the  frequency 
for  -45  deg  of  phase  lag  as  expected.  The  equivalent  Dutch-roll  mode  is 
lightly  damped,  with  a  natural  frequency  roughly  corresponding  to  the 
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Fig.  5.6.  Transfer-Function  Model  Identification  for  Roll-Rate 
Response  to  Aileron  (p/-«a).  (a)  Magnitude;  (b)  Phase. 

peak  in  the  sideslip  Bode  magnitude  plot  (Fig.  5.7a).  The  small  effec¬ 
tive  time  delays  for  the  roll  and  sideslip  responses  reflect  negligible 
high-frequency  dynamics,  and  supports  the  omission  of  corrections  for 
sideslip  sensor  dynamics. 

The  ratio  of  the  natural  frequency  of  the  numerator  complex  zero, 
to  that  of  the  denominator  complex  pole,  ud,  is  one  measure  of 
roll-yaw  coupling  in  response  to  aileron  inputs.  When  =  u>d,  the 
numerator  and  denominator  quadratic  factors  roughly  cancel,  and  the 
resulting  decoupled  roll  response  is  characterized  entirely  by  the  roll¬ 
mode  time-constant.  This  case  typically  leads  to  the  best  handling 
qualities  for  a  nominal  value  of  the  roll-mode  time-constant.  As  the 
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Fig.  5.7.  Transfer-Function  Model  Identification  for  Sideslip 
Response  to  Rudder  (6cg/-6r).  (a)  Magnitude;  (b)  Phase. 

roll-yaw  coupling  increases,  the  numerator  and  denominator  quadratic 
factors  of  Eq.  (5.1)  no  longer  cancel,  and  an  undesirable  oscillatory 
component  of  roll  rate  is  generated.  With  reference  to  Eq.  (5.3),  the 
near-unity  value  of  =  1.19  suggests  no  such  concern  for  roll-yaw 

coupling.  Therefore,  a  simple,  decoupled  first-order  roll-rate  model 
could  be  adopted  for  future  studies. 

The  poles  and  zeros  of  the  identified  roll-response  transfer  func¬ 
tion  are  shown  on  the  complex  plane  in  Fig.  5.8.  Roll-response  handling 
qualities  are  specified  in  terms  of  the  location  of  the  dominant  pole 
(1/Tr),  with  the  Level  I  (desirable  handling  qualities)  boundary 
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Fig.  5.8.  Location  of  Roll-Rate  Transfer-Function  (p/6_) 
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Parameters  in  the  s-Plane. 


(Ref.  17)  shown  in  the  figure.  These  results  indicate  very  good  lateral 
handling  qualities  for  the  cruise  flight,  in  contrast  to  the  poor  open- 
loop  lateral  handling  qualities  for  hover. 

The  identified  transfer  functions  of  Eqs.  (5.3)  and  (5.4)  are 
compared  with  the  simulation  perturbation  transfer  functions  in 
Table  5.1.  Comparison  is  seen  to  be  excellent  for  all  of  the  transfer 
function  parameters,  including  the  high-frequency  sideslip  response 
zeros.  This  comparison  is  a  reflection  of  the  very  close  frequency- 
response  results  presented  in  Figs.  5.2  and  5.4.  The  slightly  more 
shallow  phase  roll-off  for  both  the  roll  and  sideslip  responses 


'©y'-’ 1 


(Figs.  5.6b  and  5.7b)  is  reflected  in  the  Dutch  roll  damping  ratio 
discrepancy  (11?). 

5.2  Longitudinal  Dynamics 

5.2.1  Frequency-Response  Identification 

The  elevator  surface  to  pitch  rate  open-lcop  frequency  response 
q/«e  is  shown  in  Fig.  5.9-  The  short-period-mode  excitation  causes  a 
peak  in  the  pitch-rate  response  at  about  2.0  rad/sec,  with  an  associated 
phase  lag  of  -45  deg.  At  higher  frequencies,  the  magnitude  response 
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Fig.  5-9.  Pitch-Rate  Response  to  Elevator  (q/-6g).  (a)  Magnitude; 


(b)  Phase. 


rolls  off  at  -20  dB/decade  and  the  phase  shift  approaches  -90  deg  owing 
to  the  K/s  pitch-rate  characteristic  and  negligible  flexibility/ 
servolag  effects.  The  drop  in  the  magnitude  response  and  associated 
positive  phase  response  for  frequencies  below  0.3  rad/sec  are  due  to  the 
phugoid  dynamics. 

The  coherence  function  shown  in  Fig.  5.10  is  strong  over  the 
frequency  range  of  0. 2-7.0  rad/sec.  For  input  frequencies  above 


Identification. 

7.0  rad/sec,  the  coherence  function  becomes  erratic  and  the  transfer- 
function  identification  is  less  accurate.  For  low-frequency  inputs 
(less  than  0.3  rad/sec),  the  pitch-rate  response  decreases,  even  for  the 
nearly  constant  input  amplitude,  owing  to  the  effect  of  the  phugoid 
dynamics.  This  results  in  a  decrease  of  information  transfer  and  an 
associated  drop  in  the  coherence  function.  This  coherence  function 
roll-off  is  also  attributable  to  atmosphere  turbulence  effects,  which 
become  more  important  at  low  frequency  where  the  turbulence  spectrum 
peaks  (Ref.  46). 


The  pitch-rate  frequency  response  obtained  from  the  simulation 
perturbation  transfer  function  is  shown  with  the  flight  data  in 
Fig.  5.9.  The  comparison  is  excellent  in  the  frequency  range  of 
0. 3-1.0  rad/sec  showing  that  the  mathematical  model  correctly  predicts 
the  steady-state  pitch-rate  response  to  elevator  inputs.  The  differ¬ 
ences  at  very  low  frequency  are  not  conclusive  because  of  the  poor 
coherence  for  frequencies  below  0.2  rad/sec  (Fig.  5.10)  and  the  lack  of 
governor  dynamics  in  the  simulation  perturbation  transfer-function 
models.  In  the  frequency  range  of  1-10  rad/sec  there  is  a  significant 
shift  in  the  magnitude  and  phase  curves  between  the  aircraft  and  simula¬ 
tion  responses.  This  discrepancy  indicates  that  both  the  high-frequency 
gain  and  dominant  response  frequency  are  in  error  in  the  model.  How¬ 
ever,  the  fact  that  the  steady-state  gain  is  correct  suggests  a  dynamic 
response  rather  than  a  static  response  error.  This  observation  in 
addition  to  the  fact  that  very  good  static  trim  correlation  has  been 
achieved  in  cruise  (Ref.  5),  suggests  an  error  in  the  pitching  moment- 
of-inertia,  as  was  also  indicated  in  the  hover  results.  As  previously 
discussed,  ballast  which  is  used  in  the  aircraft  for  c.g.  control,  is 
not  included  in  the  present  simulator  data  base  and  is  a  likely  cause  of 
this  moment-of- inertia  error.  The  nearly  parallel  shift  between  the 
simulation  and  flight  data  curves  indicates  that  except  for  this  error, 
the  dominant  vehicle  dynamics  are  accurately  modeled. 

The  normal-acceleration  response  to  elevator,  az/6e,  is  shown  in 
Fig.  5.11.  The  response  is  dominated  by  the  classical  second-order 
short-period  mode  over  most  of  the  frequency  range.  The  magnitude  curve 
is  flat  at  mid-frequency,  indicating  a  constant  normal  acceleration 
(load  factor)  response  to  a  step  elevator  input,  with  a  roll-off  in 
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Fig.  5.11.  Vertical-Acceleration  Response  to  Elevator  (a„/60). 
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(a)  Magnitude;  (b)  Phase. 

response  for  frequencies  beyond  the  short-period  mode.  The  fall-off  in 
normal  acceleration  response  for  frequencies  below  0.3  rad/sec  is  prob¬ 
ably  a  result  of  the  dominance  of  the  phugoid  dynamics.  The  phase  curve 
exhibits  the  classical  second-order  characteristic,  that  is,  0  deg  of 
phase  lag  at  low  frequency,  180  deg  of  phase  lag  at  high  frequency,  and 
90  deg  of  phase  lag  at  the  second-order  mode  (u>  =  2.0  rad/sec).  This  is 
consistent  with  the  previous  pitch-rate  results.  The  coherence  function 
for  the  normal  acceleration  response  shown  in  Fig.  5.12  is  strong  over 
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Fig.  5.12.  Coherence  Function 
Response  Identification. 
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the  wide  frequency  range  of  0.3-10.0  rad/sec,  with  the  fall-off  at  low 
frequency  again  resulting  from  the  dominance  of  phugoid  dynamics  and 
turbulence  effects.  The  simulation  frequency-response  is  shown  with  the 
flight  data  in  Fig.  5.11.  As  in  the  pitch  case,  the  steady-state 
response  to  elevator  inputs  is  accurately  modeled  in  the  simulation 
(0. 3-1.0  rad/sec).  Also,  there  is  a  parallel  shift  at  high  frequency  in 
both  the  magnitude  and  phase  curves,  which  as  with  the  previous  pitch- 
rate  results,  indicates  an  error  in  the  assumed  pitching 
moment-of- inertia . 

5.2.2  Transfer-Function  Models 

Examination  of  the  longitudinal  frequency  responses  of  Figs.  5.9 
and  5.11  indicates  that  the  longitudinal  dynamics  for  this  flight  condi¬ 
tion  are  dominated  by  short-period  mode  (u>sp  =  2  rad/sec).  Therefore, 
the  classical  pitch-rate  and  center-of-rotation  normal-acceleration 
responses  to  elevator  are  applicable: 


(5.5) 
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where  the  q(s)/6e(s)  is  the  Laplace-transformed  pitch  rate  response  to 

elevator  surface  deflection  (deg/sec/deg-elevator) ;  M,  is  the  elevator 

°e 

pitch  sensitivity;  1/T0^  is  the  first-order  numerator  inverse  time- 
constant;  c,,n  and  ui  are  the  equivalent  short-period  mode  damping  and 
natural  frequency,  respectively;  and  t0  is  the  effective  time  delay. 
And 
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where  az(s)/<Se(s)  is  the  Laplace-transformed  vertical  (positive  down¬ 
ward)  acceleration  response  at  the  instantaneous  center-of-rotation  to 

elevator  surface  deflection  (g/deg -elevator ) ,  a„  is  the  elevator 

5e 

vertical  sensitivity  and  the  denominator  parameters  are  identical  to 

those  of  Eq.  (5.5).  The  effective  time  delay  for  the  vertical- 

acceleration  response  is  x_  .  The  low-frequency  phugoid  dynamics  are 

“z 

ignored,  since,  as  seen  in  Figs.  5.9  and  5.11,  these  are  important  only 
for  inputs  of  the  very  lowest  frequencies. 

The  transfer-function  parameters  of  Eqs.  (5.5)  and  (5.6)  are  itera¬ 
tively  varied  in  order  to  obtain  a  best  least-squares  fit  between  the 
equations  and  the  frequency-responses  of  Figs.  5.9  and  5.11  over  the 
selected  frequency  range;  commonality  of  the  denominator  parameters  is 
imposed.  For  a  first-cut  model,  the  center-of-gravity  (c.g.)  and 
longitudinal  instantaneous  center-of-rotation  (ICR)  are  assumed  to  be 
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coincident.  The  iterative  fitting  procedure  is  completed  using  the 
LONFIT  computer  program.  For  the  present  case,  the  selected  range  for 
simultaneous  fitting  of  the  pitch-rate  and  normal  acceleration  responses 
is  0.3-7  rad/sec.  In  this  frequency  range,  the  dynamics  are  clearly 
dominated  by  the  short-period  mode  and  the  coherence  is  strong  for  both 
transfer  functions.  Once  the  short-period  damping  and  natural  fre¬ 
quency  c  and  to  are  obtained  for  the  simultaneous  fit,  the  high- 

frequency  gain  for  the  normal  acceleration  response  a_  is  varied, 

% 

holding  the  damping  and  frequency  constant.  This  procedure  is  analogous 
to  that  used  in  the  lateral/directional  case  and  optimizes  the  az  fit 
over  the  frequency  range  0.3-10.0  rad/sec,  since  the  coherence  function 
for  this  measurement  remains  strong  out  to  higher  frequencies.  The 
transfer-function  parameters  for  the  pitch-rate  and  normal-acceleration 
responses  are  finally  obtained: 
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A  comparison  of  these  lower-order  models  with  the  flight-test  results 
is  presented  in  Figs.  5.13  and  5.14.  The  pitch-response  model 
[Eq.  (5.7)]  matches  the  data  very  well  over  the  selected  fitting  range 
(0. 3-7.0  rad/sec  for  q/6e),  which  shows  that  the  short-period  approxi¬ 
mation  adequately  represents  the  high-  and  mid-frequency  dynamics.  For 
low-frequency  inputs  (u>  <  0.3  rad/sec),  the  phugoid  mode  tdusea  a  drop 
in  the  pitch  response,  a  characteristic  which  is  not  "captured"  by  the 
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Fig.  5.13.  Transfer-Function  Model  Identification  for  Pitch-Rate 


Response  to  Elevator  (q/-6e).  (a)  Magnitude;  (b)  Phase. 
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Fig.  5.14.  Transfer-Function  Model  Identification  for  Vertical 
Acceleration  Response  to  Elevator  (az/6e) .  (a)  Magnitude;  (b)  Phase. 

short-period  model.  The  match  between  the  normal  acceleration  transfer 
function  [Eq.  (5.8)]  and  the  flight  data  is  not  nearly  as  good  in  the 
fitting  range,  and  shows  a  noticeable  discrepancy  in  the  magnitude 
response  for  frequencies  greater  than  2  rad/sec. 

Such  an  error  in  the  mid-  and  high-frequency  az  model  could  arise 
if  the  c.g.  and  ICR  are  not  coincident.  Then  the  vertical  acceleration 
response  at  the  c.g.  is  given  approximately  by 
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where  xCg  is  the  (unknown)  axial  location  of  the  center-of-gravity 
relative  to  the  instantaneous  center-of-rotation, 

az 

lCR 

— - — (s)  =  acceleration  response  at  the  ICR 
6e 

a 

=  -i(s)  of  Eq.  (5.6) 


§-(s)  =  Eq.  (5.5) 


The  frequency  responses  of  Figs.  5.13  and  5.14  were  refitted  with  the 
models  of  Eqs.  (5.5)  and  (5.9),  with  xno  included  in  the  parameter 

eg 

search.  The  resulting  model  of  pitch-rate  response  [Eq.  (5.5)]  is 
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The  vertical  acceleration  response  at  the  center-of-rotation  [Eq.  (5.6)] 
is 
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and  the  location  of  the  c.g.  relative  to  the  ICR  is  =  -5.73  ft 

eg 

(c.g.  behind  ICR).  The  resulting  acceleration  response  at  the  c.g.  is 
obtained  from  Eq.  (5.9): 
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which  is  a  higher  numerator-order  model  than  the  simple  model  of 
Eq.  (5.8),  and  displays  the  nonminimum  phase  (right-half  plane)  zero 
that  results  from  the  c.g.  being  located  behind  the  ICR. 

The  pitch-rate  model  of  Eq.  (5.10)  is  shown  in  Fig.  5.13  for  com¬ 
parison  with  the  previous  results.  Neither  the  frequency  response  nor 
the  transfer-function  parameters  of  the  q/5e  model  are  significantly 
altered  by  including  the  degree-of-freedom.  However,  the  c.g. 

acceleration  model  of  Eq.  (5.12)  fits  the  extracted  frequency  response 
much  more  closely,  as  shown  in  Fig.  5.14.  Thus  the  noncoincident  c.g. 
and  ICR  appear  to  be  a  likely  source  of  the  discrepancy  obtained  when 
the  simple  model  of  Eq.  (5.8)  is  used. 

The  identified  short-period  damping  and  frequency  (t  =  0.536  and 

sp 

ui  =  2.02  rad/sec)  are  consistent  with  preliminary  observations,  which 
sp 

were  based  on  the  raw  Bode  plot  information.  The  relatively  high  degree 
of  damping  is  a  reflection  of  the  small  response  peak  of  the  pitch-rate 
transfer  function.  The  small  ei  'ective  time  delays  (in  the  pitch  atti- 
tude^nd  normal  acceleration  transfer-function  fits)  indicate  that  the 
high-frequency  flexibility  effects  are  not  important  for  this  flight 
condition,  as  noted  earlier.  The  military  standards  requirements 
(Ref.  17)  for  lower-order  equivalent  system  pitch  response  are  given  in 
terms  of  the  parameters  “gpTg^*  tSp,  tQ.  The  values  of  these  parame¬ 
ters  given  in  Eq.  (5.10)  are  well  within  the  Level  I,  category  A,  indi¬ 
cating  desirable  longitudinal  handling  for  cruise  flight. 

The  identified  transfer  functions  are  compared  with  the  simulation 
perturbation  transfer  functions  in  Table  5.1.  The  simulation  transfer 
functions  contain  an  additional  complex  mode  associated  with  the  phugoid 
dynamics  and  a  nearly  canceling  pole/zero  pair  associated  with  the  rotor 


angular  velocity  (rpm)  degree-of- freedom.  The  fact  that  the  rpm  mode 
nearly  cancels  and  that  the  phugoid  mode  is  at  very  low  frequency  justi¬ 
fies  ignoring  these  dynamics  and  adopting  the  simple  short-period 
models . 

The  pole  and  zero  locations  of  the  simulation  transfer  functions 
compare  very  well  with  the  flight-extracted  results.  The  good  corre¬ 
spondence  of  the  numerator  factors  in  the  vertical  acceleration  response 
supports  the  corrections  based  on  instantaneous  center-of-rotation 
effects  [Eqs.  (5.9),  (5.12)].  The  differences  in  short-period  frequency 
and  pitch-response  sensitivity  can  both  be  explained  by  the  error  in 
pitching  moment-of- inertia.  The  pitch-axis  sensitivity  Mg  is  the 
pitching  moment  (foot-pounds)  per  inch  of  elevator  normalized  by  the 
pitching  inertia  (Iy)  (Ref.  41).  The  ratio  of  pitch  sensitivities  for 
the  simulation  and  flight  results  (Table  5.1)  gives: 
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The  short-period  frequency  w  depends  on  the  square-root  of 
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stability  derivatives  which  are  also  normalized  by  the  pitching  moment- 
of-inertia  (Ref.  41).  Based  on  the  short-period  frequency: 
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which  is  the  same  result  as  in  Eq.  (5.13).  Thus,  an  error  in  pitch 
inertia  would  explain  the  differences  in  both  pitch  sensitivity  and 
short-period  frequency,  while  not  disturbing  the  good  static  trim  cor¬ 
relation.  An  underestimation  of  pitch  inertia  in  the  simulation  is 


indicated  from  both  the  hover  and  cruise  results,  although  there  is  some 
discrepancy  in  the  error  ratio  (67$  from  the  cruise  results  and  25$  from 
the  hover  results). 


5.3  Verification  of  Transfer-Function  Models 

The  high  degree  of  inherent  bare-airframe  stability  for  the  cruise 
flight  condition  makes  open-loop  verification  inputs  much  easier  to 
complete  than  in  the  hover  condition.  Steady  initial  conditions  and 
clean  long-duration  (6-8  sec)  step  inputs  are  easily  obtained.  The 
aircraft  (filtered)  roll-rate  response  to  a  step  aileron  input  is  shown 
in  Fig.  5.15.  The  transfer-function  model  of  Eq.  (5.3)  shown  in  the 
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Fig.  5.15.  Roll-Response  Model  Verification  Using  SCAS-off  Flight 


Data,  (a)  Aileron  Input;  (b)  Roll-Rate. 


dashed  line  matches  the  flight  data  quite  well  for  the  entire  input 
duration,  except  for  some  discrepancy  in  steady-state  roll  rate.  This 


error  may  be  associated  with  the  p/6_  transfer-function  model  mismatch 
at  mid-frequency  (Fig.  5.6),  some  nonlinear  effects  associated  with  the 
large  steady-state  roll  rate,  or  simply  a  gust.  The  sideslip  response 
to  a  step  rudder  input  is  shown  in  Fig.  5.16,  with  the  appropriate 
corrections  for  the  nose-boom  position.  The  transfer-function  model  of 
Eq.  (5.4)  matches  the  flight  data  almost  identically  for  the  entire 
time-history. 
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Fig.  5.16.  Sideslip-Response  Verification  Using  SCAS-off  Flight 
Data,  (a)  Rudder  Input;  (b)  Sideslip  Response. 


For  most  of  the  flight  tests,  frequency-sweeps  and  step  inputs  were 
performed  consecutively  for  each  flight  condition,  thus  maintaining 
constant  altitude,  airspeed,  and  loading  condition.  Unfortunately,  the 
open-loop  elevator  step  response  data  for  =  170  knots  were  taken 
at  Hd  =  2800  ft  (density  altitude),  whereas  Hd  =  10,000  ft  for  the 
sweeps.  Since  the  indicated  airspeed  was  constant  and,  therefore,  also 
the  dynamic  pressure,  the  derived  pitch-rate  model  is  valid  for  the 
lower-altitude  step  responses.  However,  the  significant  difference  in 
true  airspeed  (about  20  knots)  and  associated  centripetal  acceleration 
(=  UQq)  makes  an  az  step  response  comparison  invalid  for  these  data. 
The  response  to  elevator  doublet  inputs,  obtained  from  closed-loop 
(SCAS-on)  tests,  for  =  170  knots  and  Hd  =  12,000  ft  is  used  to 
validate  the  az  model,  although  the  aircraft  loading  condition  (and 
altitude)  is  still  not  exactly  matched. 

The  SCAS-off  (open-loop)  pitch-rate  response  to  step  elevator 
deflection  (Hd  =  2800  ft)  is  shown  in  Fig.  5.17.  The  elevator  input 
corresponds  to  about  40%  of  the  maximum  control  deflection.  The  compar¬ 
ison  between  the  ICR-corrected  model  response  [Eq.  (5.10)]  and  the 
flight  data  is  seen  to  be  very  good  over  most  of  the  time-history.  The 
slight  deviations  occurring  toward  the  end  of  the  run  are  due  to  the 
inadequacy  of  the  short-period  approximation  in  modeling  the  low- 
frequency  (phugoid)  dynamics  (see  Fig.  5.13a).  Even  so,  the  short- 
period  approximation  clearly  gives  an  excellent  characterization  of  the 
important  initial  response  dynamics. 
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Fig.  5.17.  Pitch-Response  Model  Verification  Using  SCAS-off  Flight 
Data  (hd  =  2800  ft),  (a)  Elevator  Input;  (b)  Pitch  Rate. 

The  pitch-rate  model  also  compares  favorably  with  the  flight  data 

for  the  response  to  an  elevator  doublet  input,  as  shown  in  Fig.  5.18c. 

However,  the  normal-acceleration  model  of  Eq.  (5.12)  is  not  very  good. 

The  initial  reversal  predicted  by  the  model  (effect  of  the  xno  correc- 

tion)  does  not  occur  in  the  flight  response  in  Fig.  5.18c.  In  fact,  if 

the  original  center-of-rotation  model  of  Eq.  (5.8)  is  used,  the  match 

with  the  flight  data  is  much  better,  as  is  also  shown  in  Fig.  5.18c. 

These  results  indicate  that  the  a  /6_  model  which  includes  the 

eg  e 

center-of-rotation  correction  [Eq.  (5.12)]  is  highly  tuned  to  the 
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Fig.  5.18.  Pitch  and  Vertical-Acceleration  Response  Model 
Verification  Using  SCAS-on  Flight  Data  (hd  =  12000  ft),  (a)  Elevator 
Input;  (b)  Pitch  Rate;  (c)  Vertical-Acceleration  Response. 


measured  frequency  response  for  sinusoidal  inputs  (Fig.  5. I1*)  and  is  not 
robust  to  small  changes  in  flight  condition  and  input  form.  The  fact 
that  acceleration  is  a  higher-order  state  derivative  than  the  other 
angular-rate  derivatives  further  exaggerates  the  sensitivity  of  this 
model,  especially  for  the  initial  (high-frequency)  response  dynamics. 
Therefore,  the  lower-order  model  of  Eq.  (5.8)  should  be  used  for  future 
handling  qualities  and  control-system  analyses  of  this  flight  condi¬ 
tion.  This  is  a  good  example  of  the  dangers  of  model  overparameteriza- 
tion.  The  capability  to  evaluate  directly  the  effect  of  model  order  on 
response  fitting  is  also  clearly  seen  in  this  example,  again  represent¬ 
ing  a  unique  advantage  of  the  frequency-domain  approach  for  this  kind  of 
comparison. 

Overall,  these  verification  results  indicate  that  the  derived 
lower-order  transfer  functions  are  very  good  predictors  of  the  open-loop 
vehicle  response  in  cruise  over  a  wide  frequency  range. 


Chapter  6 
CONCLUSIONS 


The  following  conclusions  are  separated  into  those  specifically 
associated  with  the  analysis  of  the  XV- 15  and  those  associated  with  the 
analysis  of  the  frequency-response  identification  method  in  general. 

6.1  Analysis  of  the  XV-15 

Specific  findings  from  the  identification  of  the  open-loop  hover 
dynamics  are  as  follows. 

1 .  Heave  response  to  power  lever  and  yaw  response  to  rudder  are 
decoupled  from  the  other  degrees-of-freedom  and  are  very  lightly  damped. 

2.  Pitch  and  roll  open-loop  dynamics  exhibit  unstable  low- 
frequency  oscillations.  High  pilot  work-load  should  be  expected  in  the 
event  of  a  SCAS  failure  in  either  of  these  axes. 

3.  Significant  coupling  of  yaw  response  to  aileron  has  been  iden¬ 
tified.  Pedal  deflection  equal  to  about  50%  of  lateral  stick  inputs  is 
necessary  to  maintain  a  fixed  heading. 

4.  The  effective  high-frequency  time-delays  are  small  in  all  axes 
and  are  probably  not  important  in  determining  the  handling-qualities  of 
the  XV-15  aircraft. 

5.  The  comparison  of  piloted  simulation  and  flight-test  dynamic 
results  for  high-frequency  inputs  (u>  >  1.0  rad/sec)  is  generally  favor¬ 
able  in  all  axes,  with  larger  anomalies  existing  for  low-frequency 
inputs.  However,  underestimates  in  the  pitch  and  yaw  inertias  used  in 
the  simulation  model  caused  the  overprediction  of  pitch-  and  yaw-control 


sensitivities. 


6.  Small-perturbation  transfer  functions  obtained  from  the 
piloted-simulator  (VMS)  are  poor  representations  of  the  low-frequency 
dynamics.  However,  perturbation  transfer  functions  for  the  current 
nonreal-time  simulation  compare  more  favorably  with  the  flight  data, 
suggesting  that  the  corrections  to  the  mathematical  models  have  been 
appropriate . 

7.  The  step  responses  of  the  identified  transfer-functions  and  the 
aircraft  generally  match  quite  well,  even  for  large-amplitude  motions. 
This  verifies  that  an  accurate  representation  of  the  open-loop  hover 
dynamics  has  been  achieved,  and  that  the  hover  transfer-function  models 
derived  herein  can  be  used  for  future  XV- 15  control  system  and  handling- 
qualities  analyses. 

Specific  findings  from  the  identification  of  the  open-loop  cruise 
(V^  =  170  knots)  dynamics  are  as  follows. 

1.  The  response  characteristics  are  very  stable  and  decoupled,  and 
indicate  desirable  open-loop  handling-qualities. 

2.  Standard  lower-order  equivalent  models  adequately  match  the 
identified  dynamics  in  all  degrees-of-freedom  except  for  the  normal 
acceleration  responses  to  elevator  deflection.  Although  a  center-of- 
rotation  correction  is  needed  to  accurately  fit  the  az/6e  freQuency_ 
response  data,  the  resulting  higher-order  model  is  not  robust  to  small 
changes  in  flight  condition  and  input  form. 

3.  Small  perturbation  transfer  functions  and  associated  frequency 
responses  from  the  nonreal-time  simulation  correlate  very  well  with  the 
flight-extracted  results  except  for  an  error  in  the  simulation  pitch  and 
yaw  inertias. 


4.  The  step  responses  of  the  identified  transfer  functions  and  the 
aircraft  generally  match  very  well. 

6.2  Frequency-Response  Identification  Method 

1.  Non-parametr ic  frequency-response  identification  techniques 
produce  an  excellent  physical  understanding  of  basic  vehicle  dynamics. 

2.  Non-parametric  techniques  are  especially  useful  for  documenting 
the  response  characteristics  of  a  new  configuration  for  which  a  prior 
parametric  model  is  not  known. 

3.  Extraction  of  open-loop  responses  from  closed-loop  data  is 
feasible  with  frequency -sweep  inputs.  Closed-loop  testing  of  rotorcraft 
with  highly  unstable  bare-airframe  dynamics  allows  the  execution  of  long 
duration  inputs  needed  for  low-frequency  response  identification. 

4.  Multi-input/multi-output  frequency-domain  methods  are  suitable 
for  extracting  the  coupled  responses  of  rotorcraft. 


Chapter  7 


RECOMMENDATIONS  FOR  FUTURE  STUDY 

The  objectives  of  the  XV-15  study  (Sec.  1.1)  have  been  achieved. 

The  following  recommendations  pertain  to  future  studies  of  the 
frequency-response  identification  method. 

Future  studies  in  frequency-response  identification  should  focus  on 
advancing  the  method  in  two  major  areas: 

1.  Reconstruction  of  physical  stability-derivatives  from  the 
transfer-function  matrix 

2.  Generalization  of  two-input/single-output  frequency-response 
identification  to  multi-input/multi-output  (MIMO)  systems 

7.1  Reconstruction  of  Physical  Stability  Derivative  Parameters 

The  frequency-response  method  illustrated  in  this  study  yields  a 
transfer-function  matrix  for  the  coupled  system.  The  next  logical  step 
in  a  full  identification  procedure  is  to  reconstruct  the  physical 
stability-derivative  (state-space)  realization  rather  than  a  canonical 
realization  as  in  the  present  aproach.  The  MIMO  state-space  reconstruc¬ 
tion  approach  given  by  Chen  (Ref.  47)  is  well  suited  for  this  purpose 
because  it  takes  into  account  the  physical  redundancy  in  the  transfer- 
function  zeros.  Frequency-domain  parameter  identification  would  thus  be 
completed  in  a  three  step  procedure: 

Step  1.  Frequency-response  identification 

Step  2.  Transfer-function  modeling 

Step  3-  Stability-derivative  reconstruction 


By  breaking  the  parameter-identification  approach  into  these  three 
distinct  steps,  maximum  benefit  can  be  gained  from  the  nonparametric 
identification  before  the  parametric  step  is  attempted.  This  approach 
may  yield  both  a  more  accurate,  and  physically  meaningful  and  consistent 
state-space  model  than  can  be  identified  using  direct  time-domain  or 
frequency-domain  parameter-identification  procedures. 

Stability-derivative  reconstruction  should  be  tested  on  the  XV-15 
transfer-functions  for  the  hover  flight  condition,  since  a  complete 
state-space  model  has  previously  been  identified  for  this  condition 
using  time-domain  identification  methods  (Ref.  22).  The  frequency- 
domain  parameter-identification  procedure  should  then  be  tested  on  a 
single-rotor  helicopter  which  has  more  highly-coupled  system  dynamics 
than  the  XV-15. 

7.2  Generalization  to  MIMO  Systems 

Spectral  quality  assessment  is  conducted  on  each  input/output 
record  pair  to  determine  the  runs  suitable  for  concatenation  and  best 
window  sizes.  Although  this  assessment  is  relatively  straightforward, 
the  process  is  computer  and  labor  intensive,  and  could  be  automated 
using  a  rule-based  algorithm.  This  would  be  especially  useful  for  MIMO 
identification  where  the  number  of  required  frequency  responses 
increases  rapidly  with  additional  cross-coupling  paths. 

The  two-input  identification  software  tools  developed  (or  adapted) 
in  this  study  (TISOSA,  TFTHISTORY,  NAVFIT,  LONFIT,  LATFIT)  need  to  be 
further  improved  to  allow  MIMO  identification  and  verification  of  sys¬ 
tems  with  more  than  two  control  inputs.  A  general  matrix  formulation 
will  allow  extension  to  the  multi-input  case. 


Multi-input/multi-output  frequency-response  identification  tools 
should  be  tested  on  a  modern  single-rotor  helicopter  to  assess  their 
suitability  for  highly-coupled,  hingeless  (or  bearingless)  rotor 
systems. 
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Appendix  A 
CHIRP  z-TRANSFORM 


Existing  computer  programs  for  identifying  frequency-responses  from 
flight  data  generally  use  the  standard  (Cooley-Tukey )  FFT  procedure  to 
transform  data  from  the  time-domain  to  the  frequency-domain  (Refs.  10, 

11,  13).  Initial  versions  of  the  FRESPID  program  also  used  a  standard 
FFT  algorithm,  but  this  was  later  replaced  by  the  chirp  z-transform 
(CZT).  The  mathematical  basis  for  this  algorithm  is  given  in  Refs.  28 
and  32,  and  is  not  discussed  in  this  appendix.  The  CZT  provides  more 
flexibility  in  selection  of  computational  parameters  than  is  allowed  by 
the  standard  FFT.  However,  it  has  increased  computational  time  and 
storage  requirements. 

Three  important  differences  between  the  FFT  and  the  CZT  are 
reviewed  in  this  appendix:  (1)  Resolution,  Af;  (2)  Starting  frequency, 
fQ;  and  (3)  Window  size,  LAt.  These  parameters  define  the  discrete 
Fourier  transform  calculation  (DFT),  which  is  repeated  from  Eq.  (3-2): 

N-1 

X(f  )  =  X(kAf)  =  At  5^  x  exp[-j2*(kn)/N]  ;  k  =  0,  1,  2,...,  N-1 
k  n=0  n 

(A.1) 

where 

X(f|{)  =  Fourier  coefficients 
xn  =  x(nAt) 

At  =  time  increment 

N  =  number  of  discrete  frequency  points 


A.1  Resolution 


The  (frequency)  resolution  of  the  standard  FFT  (Af)  is  determined 
by  the  width  of  the  window  (LAt,  sec),  since  for  this  algorithm  the 
number  of  time-history  points  (L)  is  also  the  same  as  a  number  of  output 
spectral  points  (N): 

Af  =  LAT  ;  Hz  (A,2) 

which  in  terms  of  the  sampling  frequency  f_  =  1/At  is 

s 


Thus,  the  FFT  output  frequencies  are  N  evenly  spaced  points 

around  the  unit  circle.  Since  the  z-transform  is  a  mirror  image  with 

respect  to  the  Nyquist  frequency  (f  =  f  /2),  the  second  half  of  the 

c  s 

spectrum  (f  /2  <  kAf  <  f  )  is  redundant  and  is  not  actually  calculated 
s  s 

in  the  FFT  procedure  of  Eq.  (A.1).  Furthermore,  the  Nyquist  frequency 

is  generally  selected  to  be  much  greater  than  the  upper  end  of  the 

identification  frequency  range  (0.1-10.0  rad/sec  =  0.016-1.59  Hz  in  the 

present  study)  to  avoid  digital  distortion,  so  many  of  the  spectral 

calculation  points  are  not  needed  (1.59  <  kAf  <  f  /2). 

s 

The  basic  difference  in  the  CZT  algorithm  is  that  the  spectral 
calculation  is  completed  along  an  arbitrary  specified  arc  in  the 
z-plane.  When,  as  in  the  present  case,  this  arc  is  taken  as  a  segment 
of  the  unit  circle,  the  CZT  provides  the  Fourier  coefficients.  A  key 
characteristic  of  the  CZT  is  that  the  number  of  calculation  points  N 


within  this  segment  is  arbitrary;  thus,  the  frequency  resolution  is  also 


arbitrary t : 


Af 


(A.14) 


where 

fg  =  starting  spectral  frequency  (arbitrary,  Sec.  A. 2) 
fl  =  ending  spectral  frequency  (arbitrary,  Sec.  A. 2) 

A  commonly  held  misconception  is  that  arbitrary  frequency  resolu¬ 
tion  can  also  be  achieved  using  the  standard  FFT  by  simply  "padding"  the 
input  sequence  with  zeroes  (artificially  increasing  L).  However,  pad¬ 
ding  with  zeroes  only  allows  interpolation  between  true  frequency  reso¬ 
lution  points  ( Af  =  fs/b)  and  does  not  increase  the  frequency  resolu¬ 
tion.  The  only  way  to  increase  the  true  frequency  resolution  of  the 
standard  FFT  is  to  increase  the  width  of  the  window,  which  decreases  the 
amount  of  averaging  and  so  increases  the  random  error  (Sec.  3.2.3). 


A. 2  Starting  Frequency 

The  minimum  output  frequency  fg  that  can  be  identified  using  the 
standard  FFT  procedure  is  also  the  frequency  spacing: 


(A. 5) 


t  The  author  has  encountered  numerical  problems  with  the  CZT  algorithm 
when  the  number  of  output  points  (N)  is  much  greater  than  the  number  of  input 
points  (L),  N  j>  2L.  However,  since  all  the  output  points  are  concentrated 
within  the  desired  frequency  range,  the  resulting  frequency  resolution  is 
still  substantially  higher  than  is  achievable  with  the  standard  FFT  procedure. 


Reducing  the  starting  frequency  to  allow  an  identification  of  the 
lower  frequency  dynamics  means  selecting  wider  windows  and,  therefore, 
reducing  the  amount  of  spectral  averaging  [and  increasing  the  spectral 
variance,  Eq.  (3.10)].  In  the  CZT,  the  initial  (and  final)  frequency 
fg  (and  f 1 )  is  arbitrary.  Thus,  the  initial  frequency  is  independently 
selected  to  correspond  to  the  lower  end  of  the  frequency  range  of  inter¬ 
est  (fg  =  0.1/2n  =  0.016  Hz  in  this  case). 

A. 3  Window  Size 

In  both  the  standard  FFT  and  CZT  procedures,  window  size  (LAt) 
controls  the  balance  between  low-frequency  signal  content  and  spectral 
averaging  (Sec.  3.2.1).  In  the  standard  FFT,  the  number  of  time-history 
points  in  the  window  must  be  an  integer  power  of  2  (a  "composite 
number,"  e.g.,  512,  1024,  2048).  Thus,  the  allowable  adjustment  in  the 
width  of  window  is  coarse  (only  halving  or  doubling  in  size),  for  a 
fixed  sample  interval  (At).  The  chirp  z-transform  allows  the  arbitrary 
selection  of  L  (including  prime  numbers).  The  capability  to  indepen¬ 
dently  set  the  starting  frequency  (fg,  Sec.  A. 2)  and  then  finely  tune 
the  window  size  yields  both  improved  coherence  and  lower  spectral  vari¬ 
ance  over  a  wide  frequency  range  as  compared  to  the  standard  FFT. 

These  three  characteristics  of  the  chirp  z-transform  make  it  a 
significantly  more  flexible  algorithm  than  the  standard  FFT;  this  flexi¬ 
bility  can  yield  a  substantial  improvement  in  the  spectral  resolution 
and  dynamic  range,  especially  at  the  low-frequency  end.  The  author  has 
found  the  flexibility  of  the  chirp  z-transform  to  be  very  important  in 
achieving  high-quality  frequency-response  identification  from  (noisy) 
flight  data. 


Appendix  B 

LIMITATIONS  IN  CLOSED-LOOP  IDENTIFICATION 


The  identification  of  low-frequency  helicopter  dynamics  requires 
long-duration  inputs.  In  hover  flight  the  dominant  open-loop  dynamics 
are  unstable,  so  execution  of  these  inputs  with  the  stability  and  con¬ 
trol  augmentation  system  disengaged  (SCAS-off)  is  not  practical.  There¬ 
fore,  it  is  desirable  to  extract  the  open-loop  dynamics  from  flight-test 
data  obtained  with  the  augmentation  system  engaged  ("closed-loop  identi¬ 
fication"  ) . 

This  appendix  derives  the  fundamental  relationships  for  single¬ 
input/single-output  closed-loop  identification  in  the  presence  of  pro¬ 
cess  noise.  A  simple  relationship  relating  bias  errors  in  the  identi¬ 
fied  open-loop  frequency-response  to  the  input-to-process  noise  ratio  is 
developed.  A  numerical  study  illustrates  the  important  concepts  using  a 
simple  model  of  the  XV- 15  unstable  dynamics  in  hover. 

B.1  Spectral  Relationships 

Consider  the  identification  of  the  single-input/single-output  roll- 
rate  response  to  aileron  from  closed-loop  tests  as  illustrated  in 
Fig.  B.1.  The  closed-loop  system  is  the  same  as  that  previously  dis¬ 
cussed  in  Sec.  3.2.4  (Fig.  3.7),  except  in  the  present  case  the  compen¬ 
sation  element  (Gc)  is  placed  in  the  forward  loop.  This  provides  a 
steady-state  closed-loop  gain  of  approximately  unity-- independent  of  the 
compensation  gain.  The  other  elements  of  Fig.  B.1  are  as  follows. 


1.  =  pilot  lateral-stick  input 

2.  p  =  measured  roll-rate  (measurement  noise  is  neglected  here) 


Fig.  B.l.  Closed-Loop  Identification. 


3.  =  total  aileron  surface  deflection 

cl 

4.  H  =  pfi/6a  =  °Pen-^0°P  I'oll-rate  frequency-response 

5.  n  =  roll-rate  process  noise. 

Some  common  sources  of  noise  are:  (a)  turbulence,  the  auto¬ 
spectrum  of  n  is  determined  from  the  turbulence  velocity  autospectrum 
and  the  roll-rate  response  to  turbulence  velocity  inputs;  (b)  neglected 
secondary  inputs;  (c)  nonlinear  remnants  which  are  not  accounted  for  by 
the  linear  open-loop  frequency -response  H.  The  external  inputs  from 
the  pilot  and  process  noise  sources  are  assumed  to  be  uncorrelated; 

viz . ,  G.  =0. 

6Ln 

The  desired  open-loop  frequency-response  estimate  is  obtained  from 
Eq.  (3-8): 


H  = 


«  P 


6  6 
a  a 


(B.l) 


The  cross-spectrum  G&  ^  is  determined  from  Fig.  B.l  (Ref.  28): 


G  =  G  +  G 
Sap  5ap6  6an 


(B.2) 
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Substituting  Eq.  (B.2)  into  Eq.  (B.1): 


H 


6ap6 


6  n 
a 


6  6 
a  a 


H  = 


L 


6  n 
a 

'6  6 
a  a 


True  response  Bias  error 


(B.3) 


(B.4) 


which  shows  that  the  spectral  estimate  is  biased  by  correlation  between 
the  aileron  surface  and  process  noise  inputs.  This  bias  term  is  zero 
for  open-loop  systems  since  the  pilot  and  process  noise  inputs  are 
uncorrelated;  the  undesirable  correlation  results  from  the  feedback  of 
the  process  noise  around  the  loop.  Equation  (B.4)  also  shows  that  the 
true  (unbiased)  open-loop  frequency-response  is  identified  regardless  of 
the  level  of  compensation  when  process  noise  is  not  present. 

The  cross-spectrum  G6  n  is  determined  from  the  closed-loop  rela- 

cl 

tionship  between  process  noise  inputs  and  aileron  surface  response 
(Ref.  28) 


(B.5) 


t  Asterisk  (*)  denotes  complex-conjugate. 
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The  aileron  surface  input-autospectrum  is  composed  of  components 


from  the  two  external  (uncorrelated)  excitations: 

.2 


6  6 
a  a 


=  , 

6 

a 
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+  G 
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6l6l 

6L 

nn 
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nn 
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UGH 
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c 

c 

(B.6) 


Substituting  Eqs.  (B.5)  and  (B.6)  into  Eq.  (B.4) 


H  =  H  - 


G  [G  /(I  -  G  H)J 
nn  c  c 


|G  /(I  +  G  H)r(Gx  ,  ♦  G  ) 

1  c  c  1  nn 


(B.7) 


But, 


GC  | 2  i 

’  G 

c 

WGo  x 

1  -  GCH|  *  \ 

UGH 
'  c 

)  V1  +  GCH/ 

(B.8) 


so, 


H  =  H  - 


nn 


[v<’  ♦cch)'[vl*g""1 


(B.9) 


When  the  closed-loop  system  is  excited  only  by  process  noise  (i.e., 

no  pilot  inputs),  G.  =0  and 

5L6L 


H  =  H  - 


nn 


tcc/(1  *G0H)]Gnn 


(B.10) 


"  1 
H  =  -s; 
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which  is  the  limiting  condition  of  identifying  the  inverse  compensa¬ 
tion.  This  result  is  in  agreement  with  the  simple  analysis  of 
Eqs.  (3.11  and  3.12). 


When  pilot  and  process  noise  inputs  are  both  present,  we  can  define 
the  signal-to-noise  ratio: 

Gs  & 

2  L  L 

S*(o»)  E  r~v  Y  (B.  11) 

G  (w) 
nn 

which  depends  on  frequency. 

Substituting  Eq.  (B.11)  into  Eq.  ( B . 9 )  and  normalizing  the  estimate 
by  the  true  frequency-response  yields  the  following  normalized  bias 
error : 


e 


b 


bill 

H 


_ -_1 _ 

[G  H/(  1  +  G  H)][S2(u>)  ♦  1] 
c  c 


(B. 12) 


where  it  should  be  recognized  that  the  denominator  term  in  Eq.  (B.12) 
GCH/[1  +  GCH]  is  the  closed-loop  transfer-function.  The  open-loop 
unstable  modes  of  rotorcraft  are  generally  at  frequencies  well  below  the 
open-loop  crossover  frequency  wc  (defined  as  frequency  at  which 
|GCH|  =  1)  for  hovering  vehicles.  In  the  case  of  the  XV-15,  the  domi¬ 
nant  open-loop  unstable  modes  are  at  a  frequency  of  about 
u>  =  0.5  rad/sec  [Eq.  (4.4)],  compared  to  a  crossover  frequency  of 
u>c  =  2.5  rad/sec  (closed-loop  time-constant  =  0.4  sec,  Sec.  2.1.4).  At 
low  frequency  (w  <  (dc),  the  closed-loop  transfer-function  becomes 


1 


(B. 13) 


G  H 
c 


1  +  G  H 
c 


which  when  substituted  into  Eq.  (B.12)  gives  the  following  simple 
expression  for  the  normalized  bias  error  in  the  frequency  range  of  the 
open-loop  dynamics: 


1 


b  S2(w)  +  1 


(B.  1*1) 


Equation  (B.14)  indicates  that  the  normalized  bias  error  can  be 
kept  to  within  1 0%  (eb  <  0.10)  for  a  signal-to-noise  ratio: 


SU)  >  3 


(B. 15) 


In  terms  of  the  noise-to-signal  ratio  o(w)  =  1/S(u>),  the  requirement  is 


5(u)  <  0.33 


(B. 16) 


which  was  quoted  in  Sec.  3.2.4. 

This  analysis  shows  that  the  pilot's  input  must  exceed  the  input 
from  process  noise  by  at  least  a  factor  of  3  across  the  frequency  range 
of  interest  to  reduce  the  normalized  bias  error  to  an  acceptable 
value.  The  use  of  the  frequency-sweep,  which  generates  a  persistent  and 
relatively  large  vehicle  response  over  the  entire  identification  fre¬ 
quency  range,  ensures  that  this  condition  is  maintained  and  that  the 
bias  errors  are  small.  Multi-step  inputs  generally  do  not  have  uniform 
power  distribution.  If  process  noise  inputs  are  large  in  the  frequency 
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range  where  the  multi-step  signal  content  is  small,  large  bias  errors 
may  result,  making  such  inputs  poorly  suited  for  closed-loop 
identification. 

B.2  Numerical  Study 

A  numerical  study  was  conducted  to  illustrate  the  concepts  dis¬ 
cussed  in  the  preceding  section  and  to  validate  the  expression  for 
normalized  bias  error  given  in  Eq.  (B.14).  A  simple  model  of  the  XV-15 
dynamics  in  hover  is  used  to  facilitate  a  clear  physical  understanding 
of  the  limitations  of  closed-loop  identification,  which  is  important  for 
future  applications. 

B.2.1  Simplified  Closed-Loop  Model 

The  open-loop  roll-rate  response  in  hover  [ Eq .(4.4) ]  is  represented 
by  a  simple  (normalized)  first-order  unstable  transfer-function: 

-0.005s 

H(s)=rr-oX  <B'17) 

where  the  5  msec  time-delay  is  included  to  account  for  numerical 
integration. 

The  compensation  transfer-function  is  assumed  to  be  a  pure  gain: 

Gc  =  K  (B. 18) 

The  closed-loop  time-constant  of  the  XV-15  roll-response  in  hover 


is  0.4  sec  (Sec.  2.1.4),  which  implies  a  closed-loop  first-order  roll¬ 
mode  at  s  =  -2.5  rad/sec,  yielding  a  compensation  gain  of: 


K  =  3.0 


(B. 19) 


The  pilot  lateral-stick  excitation  was  derived  by  passing  a 
Gaussian  random  sequence  (200-sec  duration)  through  a  low-pass  filter 
(1.0  rad/sec  roll-off  frequency).  This  yields  an  input-autospectrum 
(G  )  approximating  that  of  the  lateral-stick  frequency-sweep  of 

6l6l 

Fig.  4.2.  For  illustrative  simplicity,  the  process  noise  spectrum  was 
also  obtained  by  passing  a  Gaussian  random  sequence  through  the 
1.0  rad/sec  low-pass  filter.  This  yields  a  signal-to-noise  ratio  S(u>) 
which  is  independent  of  frequency,  and  is  Just  equal  to  the  stick  input- 
to-noise  rms  ratio. 

A  simulation  program  was  developed  to  numerically  integrate  the 
governing  equations  of  Fig.  B.1  and  provide  time-histories  for 
frequency-response  identification  using  the  FRESPID  program.  Results 
were  obtained  as  a  function  of  the  noise-to-signal  ratio  o  =  1/S  for 
the  nominal  compensation  gain  of  the  XV- 15  (K  =  3.0)  and  a  gain  typical 
of  a  high-bandwidth  control  system  (K  =  6.5  rad/sec). 

B.2.2  Results  of  Numerical  Study 

The  open-  and  closed-loop  frequency-responses  (p/6a  and  p/6L, 
respectively)  for  the  nominal  gain  (K  =  3.0)  are  shown  in  Figs.  B.2 
and  B.4,  along  with  the  ideal  transfer-functions.  The  associated  coher¬ 
ence  functions  shown  in  Figs.  B.3  and  B.5  indicate  nearly  perfect  iden¬ 
tification  over  the  entire  frequency  range,  with  a  slight  random  error 
for  the  open- loop  response  in  the  frequency  range  near  the  dominant 
open-loop  mode  (w  =  0.5  rad/sec).  The  open-loop  frequency-response  is 
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Fig.  B.2.  Identification  of  Open-Loop  Frequency  Response  for  the 
Nominal-Gain  System,  K  =  3.0.  (a)  Magnitude;  (b)  Phase. 


FREQUENCY,  rad/sec 

Fig.  B.3.  Coherence  Function  for  Open-Loop  Roll-Rate 
Identification,  Nominal-Gain  System,  K  =  3.0. 
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Fig.  B.4.  Identification  of  the  Closed-Loop  Frequency -Response  for 
the  Nominal-Gain  System,  K  =  3.0.  (a)  Magnitude;  (b)  Phase. 
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Fig.  B.5.  Coherence  Function  for  Closed-Loop  Response 
Identification,  Nominal-Gain  System,  K  =  3.0. 


also  identified  with  almost  no  error  for  the  high-gain  case  (Figs.  B.6, 
B.7);  the  higher-bandwidth  closed-loop  frequency-response  matches  the 
true  response  (Figs.  B.8,  B.9). 


FREQUENCY,  rad/sec 


Fig,  B.6.  Identification  of  the  Open-Loop  Frequency-Response  for 
the  High-Gain  System,  K  =  6.5.  (a)  Magnitude;  (b)  Phase. 


FREQUENCY,  rad/sec 

Fig.  B.7.  Coherence  Function  for  Open-Loop  Response 
Identification,  High-Gain  System,  K  =  6.5. 
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FREQUENCY,  rad/sec 

Fig.  B.8.  Identification  of  the  Closed-Loop  Frequency-Response  for 
the  High-Gain  System,  K  =  6.5.  (a)  Magnitude;  (b)  Phase. 


FREQUENCY,  rad/sec 

Fig.  B.9.  Coherence  Function  for  Closed-Loop  Response 
Identification,  High-Gain  System,  K  =  6.5. 


The  effect  of  increasing  the  noise-to-signal  ratio  a  on  the 
open-loop  identification  for  the  nominal  gain  system  is  shown  in 
Figs.  B.10  and  B . 1 1 .  Very  little  error  is  apparent  for  noise-to-signal 
ratios  of  o  <  0.3.  At  moderate  noise-to-signal  ratios  (o  =  0. 3-1.0), 
the  identification  errors  (associated  with  both  bias  and  random  error 


FREQUENCY,  rad/sec 


Fig.  B.10.*  Effect  of  Noise-to-Signal  Ratios  o  =  0,  0.1,  0.3  on 
the  Identification  of  the  Open-Loop  Frequency-Response  for  the  Nominal- 
Gain  System,  K  =  3.0.  (a)  Magnitude;  (b)  Phase;  (c)  Coherence. 


FREQUENCY,  rad/sec 


Fig.  B.11.  Effect  of  Noise-to-Signal  Ratios  o  =  0,  1.0,  10.0  on 
the  Identification  of  the  Open-Loop  Frequency-Response  for  the  Nominal- 
Gain  System,  K  =  3.0.  (a)  Magnitude;  (b)  Phase;  (c)  Coherence. 
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contributions)  are  significant  and  are  reflected  by  the  drop  in  the 
coherence  function.  Thus,  for  low  and  moderate  noise  levels,  the  coher¬ 
ence  function  remains  a  good  indicator  of  spectral  accuracy  for  closed- 
loop  identification.  The  case  of  nearly  infinite  noise-to-signal  ratio 
(a  -  10)  shows  that  the  open-loop  frequency-response  identification 
approaches  the  inverse  compensation  transfer-function  as  indicated  in 
Eq.  (B. 10) : 


G_ 


-  £  =  0.333  =  -9.54  dB 


(-  g;)  *  4  if)  *  -’80 


(B.20) 


and  the  associated  coherence  function  approaches  unity  (Fig.  B.llc). 

A  better  quantitative  appreciation  for  the  effect  of  process  noise 
on  the  identified  open-loop  frequency-response  is  achieved  by  looking  at 
two  measures  of  error.  The  first  measure  is  the  total  rms  error.  This 
relates  the  deviation  between  the  estimated  open-loop  frequency- 
response  H  and  the  true  open-loop  frequency -response  H  [Eq.  (B. 17) ] : 


'  20  a  ?11/2 

20  [HdB(u)  '  HdB(u)]  J 

’  20  V/2 

so  £  'Wm} '  W“,]2J 


(B.21 ) 


These  calculations  are  performed  by  the  NAVFIT  program,  which 


selects  20  linearly  spaced  frequency  points  across  the  logarithmic  scale 
(Sec.  3.3.2).  The  total  rms  magnitude  and  phase  errors  are  plotted  as  a 
function  of  noise-to-signal  ratio  5  in  Figs.  B.12a  and  B.12b  for  the 
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NOMINAL  XV-15  GAIN,  K  =  3 
HIGH  GAIN,  K  =  6.5 


Fig.  B.12.  Variation  in  the  Total  rms  Errors  with  Noise-to-Signal 
Ratio  for  the  Identification  of  the  Open-Loop  Frequency-Response. 

(a)  Magnitude  Error;  (b)  Phase  Error. 


nominal  and  high-gain  cases.  Based  on  an  allowable  magnitude  error  of 

10H  (o  =0.92  dB),  the  results  of  Fig.  B.12a  indicate  that  the  noise- 

g 

to-signal  ratio  should  be  0  <  0.11  for  the  low-gain  case  and 


3  <  0.20  for  the  high-gain  case.  Adopting  the  relative  error  weighting 


of  1  dB  magnitude  error  to  7  deg  phase  error  (Sec.  3.3.2)  yields  an 

allowable  rms  phase  error  o  =6.4  deg  (=  7  *  0.92  dB).  Based  on  this 

P 

phase  error  limit,  the  results  of  Fig.  B.12b  indicate  an  acceptable 
noise-to-signal  ratio  of  a  <  0.22  for  the  low-gain  case  and 
a  <  0.25  for  the  high-gain  case.  Thus,  the  noise  limits  based  on  the 
magnitude  criteria  are  more  restrictive  for  this  study.  The  larger 
allowable  levels  of  noise  for  the  high-gain  case  occur  because  the 
resulting  closed-loop  system  is  more  resistant  to  excitation  from  pro¬ 
cess  noise,  thereby  reducing  the  total  level  of  random  error.  For  large 
noise-to-signal  ratios  5  >>  1,  the  magnitude  errors  asymtotically 
approach  different  levels  for  the  two  feedback  gain  cases  because  the 
associated  inverse  compensation  gain  values  -1/K  are  different.  The 
rms  phase  errors  approach  a  common  value  because  the  open-loop  identifi¬ 
cation  tends  toward  9  =  - 1 80  deg,  independent  of  the  gain  level 
Eq.  (B .20). 

A  direct  measure  of  parameter  bias  error  is  obtained  by  fitting 
each  of  the  frequency-responses  in  Fig.  B.10  (and  similarily  for  the 
high-gain  case)  with  the  first-order  transfer-function  model: 

*  k  p_TS 

H  =  — -  ( B .  22 ) 

s  -  a 

and  tabulating  the  bias  errors  normalized  with  respect  to  the  true  open- 
loop  transfer-function  parameters  of  Eq.  (B.17): 


The  normalized  bias  errors  for  transfer-function  gain  eb[K]  and 
mode  location  eb[a]  are  plotted  in  Figs.  B.13a  and  B.13b.  The  small 
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Fig.  B . 1 3 -  Variation  in  Normalized  Bias  Errors  with  the  Noise-to- 
Signal  Ratio  for  the  Identification  of  the  Open-Loop  Frequency  Response, 
(a)  Normalized  Gain  Error;  (b)  Normalized  Mode  Location  Error. 


bias  error  for  no  process  noise  is  associated  with  the  overall  spectral 
identification  procedure  (Ref.  28).  The  increase  in  bias  error  for  the 


high-gain  case  at  low  process  noise  levels  (o  <  0.1)  probably  results 
from  the  integration  errors  which  become  more  significant  at  the  higher 
closed-loop  natural  frequencies.  Restricting  the  maximum  bias  error 
(gain  and  mode)  to  within  1 0%  indicates  a  requirement  for  noise-to- 
signal  ratios  of  3  <  0.3,  which  is  in  good  agreement  with  the  theoreti¬ 
cal  results  of  Eq.  (B.16).  The  fact  that  the  bias  errors  are  very  small 
for  process  noise  ratios  of  o  <  0.10  indicates  that  the  random  error 
contribution  is  the  significant  factor  in  the  total  rms  error  results  of 
Fig.  B.12a  and  B.12b.  By  limiting  the  noise  ratio  to  a  <  0.10  in 
order  to  maintain  acceptable  levels  of  total  rms  error,  negligible  bias 
errors  are  ensured. 


B.3  Conclusions  of  Closed-Loop  Identification  Study 

Two  major  conclusions  can  be  drawn  from  this  study  of  limitations 
on  closed-loop  identification: 

1.  Identifying  the  open-loop  frequency-response  from  closed-loop 
flight  data  results  in  a  biased  estimate  whenever  process  noise  is 
present.  Common  sources  of  process  noise  which  must  be  considered  are 
turbulence,  extraneous  unaccounted  for  inputs,  and  nonlinearities. 

2.  Process  noise  causes  random  and  bias  errors  in  the  estimates. 
Numerical  studies  based  on  the  simplified  model  of  XV-15  dynamics 
suggest  that  allowable  levels  of  noise  are  determined  more  from  the 
effects  of  random  error  than  from  those  of  bias  error.  Total  error  can 
be  kept  to  an  acceptable  level  when  the  noise-to-signal  rms  ratio  is 

o  <  0.10.  Bias  errors  do  not  become  important  until  process  ncise 
ratios  reach  o  >  0.3.  The  coherence  function  is  a  good  indicator  of 
the  level  of  total  random  error  for  low  and  moderate  process  noise 


levels,  and  can  be  used  to  assess  the  quality  of  the  open-loop 
frequency-response  identification  from  closed-loop  test  data. 

This  relatively  simple  numerical  simulation  provides  a  good  appre¬ 
ciation  for  the  effects  of  process  noise  and  closed-loop  compensation. 
Similar  studies  should  be  conducted  (perhaps  using  more  sophisticated 
models  of  vehicle  response  and  noise  environment)  to  support  future 
closed-loop  identification  experiments  on  rotorcraft. 
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